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Summary 



Non-equilibrium critical phenomena have attracted a lot of research interest in the 
recent decades. Similar to equilibrium critical phenomena, the concept of universality 
remains the major tool to order the great variety of non-equilibrium phase transitions 
systematically. All systems belonging to a given universality class share the same set 
of critical exponents, and certain scaling functions become identical near the critical 
point. It is known that the scaling functions vary more widely between different uni- 
versality classes than the exponents. Thus, universal scaling functions offer a sensitive 
and accurate test for a system's universality class. On the other hand, universal scaling 
functions demonstrate the robustness of a given universality class impressively. Unfor- 
tunately, most studies focus on the determination of the critical exponents, neglecting 
the universal scaling functions. 

In this work a particular class of non-equilibrium critical phenomena is considered, 
the so-called absorbing phase transitions. Absorbing phase transitions arc expected 
to occur in physical, chemical as well as biological systems, and a detailed introduc- 
tion is presented. The universal scaling behavior of two different universality classes is 
analyzed in detail, namely the directed percolation and the Manna universality class. 
Especially, directed percolation is the most common universality class of absorbing 
phase transitions. The presented picture gallery of universal scaling functions includes 
steady state, dynamical as well as finite size scaling functions. In particular, the effect 
of an external field conjugated to the order parameter is investigated. Incorporating 
the conjugated field, it is possible to determine the equation of state, the suscepti- 
bility, and to perform a modified finite-size scaling analysis appropriate for absorbing 
phase transitions. Focusing on these equations, the obtained results can be applied to 
other non-equilibrium continuous phase transitions observed in numerical simulations 
or experiments. Thus, we think that the presented picture gallery of universal scaling 
functions is valuable for future work. 

Additionally to the manifestation of universality classes, universal scaling functions 
are useful in order to check renormalization group results quantitatively. Since the 
renormalization group theory is the basis of our understanding of critical phenomena, 
it is of fundamental interest to examine the accuracy of the obtained results. Due 
to the continuing improvement of computer hardware, accurate numerical data have 
become available, resulting in a fruitful interplay between numerical investigations and 
renormalization group analyzes. 
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Introduction 



1.1 Opening remarks and outline 

One of the most impressive features of continuous phase transitions is the concept of 
universahty, that allows to group the great variety of different critical phenomena into a 
small number of universality classes (see U for a recent review). All systems belonging 
to a given universality class have the same critical exponents, and certain scaling func- 
tions (equation of state, correlation functions, etc.) become identical near the critical 
point. The universality has its origin in the long range character of the fluctuations. 
Close to the transition point, the corresponding correlation length becomes much larger 
than the typical range of molecular interactions. Then the behavior of the cooperative 
phenomena becomes independent of the microscopic details of the considered system. 
The concept of universality is well established for equilibrium phase transitions where 
a unifying theoretical framework exists. In that case, the universal behavior of short- 
range interacting systems depends only on few fundamental parameters, namely the 
dimensionality of space and the symmetry of the order parameter |2|. Classical ex- 
amples of such universal behavior in critical equilibrium systems are for instance the 
coexistence curve of liquid- vapor systems and the equation of state in ferromagnetic 
systems (e.g. 011]). A complete understanding of the critical behavior of a given sys- 
tem would require to calculate the critical exponents and the universal scaling functions 
exactly. In general, this is only possible above the upper critical dimension where mean 
field theories apply. But the universality ensures that even rather crude modeling of 
complicated microscopic behavior provides quantitatively many essential features of the 
critical behavior. Therefore, highly accurate estimates of critical exponents of various 
universality classes are known (see e.g. jSj). 

In contrast to equilibrium critical phenomena, less is known in case of non-equi- 
librium phase transitions. Here, a generalized treatment is not known, lacking an 
analog to the equilibrium free energy. Thus the rich and often surprising variety of 
non-equilibrium phase transitions (see for example [01 El IHl El ^1 ) observed in physical, 
chemical, biological, as well as socioeconomic systems, has to be studied for each system 
individually. But similar to equilibrium systems it is believed that non-equilibrium crit- 
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ical phenomena can be grouped into universality classes and the concept of universality 
plays again a central role in theoretical and numerical analysis. The universality implies 
in turn that oversimplified representations or caricatures of nature provide quantita- 
tively correct results, if the essential features are captured which are responsible for 
non-equilibrium ordering. Partial differential equations as well as interacting lattice 
models are two established approaches to study non-equilibrium systems. In the first 
case a set of partial differential equations is usually constructed on a mean field level by 
directly translating the reaction scheme (e.g. of a chemical reaction) into equations for 
gain and loss of certain quantities. The typically non-linear dynamics is described by 
deterministic equations, and phase transitions are related to bifurcations fS^. Adding 
suitably chosen noise functions, improved results can be obtained within a Langevin 
approach or a Fokker-Planck description (see e.g. ^DE])- In that case, field theoreti- 
cal approaches assisted by renormalization group techniques are successfully applied to 
obtain results beyond the mean field level. On the other hand, microscopic interacting 
particle systems like lattice-gas models or cellular automata provide another in- 
sight into non-equilibrium critical phenomena. Although an exciting development has 
been seen in the last decade leading to a series of exact solutions of interacting particle 
systems (see e.g. |14|'). most models are not accessible to exact mathematical treatment, 
in particular in higher than one dimension. Thus numerical simulations on increasingly 
powerful computers are widely used in order to obtain quantitative results. 

As pointed out, a full classification of the universality classes of non-equilibrium 
phase transitions is still lacking, i.e., neither the universality classes nor their defining 
fundamental parameters are known. Therefore, numerous (mostly phenomenologically 
motivated) classifications schemes are discussed in the literature. These universality 
hypotheses have to be checked model by model. Due to a lack of analytical solutions, 
numerical simulations or renormalization group treatments are used (often in a fruitful 
interplay) to identify a system's critical behavior, i.e., to specify the order parameter, 
predicting the order of the transition, and describing the scaling behavior in the vicin- 
ity of the transition point via critical exponents and scaling functions. Unfortunately, 
most work focuses on the determination of the critical exponents only, neglecting the 
determination of the universal scaling functions. It turns out that checking the univer- 
sality class it is often a more exact test to consider scaling functions rather than the 
values of the critical exponents. While for the latter ones the variations between dif- 
ferent universality classes are often small, the scaling functions may differ significantly. 
Thus the agreement of universal scaling functions provides not only additional but also 
independent and more convincing evidence in favor of the conjecture that the phase 
transitions of two models belong to the same universality class. 

It is the aim of this work to demonstrate the usefulness of universal scaling functions 
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for the analysis of non-equilibrium phase transitions. In order to limit the coverage 
of this article, we do not present an overview of non-equilibrium phase transitions 
as it was done for example in recent review articles [7| |H]. Instead we focus on a 
particular class of non-equilibrium critical phenomena, the so-called absorbing phase 
transitions. These phase transitions arise from a competition of opposing processes, 
usually creation and annihilation processes. The transition point separates an active 
phase and an absorbing phase in which the dynamics is frozen. A systematic analysis 
of universal scaling functions of absorbing phase transitions is presented, including 
static, dynamical, and finite-size scaling measurements. As a result a picture gallery 
of universal scaling functions is presented which allows to identify and to distinguish 
universality classes. 

The outline of this work is as follows: in the remaining part of the introduction 
a transparent formulation of scaling and universality is presented, since both notions 
are central to the understanding of critical phenomena. Therefore, we follow the his- 
toric perspective and survey phenomenologically the basic ideas of equilibrium critical 
phenomena, and discuss the concepts of critical exponents, generalized homogenous 
functions, scaling forms, as well as universal amplitude combinations. A foundation 
for an understanding of scaling and universality is provided by Wilson's renormaliza- 
tion group theory |15l I16j . which is a topic on its own and not presented in this work. 
Instead we focus on the implications on universal scaling and just illustrate the main 
results, e.g. we sketch how the renormalization group theory allows to decide on the rel- 
evant system parameters determining the universality class. The renormalization group 
also provides a tool for computing critical exponents as well as the universal scaling 
functions, and it explains the existence of an upper critical dimension. For a rigorous 
substantiation of scaling and universality the interested reader is referred to established 
textbooks (e.g. HZl UHl HH] ) and review articles [201 HH |12] • While we have attempted 
to present the introduction as general as possible, we use for the sake of concreteness 
the language of ferromagnetic phase transitions. In conclusion, the subsequent sections 
can be used as an introduction to the theory of phase transitions for readers unfamiliar 
with the concepts of scaling and universality. Readers familiar with equilibrium critical 
phenomena and interested in non-equilibrium systems are recommended to skip the 
introduction. 

In chapter |2l we introduce the definitions and notations of absorbing phase transi- 
tions. To be specific, we present a mean field treatment of the contact process providing 
a qualitative but instructive insight. Various simulation techniques are discussed for 
both steady state and dynamical measurements. In particular, the mean field results 
are compared briefiy to those obtained from simulations of the high-dimensional con- 
tact process. The chapters |21 and |1] are devoted to the investigation of two different 
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universality classes: the directed percolation universality class and the so-called Manna 
universality class. According to its robustness and ubiquity the universality class of 
directed percolation (DP) is recognized as the paradigm of the critical behavior of sev- 
eral non-equilibrium systems exhibiting a continuous phase transition from an active 
to an absorbing state. The widespread occurrence of such models describing critical 
phenomena in physics, biology, as well as catalytic chemical reactions is reflected by 
the universality hypothesis of Janssen and Grassberger: Short-range interacting mod- 
els, which exhibit a continuous phase transition into a unique absorbing state, belong 
to the directed percolation universality class, provided they are characterized by a 
one-component order parameter [231 124j . Different universality classes are expected to 
occur e.g. in the presence of additional symmetries or quenched disorder. The univer- 
sality class of directed percolation is well understood. In particular, a field theoretical 
description is well established, and renormalization group treatments provide within 
an e-expansion useful estimates of both the critical exponents and the scaling func- 
tions. But it should be remembered that the renormalization group rarely provides 
exact results, i.e., it is essential to have an independent (usually numerical or exper- 
imental) check. These checks are of fundamental interest since the renormalization 
group theory is the basis of our understanding of critical phenomena. Remarkably, a 
detailed analysis shows that renormalization group approximations reveals more accu- 
rate results for directed percolation than for standard equilibrium models. In contrast 
to directed percolation, less is known in case of the Manna universality class. For 
example, field theoretical approaches using renormalization group techniques run into 
difficulties. Thus a systematic e-expansion is still lacking, and most quantitative results 
are obtained from numerical simulations. The Manna universality class is of particular 
interest since it is related to the concept of self-organized criticality. 

For both universality classes considered several lattice models are investigated. A 
systematic analysis of certain scaling functions below, above as well as at the upper 
critical dimension is presented. At the upper critical dimension, the usual power-law 
behavior is modified by logarithmic corrections. These logarithmic corrections are 
well known for equilibrium critical phenomena, but they have been largely ignored 
for non-equilibrium phase transitions. Due to the considerably high numerical effort, 
sufficiently accurate simulation data have become available recently, triggering further 
renormalization group calculations and vice versa. The results indicate that the scaling 
behavior can be described in terms of universal scaling functions even at the upper 
critical dimension. Independent of the dimension, we investigate steady state scaling 
functions, dynamical scaling functions, as well as finite-size scaling functions. Focusing 
on common scaling functions such as the equation of state our method of analysis can be 
applied to other non-equilibrium continuous phase transitions observed numerically or 
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experimentally. Thus we hope that the presented gallery of universal scaling functions 
will be useful for future work where the scaling behavior of a given system has to 
be identified. Additionally, certain universal amplitude combinations are considered. 
In general, universal amplitude combinations are related to particular values of the 
scaling functions. Systematic approximations of these amplitude combinations are often 
provided by renormalization group treatments. Furthermore, the validity of certain 
scaling laws connecting the critical exponents is tested. 

Crossover phenomena between different universality classes are considered in chap- 
ter 13 Although crossover phenomena are well understood in terms of competing fixed 
points, several aspects are still open and are discussed in the literature. For example, 
the question whether the full crossover region, that spans usually several decades in 
temperature or conjugated field, can be described in terms of universal scaling func- 
tions attracted a lot of research activity. Here, we consider several models belonging 
to the Manna universality class with various interaction range. Our analysis reveals 
that the corresponding crossover from mean field to non-mean field scaling behavior 
can be described in terms of universal scaling functions. This result can be applied to 
continuous phase transitions in general, including equilibrium crossover phenomena. 

Concluding remarks are presented in Chapter|Hl Furthermore, we direct the reader's 
attention to areas where further research is desirable. The work ends with appendices 
containing tables of critical exponents as well as universal amplitude combinations. 

A comment is worth making in order to avoid confusion about the mathematical 
notation of asymptotical scaling behavior. Throughout this work the Landau notation 
[ini is used, i.e., the symbols oc , ~, and O denote that two functions f{x) and 
g{x) are 

fix) 

f{x) oc g{x) <^=^> = c, Vx (proportional) , 

9{x) 



f{x)=g{x) + 0{x^) 



f{x) -g{x) 



< c , V X > Xo (of the order of). 



f fx) 

/(x) ~ ^(x) <^=^> lim ^ ' = 1 (asymptotically equal) 



X^Xc 



In the following, the mathematical limit x ^ x^. corresponds to the physical situation 
that a phase transition is approached (usually x — > or x — > oo). Beyond this stan- 
dard notation the symbol oc is used to denote that the functions are asymptotically 
proportional, i.e., 

f(x) 

f(x)6cg(x) lim — —— = c (asymptotically proportional). 

— c g{x) 
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1.2 Scaling theory 

Phase transitions in equilibrium systems are characterized by a singularity in the free 
energy and its derivatives [171 |2H1 US] • This singularity causes a discontinuous behavior 
of various physical quantities when the transition point is approached. Phenomenolog- 
ically the phase transition is described by an order parameter, having non-zero value 
in the ordered phase and zero value in the disordered phase [SO]. Prototype systems 
for equilibrium phase transitions are simple ferromagnets, superconductors, liquid-gas 
systems, ferroelectrics, as well as systems exhibiting superfluidity. 

In case of ferromagnetic systems the transition point separates the ferromagnetic 
phase with non-zero magnetization from the paramagnetic phase with zero magne- 
tization. The well-known phase diagram is shown in Figure^ Due to the reversal 
symmetry of ferromagnets all transitions occur at zero external field h. The phase di- 
agram exhibits a boundary line along the temperature axis, terminating at the critical 
point . Crossing the boundary for T < the magnetization changes discontinuously, 
i.e. the systems undergoes a first order phase transition. The discontinuity decreases 
if one approaches the critical temperature. At the magnetization is continuous but 
its derivatives are discontinuous. Here the system undergoes a continuous or so-called 
second order phase transition. For T > T^, no singularities of the free energy occur and 
the systems changes continuously from a state of positive magnetization to a state of 
negative magnetization. 

In zero field the high temperature or paramagnetic phase is characterized by a van- 
ishing magnetization. Decreasing the temperature, a phase transition takes place at 
the critical temperature and for T < one observes an ordered phase which is sponta- 
neously magnetized (see Figure^. Following Landau, the magnetization m is the order 
parameter of the ferromagnetic phase transition (HO] • Furthermore the temperature T 
is the control parameter of the phase transition and the external field h is conjugated to 
the order parameter. As well known, critical systems are characterized by power-laws 
sufficiently close to the critical point, e.g. the behavior of the order parameter can be 
described by 



with the reduced temperature 5T = (T — T^)/T^ and the exponent (3. 

For non-zero field the magnetization increases smoothly with decreasing tempera- 
ture. At the critical isotherm {5T = 0) the magnetization obeys another power-law for 





h 







(1.2) 



Further exponents are introduced to describe the singularities of the order parameter 
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Figure 1: Sketch of the ferromagnetic phase diagram (left) and the corresponding order parameter 
behavior (right). The phase diagram contains a line of first-order transitions at zero field (dashed line) 
that ends in a critical point at temperature T^. The temperature dependence of the magnetization is 
sketched for zero-field (thick line) and three different values of the conjugated field. 



x{T,h 


= 0) 


6c i^rr^ 


C{T,h 


= 0) 


6c i^rr" 


mh 


= 0) 


6c i^rr'^ 



susceptibility X; the specific heat C, as well as the correlation length ^ 

(1.3) 
(1.4) 
(1.5) 

Another exponent 77 describes the spatial decay of the correlation function at criticality 

r(r) 6c , (1.6) 

where D denotes the dimensionality of the system. 

A comment about vanishing exponents is worth making. A zero exponent corre- 
sponds either to a jump of the corresponding quantity at the critical point or to a 
logarithmic singularity since 

-ln\6T\ = lim-^ . (1.7) 

S-+0 S 

Often, it is notoriously difficult to distinguish from experimental or numerical data a 
logarithmic singularity from a small positive value of the exponent. 

The exponents a, (3, 7, 5, t], and u are called critical exponents. Notice in Eqs. I\l.'6i 
II .5)1 the equality of the critical exponents below and above the critical point, which is an 
assumption of the scaling theory. The phenomenological scaling theory was developed 
by several authors in the 1960s (e.g. EH ESI 133 ESI ES] ) and has been well verified 
by experiments as well as simulations. In particular, the scaling theory predicts that 
the critical exponents mentioned above are connected by the scaling laws 

a + 2/3-^7 = 2 (Rushbrooke), (1.8) 
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7 = /3(6-l) (Widom), (1.9) 

7 = {2-r])u (Fisher), (1.10) 

2 — a = u D (Josephson) . (l-H) 



The Josephson law includes the spatial dimension D of the system and is often termed 
the hyperscaling relation. In this way the critical behavior of an equilibrium system is 
determined by two independent critical exponents. 

The scaling theory rests on the assumption that close to the critical point the singu- 
lar part of a particular thermodynamic potential is asymptotically a generalized homo- 
geneous function (see e.g. jHZI)- A function /(xi,X2, . . .) is a generalized homogeneous 
function if it satisfies the following equation for all positive values of A 

Xf{x„x,...) = f{X''x„X'H„...). (1.12) 

The exponents Si, S2, . . . are usually termed scaling powers and the variables Xi, X2, . . . 
are termed the scaling fields. In case of ferromagnetism the singular part of the Gibbs 
potential per spin gsing{ST,h) is assumed to scale asymptotically as 

gsU^T,h) ~ A"-2 ~g{6TX,hX''^), (1.13) 

with the scaling function g{x, y). The scaling power of the conjugated field is often de- 
noted as the gap-exponent A = P 6. Strictly speaking, Eq. H1.13() is only asymptotically 
valid, i.e., only when 5T and h tend to zero. Corrections occur away from this limit. It 
turns out that all Legendre transformations of generalized homogenous functions and 
all partial derivatives of generalized homogenous functions are also generalized homoge- 
nous functions. Thus with the Gibbs potential all thermodynamic potentials and all 
thermodynamic functions that are expressible as derivatives of thermodynamic poten- 
tials, like the magnetization, specific heat, etc., are generalized homogeneous functions. 

Consider, for example, the magnetization and the corresponding susceptibility. Dif- 
ferentiating Eq. (|1.13|) with respect to the conjugated field we find 

m{6T,h) = -(^^] ~ A"-2+/35 ^(5rA,/iA^^), (1.14) 

\ oh / 2" 

x{ST,h) = ~ A--2+2/5^ x(5TA,/iA^'^), (1.15) 
where the scaling functions are given by 

-(x,y) = -(^^^j^, X{x,y) = -^-^^j^. (1.16) 

Since the Eqs. (11.14111.15]) are valid for all positive values of A they hold for A = 1/|(5T|, 

hence one obtains at zero-field 

m{6T,0) ~ (_5r)-("-2+/35) rh{-l,0), (1.17) 
X{ST,0) ~ |jr|-(-2+2;35) x(±l,0), (1.18) 
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where the magnetization is defined for 6T < only. Notice that one expects x{+l, 0) ^ 
x(— 1, 0) in general, i.e., the amplitudes of the susceptibility are different below (T < TJ 
and above {T > T^) the transition. Comparing with Eq. and Eq. (|1.3() we find 

f3 = -a + 2-p6, 7 = a-2 + 2/35, (1.19) 

which leads directly to the Rushbrook [Eq. H1.8|) ] and Widom [Eq. ()1.9|) ] scaling law. The 
Fisher and Josephson scaling laws can be obtained in a similar way from the scaling 
form of the correlation function (see e.g. |19j). In order to obtain the Josephson law 
both thermodynamic scaling forms and correlation scaling forms have to be combined. 
Scaling relations obtained in this way are usually termed hyperscaling laws and do not 
hold above the upper-critical dimension. 

The scaling theory implies still more. Consider for instance the M-h-t equation of 
state [Eq. (fTT^ ] . Choosing hX^^ = 1 in Eq. (fTTI)) we find 

m{6T,h) ~ A"^ ^i^'^^^^^^^)\^_i^-ips = ^^'^ m{5Th-^/^\l). (1.20) 

At the critical isotherm {5T = 0) we recover Eq. (|1.2[) . Furthermore, the equation of 
state may be written in the rescaled form 

~ m(5r,.,l) with m„ = m/^"^/^ = 5T h'^/^^ . (1.21) 

In this way the equation of state is described by the single curve m(x, 1) and all M-h-t 
data points will collapse onto the single curve m(x, 1) if one plots the rescaled order 
parameter as a function of the rescaled control parameter 6T^ . This data-collapsing 
is shown in FigureE] for the magnetization curves of Figure^ 

A different scaling form is obtained if one considers instead of the Gibbs potential 
the Helmholtz potential a = g + mh. Since Legendre transforms of generalized homo- 
geneous functions are also generalized homogeneous functions the singular part of the 
Helmholtz potential obeys 

a,i„g(<5r,m) ~ A""2 d{6TX,mX'^). (1.22) 

This equation leads to the scaling form of the conjugated field 

hi6T,m) = (^^) ^ X-l^^ h{5TX,mX^). (1.23) 
\ dm J rjp 

Choosing m A'^ = 1 we find 

/i™ ~ /i((5T„.,l) with h^ = hm^\ 5T^ = 5T m'^l'^ . (1.24) 

Both, h{x, 1) and m(x, 1) are analytic in the neighborhood of x = 0, i.e., at the critical 
temperature. The function /i(x, 1) is often called the Widom-Griffiths scaling func- 
tion [331 136j whereas we term in the following 77i(x, 1) as the Hankey-Stanley scaling 
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-10 



-5 C 
a^6r h 



-1/ps 



5 



10 



2 



Figure 2: The data-collapsing of the M-h-t equation of state. The data correspond to the three different 
magnetization curves of Figure^ All data points collapse onto the curves rh{x, 1) (left, Hankey-Stanley 
scaling form [Eg. I|1.21^ ]'l and h(x,l) (right, Widom-GrifBths scaling form [Eg. I|1.21|l ]'). respectively. 
The so-called metric factors are defined in section IT31 The circles mark the normalization conditions 
Eq. Jon and Eg. iP^ . 

function [HZI- The corresponding data-collapses for both scahng forms are presented 
in FigureEJ The Hankey-Stanley scaling form is just the order parameter curve as a 
function of the control parameter in a fixed conjugated field. It is therefore the natural 
and perhaps more elegant way to present data-collapses of the equation of state. But 
often the mathematical forms of the Hankey-Stanley functions are rather complicated 
whereas the Widom- Griffiths scaling forms are analytically tractable. Therefore h{x, 1) 
is often calculated within certain approximation schemes, e.g. e- or 1/n-expansions 
within a renormalization group approach. 

Notice that the definition of a generalized homogeneous function [Eq. H1.12() ] and 
the data-collapse representations 



are mathematically equivalent, i.e., a function f{xi,X2) is a generalized homogeneous 
function if and only if Eq. H1.25j) . or equally Eq. 1)1. 26|) . is fulfilled ^7\. This was impor- 
tant for the phenomenological formulation of the scaling theory since it has been borne 
and has been confirmed by numerous data-collapses of experimental and numerical 
data. 






(1.26) 



1.3 Universality 



According to the scaling laws Eqs. ()1. 8111.111) equilibrium phase transitions are charac- 
terized by two independent critical exponents. In the 1950s and 1960s it was experi- 
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mentally recognized that quantities like r^, depend sensitively on the details of the in- 
teractions whereas the critical exponents are universal, i.e., they depend only on a small 
number of general features. This led to the concept of universality which was first clearly 
phrased by Kadanoff [2], but based on earlier works including e.g. [HH l IHI H HH I HT j H^. 
The universality hypothesis reduces the great variety of critical phenomena to a small 
number of equivalence classes, so-called universality classes, which depend only on few 
fundamental parameters. All systems belonging to a given universality class have the 
same critical exponents and the corresponding scaling functions become identical near 
the critical point. For short range interacting equilibrium systems the fundamental 
parameters determining the universality class are the symmetry of the order parameter 
and the dimensionality of space pi I41j. The specific nature of the transition, i.e. the 
details of the interactions, such as the lattice structure and the range of interactions (as 
long as it remains finite) do not affect the scaling behavior. For instance, ferromagnetic 
systems with one axis of easy magnetization are characterized by a one component order 
parameter (n = 1) and belong to the universality class of the Ising ferromagnet. Fur- 
thermore, liquid-gas transitions |^ El 1^ El EH] > binary mixtures of liquids jlH], 
and systems exhibiting an order-disorder transition in alloys such as beta-brass are 
described by a scalar order parameter, too, and belong therefore to the same univer- 
sality class (see e.g. jSOl)- Even phase transitions occurring in high-energy physics are 
expected to belong to the Ising universality class. For example, within the electroweak 
theory the early universe exhibits a transition from a symmetric (high temperature) 
phase to a broken so-called Higgs phase jSJ. The occurring line of first order phase 
transitions terminates at a second order point which is argued to belong to the Ising 
class. 

Ferromagnetic systems with a plane of easy magnetization are characterized by a 
two-component order parameter (n = 2). Representatives of this universality class 
are the XY-model, superconductors, as well as liquid crystals that undergo a phase 
transition from a smectic-A to a nematic phase |52| I54j . But the most impressive 
prototype is the superfluid transition of ^He along the A-line. Due to its characteristic 
features like the purity of the samples, the weakness of the singularity in the compress- 
ibility of the fluid, as well as the reasonably short thermal relaxation times, superfluid 
^He is more suitable for high-precision experiments than any other system |55| . For 
example, orbital heat capacity measurements of liquid helium to within 2 nK from the A- 
transition provide the estimate a = —0.01285 it 0.00038 j56ll57j . Thus, the A-transition 
of "^He offers an exceptional opportunity to test theoretical predictions, obtained in 
particular from renormalization group theory (see e.g. jSH] and references therein). 

The well known Heisenberg universality class describes isotropic ferromagnetic sys- 
tems that are characterized by a three component order parameter (n = 3). Despite of 
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the Ising, XY, and Heisenberg universality classes other classes with n > 4 are discussed 
in the literature. For instance, the n = 5 universality class is expected to be relevant 
for the description of high-Te superconductors |591 l6Uj whereas n = 18 is reported in 
the context of superfluid ^He |Hlj . Furthermore, the limit n ^ corresponds to the 
critical behavior of polymers and self-avoiding random walks |62j . The other limiting 
case n — > oo corresponds to the exactly solvable spherical model 1641 165j . 

In the scaling theory presented above, the exponents are already universal but the 
scaling functions are non-universal. Thus two systems are characterized by e.g. two 
different scaling functions rh although they are in the same universality class. The 
non-universal features can be absorbed into two non-universal parameters leading to 
the universal scaling function M 

c^m{5T,h) ~ A"^ M{cT6TX,c,,hX'^^). (1.27) 

In this way, the universal scaling function M is the same for all models belonging 
to a given universality class whereas all non-universal system-dependent features are 
contained in the metric factors c^, c^, and |66j . Since this scaling form is valid 
for all positive values of A, the number of metric factors can be reduced by a simple 
transformation. For the sake of convenience we choose c^^^A i-^ A, yielding 

m{6T,h) ~ A"'^ M{aT6TX,a,,hX^^), (1.28) 

with and a,j = c^c^ , respectively. The universal scaling function M is 

normed by the conditions 

M(-1,0) = 1, M(0,1) = 1, (1.29) 

and the non-universal metric factors are determined by the amplitudes of 

m{5T,h = Q) - A"'^ M(aT(5rA,0) = {-a^STf , (1.30) 

ax&T A=— 1 

m((5r = 0,/i) ~ A-'^ M(0,a„/iA^^) = {a^h)^'^ . (1.31) 

Analogously, the universal Widom- Griffiths scaling form is given by 

ahh{5T,m) ~ A^^^^ ^(aT^r A, m A^) . (1.32) 

The universal scaling function H is usually normed by the conditions 

^(-1,1) = 0, F(0,1) = 1, (1.33) 

which correspond to Eq. (|1.29|) . i.e., the metric factors are again determined by the 
amplitudes of the power-laws m{5T, = 0) ~ (— Ot 

6Tf and m{6T = 0,h) {a^hy/^ 

respectively. 
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For example, let us consider a mean field theory of a simple ferromagnet. Following 
Landau, the free-energy is given by 30 

F{6T,h,m) - Fo = ^biSTm'^ + ^hm^ - hm, (1.34) 

where the positive factors 62 and 64 are system dependent non-universal parameters. 
Variation of the free energy with respect to the magnetization yields the equation of 
state 

h = b^STm + Km^ . (1.35) 

At zero-field we find 



m = or m = \ —-r- 5T P = 7: : = > (1.36) 



whereas at the critical isotherm 



m=^ift ^ J = 3, a,. = i, (L37) 

From Eq. (|1.35)) it follows directly 

H{x,y) = xy + (1.38) 

and in particular the universal Widom-Griffiths scaling form 

a,,hm-^ = HiarSTm-^/^,!) with H{x,l) = 1 + x. (1.39) 

On the other hand the form h{x, 1) =64-1- b2X is obviously non-universal. Since the 
magnetization is a cube root function (6 = 3) the universal Hankey-Stanley scaling 
form is more complicated than the corresponding Widom-Griffiths form 

1/3 



^^-'y>'^ = - 7Z — ...... _,xi/3 - (i-^o) 



27y + y/W8x^ + 7292/2 ) 

(27y + V108x3T729r 

The data presented in Figure^ and Figure|2l correspond to the mean field treatment 
presented above. In addition to the universal scaling functions M{x, 1) and H{x, 1) the 
normalizations Eq. (|1.29j) and Eq. I|1.33|) are marked in Figure|2l 

Celebrated examples of universal scaling plots, providing striking experimental ev- 
idences for the concept of universality at all, are shown in FigureOland Figure|lJ First, 
Figure|31 presents the coexistence curves of eight different fluids characterized by differ- 
ent interatomic interactions. This plot was published by Guggenheim in 1945 and is 
one of the oldest scaling plots in history 0. The corresponding Widom-Griffiths scaling 
function is also shown in FigureEl Here, the rescaled data of the chemical potential of 
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Figure 3: The universal scaling plot of the coexistence curve of a eight different fluids, undergoing a 
liquid-gas phase transition (left). The solid line corresponds to a fit assuming the exponent /? — 1/3. 
The figure is taken from E.A. Guggenheim, J. Chem. Phys. 13, 253 (1945). The right figure displays 
the corresponding universal Widom-Griffiths scaling function H{x, 1). The scaling variable is defined 
as a; = ATlApl-i/" and Xq is related to the amplitude B of the power-law for the coexistence curve 
Ap = bat". The fi gure is taken from J. V. Sengers and J.M.H. Levelt Sengers, Critical phenomena 
in classical fluids in Progress in liquid physics, edited by C. A. Croxton, (John Wiley&Sons, New York, 
1978), 

five different gases collapse onto a single universal curve. This figure was published by 
Sengers and Levelt Sengers |17| . 

Second, Figure^] displays the universal scaling function M{x, 1) of the three-dimen- 
sional Heisenberg universality class. Here the data of five different magnetic materials 
(note that none of these materials is an idealized ferromagnet PP) agree with the cal- 
culated curve of the Heisenberg model, obtained from a series expansion j(i7| Both 
universal scaling plots demonstrate perfectly the robustness of each universality class 
with respect to variations of the microscopic interactions. 

In addition to the critical exponents and scaling functions one often refers the 
universality also to certain amplitude combinations (see [BHl for an excellent review). 
These amplitude combinations are very useful in order to identify the universality class 
of a phase transition since the amplitude combinations vary more widely than the cor- 
responding exponents. Furthermore, the measurement of amplitude combinations in 
experiments or numerical simulations yields a reliable test for theoretical predictions, 
obtained e.g. from renormalization group approximations. Consider the singular be- 
havior of the susceptibility by approaching the transition point from above and below 



x{T>T^,h = 0) ~ a^,+ 6T-\ 



(1.41) 
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Figure 4: The universal scaling plot of the equation of state of the three-dimensional Heisenberg uni- 
versality class {n — 3). The experimental data of five different magnetic materials are shown: CrBra 
(characterized by a lattice anisotropy), EuO (second neighbor interactions), Ni (itinerant-electron fer- 
romagnet), YIG (Yttrium Iron Garnet, ferrimagnet), PdsFe (ferromagnetic alloy). Despite of the 
different interaction details all data collapse onto the universal scaling function. Furthermore the solid 
line was obtained from a series expansion of the Heisenberg model (see [H I67p. The figure is taken 
from H.E. Stanley, Rev. Mod. Phys. 71, S358 (i999). 



x{T<T^,h = 0) ~ a^^_{-6T)-^. (1.42) 

The amplitudes a^_+ and a^^_ are related to the non-universal metric factors and to 
particular values of universal scaling functions. The universal scaling form of the sus- 
ceptibility is obtained from Eq. (|1.28|) 

x{5T,h) = f^JJ^^lS.) ^ a,X^ X{arSTX,a,h\^'), (1.43) 



dh 



T 



with X(a;,y) = dyM{x,y). Setting a^lf^rlA = 1 we obtain for the amplitudes 

a^_+ = a„ar^X(+l,0) (1.44) 

_ = ai^Ur'^ X{-1,0) . (1.45) 

Obviously, the amplitudes a^_+ and a^__ are non-universal, but the ratio 

a^,+ X(+1,0) 



X(-1,0) 



:i.46) 



is a universal quantity. For example, the mean field behavior of the susceptibility takes 
the form 
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leading to [see Eq. (|1.36|) ] 

Similar to the amplitude ratio of the susceptibility other universal combinations can 
be defined. Well known and experimentally important is the quantity (see ^]) 

= TD,B^-^. (1.49) 

Here, T, D^, B are the traditional, but unfortunately unsystematical, notations for the 
amplitudes of 



Xr^T6T-^ , h^D.m^ , m ^ B {-STf . (1.50) 



T>Tc 



T = Tr 



These amplitudes correspond to the values [see Eq. H1.28() ] 

r = a^,+ = aha5:^X(l,0) , A = 0^7^/(0, 1)^^ B = a^M{-l,0). (1.51) 

Using the normalizations M(— 1, 0) = M(0, 1) = 1 we find for the amplitude combina- 
tion 

R^ = X(1,0) (1.52) 

which is obviously a universal quantity. In case of the mean field theory we find R^ = 1. 
These two examples show how the universality of amplitude combinations emerges nat- 
urally from the universality of the scaling functions, i.e., universal amplitude combina- 
tions are just particular values of the universal scaling functions. 

The above presented phenomenological concepts of scaling and universality have 
been tested in a variety of systems with remarkable success. Nevertheless, they have 
certain shortcomings. For example, there is no way of determining the critical exponents 
and scaling functions explicitly. Furthermore, no mathematical substantiation for the 
underlying scaling form of thermodynamic potentials is provided. It requires Wilson 
renormalization group theory to remedy these shortcomings. This is sketched in the 
next section. 



1.4 Remarks on renormalization group theory 

A foundation for the understanding of the scaling theory and the concept of univer- 
sality has been provided by Wilson's renormalization group (RG) theory |151 116j . In 
equilibrium systems the RG theory maps the critical point onto a fixed point of a 
certain transformation of the system's Hamiltonian (introductions are presented in 
e.g. [22 El El ^1)- In case of the instructive real space RG the transformation 
contains a rescaling of a microscopic length scale, e.g. the lattice spacing a, by a factor 
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b (a I— > 6 a) and the ehmination of those degrees of freedom that correspond to the 
range between a and ba. This rescahng will change the system's properties away from 
the critical point where the system exhibits only finite characteristic length scales. But 
at criticality the infinite correlation length [Eq. (|1.5j) ] determines the physical behavior 
and the properties of the system remain unaffected by the rescaling procedure. In this 
way, criticality corresponds to a fixed point of the renormalization transformation. 

Denoting the system's Hamiltonian by 7i and the rescaled Hamiltonian by 7i' the 
renormalization transformation is described by an appropriate operator TZ, 

n' = -Ron. (1.53) 

Fixed point Hamiltonians 7i* satisfy 

n* = -Ron* (1.54) 

and it turns out that different fixed point Hamiltonians are related to different univer- 
sality classes |21| . For the sake of concreteness consider the reduced Ising Hamiltonian 

n = = ^E^^^. + ^E^" (1-55) 

with the spin variables S = ±1, the nearest neighbor interaction coupling J, and the 
homogeneous external field h, and where the first sum is taken over all pairs of neigh- 
boring spins on a given D-dimensional lattice. The partition function for a system of 
spins is 

{S} 

where the sum is taken over all possible spin configurations and where we introduce 
the couplings 

^1 = t4^, K., = (1.57) 
In general the Hamiltonian can be written as a function of the couplings 

nm = (1.58) 

i 

where Oi are the operators appearing in the Hamiltonian. Performing a renormalization 
transformation reduces the spin numbers {N ^ N' = N/b^) and leads to a rescaled 
Hamiltonian, characterized by the couplings 

K[ = K[{K„K,), Ki = Ki{K„K,), K', = K',{K„K,) , . . . . (1.59) 

Here, K[^^ account for additional coupling terms of the renormalized Hamiltonian which 
may appear as a result of the renormalization transformation even if they are not present 
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in the initial Hamiltonian. These so-cahcd RG recursion relations generate trajectories 
in the space of couplings, i.e., the couplings K_ flow under successive renormalizations 

K — ^ K^^^ — ^ K^^^ — ^ ... — ^ — ^ ^("+1) — ^ ... (1.60) 
n n n n n n 

along the RG trajectories towards a certain fixed point K_*. If the system is not initially 
at criticality the couplings will flow towards a trivial fixed point, e.g., a fixed points 
that corresponds to zero or infinite temperature. 

Linearizing the problem close to a fixed point yields 

non{K* + SK) = n* + noJ2 SKiOi 

i 

i 3 

= ii* + (1.61) 

3 

Assuming that the diagonalized operator C has the eigenoperators and eigenval- 
ues Aj such that Cili = AiO^ wc find that the couplings transform in the diagonal 
representation (Sk) according to 5k[ = AjJ^j, thus 

7^ o ^{K* + SK) = n* + ^i^i ■ (1-62) 

i 

The couplings SKi are termed scaling fields and their recursion relation can be expressed 
by the rescaling factor b 

5k', = b^^ 5Ki with Ai = b^^ . (1.63) 

Successive renormalizations correspond to 

— > b^' 5Ki — > b^y^ — > b^y^ 5k, — > ... . (1.64) 
n n n n 

Thus, the renormalization flow depends in the vicinity of a given fixed point H* on the 
exponents y,. For j/j > (A^ > 1) the corresponding scaling field 5Ki is called relevant 
since successive renormalization transformations drive the system away from Ti.*. In 
case of a negative exponent < (A^ < 1) the system approaches the fixed point 
under repeated transformations and the scaling field is termed irrelevant. Marginal 
scaling fields correspond to = (A^ = 1) and require higher than linear order in the 
expansion. In this way each fixed point is characterized by its associated scaling fields 
and by a domain of attraction which corresponds to the set of points fiowing eventually 
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to the fixed point. This set forms a hypersurface in the space of couplings and is termed 
the critical surface. 

In summary, a fixed point 7i* is approached if all associated relevant scaling fields 
are zero, otherwise the system flows away from 7i*. Examples for relevant scaling fields 
in ferromagnetism are the temperature 6T and the conjugated field h. Criticality is 
only achieved for 5T ^ and h ^ 0, therefore we may identify 

6ki = 6kt = UtST , yi = y-r > , (1.65) 
6k2 = 5ku = a^h, y2 = Vh > , (1.66) 

and yi>2 < 0. On the other hand, all Hamiltonians that differ from the fixed point 
7i* only by irrelevant scaling fields fiow towards TC*. For example the five magnetic 
materials presented in Figure|3 differ only by irrelevant scaling fields. 

Although the above linear recursion relations [Eq. (|l.(i:-{|) ] describe the RG trajec- 
tories only in the vicinities of fixed points they provide some insight in the topology 
of the entire RG fiow (see FigureEl). These RG fiow diagrams are useful in order to 
illustrate the RG transformations schematically and present a classification scheme in 
terms of fixed point stability. The stability of a given fixed point is determined by 
the number of relevant and irrelevant scaling fields. Instable fixed points are char- 
acterized by at least one relevant scaling field since Hamiltonians arbitrarily close to 
the fixed point will flow away under successive RG iterations. Ordinary critical points 
correspond to singly unstable fixed points, i.e., unstable with respect to the control pa- 
rameter (e.g. temperature) of the phase transition. Tricritical points are characterized 
by a second instability. An applied external field conjugated to the order parameter 
implies an additional instability of the fixed point. 

Furthermore, the stability of fixed points depends on the spatial dimensionality D 
of a system. Above a certain dimension the scaling behavior is determined by a 
trivial fixed point with classical mean field exponents, whereas a different fixed point 
with non-classical exponents determines the scaling behavior below D^. This exchange 
of the scaling behavior is caused by an exchange of the stability of the corresponding 
fixed points below and above j^|^[73j. At the so-called upper critical dimen- 
sion both fixed points are identical and marginally stable. The scaling behavior 
is characterized by mean field exponents modified by logarithmic corrections |74| I2()j. 
We will discuss the scaling behavior of certain non-equilibrium phase transitions at the 
upper critical dimension in detail in the following chapters. 

Let us now consider how scaling emerges from the renormalization transformation. 
It is essential that the partition function is invariant under the renormalization opera- 
tion 7^ m] 
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Figure 5: Sketch of the renormahzation group flow of an Ising ferromagnet (see [Eq. (II. 551 1*) on a 
hierarchical lattice. Hierarchical lattices are iteratively constructed lattices |69l I70| where the so-called 
Migdal-Kadanoff scheme of renormahzation is exact |71l I72| . The flowlines show the motion of the 
coupling constants x = exp{2Ki) (temperature-like variable) and y — exp{K2) (field- like variable) 
under successive iterations of the recursion relations x' — (x^ -I- y'^){x~'^ + j/^)(l -I- y^)~^ and y' = 
3/^(1 -I- x^y^){x^ + y^)~^ . The flowlines are attracted and repelled by the fixed points (1, 1), (l,oo), 
(oo, 1), and {Xa,l) with x^ = 3.38298.... The phase transition corresponds to the non-trivial fixed 
point (Xc, 1) with the eigenvalues Ai ~ 1.6785 and A2 ~ 3.6785, leading to /3 « 0.16173 and S ~ 15.549. 



Therefore, the free energy per degree of freedom transforms according to 

f{6KT,6Kh,SKs,...) = b^^ f{5K'^,6K[,6K'^, . . .) . (1.68) 

Combining this equation with Eq. (ILfjl-^j) and using the identifications Eq. (|l.(i5j) and 
Eq. (|1.66|) yields the scaling form 

f{ar5T,a,h,5K,,...) ~ b-'^ f{by^ a^ST^by^ a^h,hy^ 5k,, . . .) . (1.69) 
Introducing A = 6^^ we obtain the scahng form of the free energy [see Eq. H1.13() ] 

f{ar5T,a^,h,6K,,...) ~ A"-^ f{ar5T\a^hX^\5K,X^\...), (1.70) 
where we have identified the exponents 

Vt = - , Vh = — , y,>2 = • (1-71) 

V V V 

The possible additional scaling fields 5Kj>2 deserve comments. Irrelevant scaling fields 
('/'i>2 < 0) may cause corrections to the asymptotic scaling behavior For example, 
choosing A = l/|aT(5T| we obtain at zero field 

/(aT(5r,0,5K3) ~ |aT(5r|2-° /(±1, 0,5^3 |aT(5r|l'^3l) (1.72) 



.<5r|2-" /(±1, 0, 0) + la^^Tp-^+l-^al q, : 



... . 

a; = 
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The non-universal corrections to the leading order |(5Tp~" are termed confluent singu- 
larities and they determine the size of the critical region. Often confluent singularities 
have to be taken into account in order to analyze high precision data. Impressive exam- 
ples of confluent singularity effects of super fluid Helium are reviewed in [HSl- The above 
expansion of /(±1,0, x) implies that the free energy is an analytic function in 6K.3. If 
the free energy is non-analytic the scaling field Sk^ is termed a dangerous irrelevant 
variable |76l I77j . In that case, the free energy exhibits e.g. a power-law divergence 

f{x,y,z) = z-''f{x,y) (1.73) 

characterized by the exponent /x > 0. Singularities of this type occur for example in the 
mean field regime of the well-known Landau-Ginzburg- Wilson Hamiltonian for short 
range interacting ferromagnets (see e.g. j2E])- There, the dangerous irrelevant variable 
corresponds to the coupling constant of the (j)^ interactions. The nonanalytic behavior 
leads to the modified scaling form of the free energy 

= X-''^-^"f'^ Sk-^" fiarST A, a„/i A^^) 

= \aT6T\''^+f"''' (5/^3-^ /(±1, 0) . (1.74) 

Compared to the standard behavior / 6c |5Tp~" the above result refiects the breakdown 
of the hyperscaling law 2 — a = I'D. Additionally to the violation of scaling laws, 
dangerous irrelevant variables also cause the breakdown of common finite-size scaling 
within the mean field regime |79| I77j . This is well established in equilibrium e.g. in the 
n ^ 00 limit. We will address this point in detail in Chapter El where high-dimensional 
non-equilibrium phase transitions are considered. 

The situation is different when the scaling field 5k^ is relevant, i.e., (j)3 > 0. In that 
case the free energy at zero field is given by 

f{aT6T,0,5K3) ~ |aT<5r|2-°/(±l,0,(5K3|aT<5r|-'^3). (1.75) 

For sufficient small arguments (|(5ft;3|(5T|~'^^ | <^ 1) the relevant scaling field leads again 
to corrections to the asymptotic scaling behavior. But approaching the transition point 
(5T — > 0) the scaling argument diverges and gives rise to a different critical behavior, 
i.e., the system crosses over to a different universality class. We will discuss crossover 
phenomena in detail in ChapterEJ Eventually a marginal scaling field causes logarithmic 
corrections via 

(^^3 lor^rl"'^^ = 6K3 exp {—(p3 ln\aT5T\) 

— > 5k3 (l-(/)3 ln|aT(5r| ...). (1.76) 

03^0 
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Often, these logarithmic contributions mask the power law singularities and make the 
analysis of experimental or numerical data notoriously difficult. 

Analogous to the free energy, the renormalization group implies the scaling form of 
the correlation length ^. Performing a renormalization transformation, the correlation 
length ^, like all length scales, is decreased by the factor 6, 

e' = b-'C. (1.77) 

It is essential for the understanding of phase transitions that fixed points are char- 
acterized by an infinite (or trivial zero) correlation length since ^ satisfies at a fixed 
point 

e' = c. (1.78) 

In this way, a singular correlation length, or in other words, scale invariance is the 
hallmark of criticality. The scaling form of the correlation length is obtained from 
Eq.dEZZl) 

^{6KT,6Kh,SK3,...) = b^i6K'j,,6K[,6K'^,...), (1.79) 
yielding with A = and Eq. (fT^ 

^{aTdT,a,,h,6Ks,...) ~ X'" ^ardT X^a^hX^^dK-.X'^^ , . . .) . (1.80) 

Instead of the real space renormalization considered so far, it is convenient to work 
in momentum space. This can be achieved by reformulating the above derivations in 
terms of Fourier transforms. We refer the interested reader to the reviews of |16[ I38 | [T7]. 
The momentum space formulation allows a perturbative definition of the RG which 
leads technically to a formulation in terms of Feyman graph expansions. The appro- 
priate small parameter for the perturbation expansion is the dimensionality difference 
to the upper critical dimension e = — D [HOI, i-e., the e-expansion gives systematic 
corrections to mean field theory in powers of e. The e-expansion provides a powerful 
tool for calculating the critical exponents and the scaling functions. For example, the 
exponent 7 for n-component magnetic systems with short range interactions is given 
in second order by (see e.g. |73l IHT]) 

n + 2 2 + 2) 

2(n + 8) ^ ^ 4(n + 

Furthermore the Widom- Griffiths scaling function can be written as a power series in e 

H{x,l) = 1 + X + eH^{x,l) + e"^ H^ix,!) + 0{e^) . (1.82) 

Obviously the mean field scaling behavior [Eq. (|l..S9j) ] is obtained for e = 0. The 
functions become complicated with increasing order and we just note jS2] 

H,{x,l) = [3(x + 3)ln(x + 3) + (n-l)(x + l)ln(x + l) 

2(n + 8j 

+ 6xln2 - 9(x + l)ln3]. (1.83) 



J = i + e_lZ^ + e2_vii^(^2^22n + 52) + 0(e3). (1.81) 
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For H2{x, 1) we refer to the reviews of [ZSII^- Thus the e-expansion provides estimates 
of almost all quantities of interest as an asymptotic expansion in powers of e around the 
mean field values. Unfortunately it is impossible to estimate within this approximation 
scheme the corresponding error bars since the extrapolation to larger values of e is 
uncontrolled. Detailed analyses turn out that the critical exponents are more accurately 
estimated than the scaling functions and therefore the amplitudes. For example, the e^- 
approximation [Eq. (|1.81|) ] for the susceptibility exponent of the two-dimensional Ising 
model (n = 1, e = 2) yields 7 ~ 1.642. This value differs by 6% from the exact value 
7 = 7/4 [HHl- The amplitude ratio of the susceptibility [Eq. (|1.46j) ] can be expanded 
as [Hil 

x(-i,o) m ' 

proposing the estimate 81.14 . . . for e = 2. This result differs significantly (115%) from 
the exact value 37.69 . . . [HSllHni- The different accuracy reflects a conceptual difference 
between the universality of critical exponents and the universality of scaling functions. 
As pointed out clearly in |68j . the universality of exponents arises from the linearized 
RG flow in the vicinity of the fixed point, whereas the scaling functions are obtained 
from the entire, i.e., non-linear, RG flow. More precisely, the relevant trajectories 
from the fixed point of interest to other fixed points determine the universal scaling 
functions. This illuminates also why exponents between different universality classes 
may differ slightly while the scaling functions and therefore the amplitude combinations 
differ significantly. Thus deciding on a system's universality class by considering the 
scaling functions and amplitude combinations instead of critical exponents appears to 
be more sensitive. Hence, nothing demonstrates universality more convincing than the 
universal data-collapse of various systems (as shown in FigureOl and Figure0J) . 
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The minimization of the free energy is the governing principle of equihbrium statistical 
physics. In non-equilibrium situations no such framework exists. But nevertheless, the 
concepts and techniques which were developed for equilibrium critical phenomena can 
also be applied to non-equilibrium phase transitions. In particular, the order parameter 
concept, scaling, universality as well as renormalization group analyses are as powerful 
as in equilibrium, i.e., although they were developed in equilibrium theories their appli- 
cability is far beyond. Therefore, numerous non-equilibrium critical phenomena were 
successfully investigated by such methods, including the geometrical problem of perco- 
lation [HSl EZl EH] ) interface motion and depinning phenomena [531 IHOl 1^ > irreversible 
growth processes like diffusion-limited aggregation [HI] , as well as turbulence [HSl EU • 
In the following we focus our attention on a certain class of non-equilibrium phase 
transitions, the so-called absorbing phase transitions (APT) [HS]- These transitions 
were investigated intensively in the last two decades (see [71 [S] for recent reviews), 
since they are related to a great variety of non-equilibrium critical phenomena such as 
forest fires [HHl IHZI , epidemic spreadings in biology [HHl , catalytic chemical reactions in 
chemistry , spatio-temporal intermittency at the onset of turbulence |lfln| 11011 1102) , 
scattering of elementary particles at high energies and low-momentum transfer |103j . 
interface growth jl04j , self-organized criticality jl05l 11061 11071 1108j , damage spread- 
ing [3 1109[ lllOj , as well as wetting transitions [3 llllj . But like in equilibrium critical 
phenomena, most of the work focuses onto the determination of the critical exponents, 
neglecting the determination of the universal scaling functions. 



2.1 Definitions and the contact process at first glance 

Absorbing phase transitions occur in dynamical systems that are characterized by at 
least one absorbing state. Any configuration in which a system becomes trapped forever 
is an absorbing state. The essential physics of absorbing phase transitions is the com- 
petition between the proliferation and the annihilation of a quantity of interest A, for 
example particles, energy units, molecules in catalytic reactions, viruses, etc. Often the 
proliferation and the annihilation processes are described in terms of reaction-diffusion 
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schemes, e.g., 

prohferation : A — >2A, 2A — >3A, A — >3A,... (2.1) 

annihilation: A — >0, 2A — >0, 3A — >2A,.... (2.2) 

Obviously any configuration with zero density {pa = 0) is an absorbing state if spon- 
taneous particle generation processes 

— ^ A, — >2A, — >3A,... (2.3) 

do not take place. When the annihilation processes prevail the proliferation processes, 
the system will eventually reach the absorbing state after a transient regime. On the 
other hand, the system will be characterized by a non-zero steady-state density Pa 
in the thermodynamic limit if the proliferation processes outweigh the annihilation 
processes. In the latter case the system is said to be in the active phase whereas 
the absorbing phase contains all absorbing configurations. Let us assume that the 
competition between the proliferation and annihilation is described by a single rate r. 
An absorbing phase transition, i.e., a transition from the active phase to the absorbing 
phase takes place at the critical rate if the steady state density Pa vanishes below 

Pa = for r < r, , 

(2.4) 

Pa > for r > . 

In other words, the density Pa is the order parameter and the rate r is the control 
parameter of the absorbing phase transition. Analogous to equilibrium phase transi- 
tions, an order parameter exponent /? can be defined if the order parameter vanishes 
continuously 

Pa^ (r- ry . (2.5) 

It is worth mentioning that absorbing phase transitions have no equilibrium counter- 
parts since they are far from equilibrium per definition. In equilibrium, the associated 
transition rates between two states satisfy detailed balance |112| . In case of absorbing 
phase transitions the rate out of an absorbing state is zero. Therefore, an absorbing 
state can not obey detailed balance with any other active state |95j . 

For the sake of concreteness, let us consider the so-called contact process (CP) on a 
D-dimensional lattice. The contact process was introduced by Harris jll3j in order to 
model the spreading of epidemics (see for a review jll4j ). A lattice site may be empty 
(n = 0) or occupied (n = 1) by a particle representing a healthy or an infected site. 
Infected sites will recover with probability 1 — p corresponding to a process of particle 
annihilation. With probability p an infected site will create an offspring at a randomly 



2.1 Definitions and the contact process at first glance 



27 



chosen vacant nearest neighbor site reflecting particle propagation. Thus the dynamics 
of the contact process is described by the reaction scheme 

A — > 0, A — > 2A, (2.6) 

where the quantity A represents a particle. Particles are considered as active in the 
sense that they proliferate with rate p and annihilate with rate 1 — p. Empty lattice 
sites remain inactive and the empty lattice {p{A) = 0) is the unique absorbing state. 
In the following we denote the density of active sites, i.e., the order parameter of the 
absorbing phase transition with p^. 

In order to provide a deeper insight into absorbing phase transitions, we present 
a simple but intriguing mean field treatment. The probability that a given lattice is 
occupied is p^. This particle will be annihilated with probability 1 — p. Thus in an 
elementary update step the number of particles n is decreased with the probability 
p{An = — 1) = /9a(l — p)- On the other hand, a new particle is created with rate p at a 
randomly chosen vacant neighbor. The associate probability is p{An = 1) = p^p{\—p^). 
The number of particles remains unchanged with probability p(A?i = 0) = {1— p^)+ pl^p. 
Since we have neglected spatial fluctuations and correlations of p^, the above reaction 
scheme describes the simplest mean field theory of the contact process. This reaction 
scheme leads to the differential equation 

dtP. = PP.iX-p^) -Pa(l-p) = {2p-l)p, - ppI, (2.7) 
with the steady state solutions {dtP^ = 0) 

2p - 1 

Pa = V p, = — for p>0. (2.8) 

P 

The first equation corresponds to the absorbing state and is unstable for p > 1/2. The 
second equation yields unphysical results (p^ < 0) for p < 1/2 but describes the order 
parameter for p > 1/2. The behavior of the order parameter is sketched in FigureEl 
In the vicinity of the critical value p^ = 1/2 we find 

p. = 2 5p + 0{5p^) (2.9) 

with the reduced control parameter 5p = {p—Pc)/Pc- Thus the order parameter exponent 
is 13 = 1. 

Consider now the dynamical behavior of the order parameter. Solving Eq. (|2.7|) we 
obtain 

pM = , — (2.10) 

p + [bp/p^^t^o - p) exp [-dp t) 
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Figure 6: Sketch of the steady state behavior (left) of the order parameter of the contact process within 
the mean field approximation [Eq. (EHJ]- The dashed lines mark unphysical or unstable solutions. The 
right figure displays the dynamical behavior of the order parameter. The parameters of the presented 
mean field curves are Pa,t=o — 1, Sp — and Sp — ±0.02, respectively. 



Asymptotically {t — > oo) the order parameter behaves as 

6p 



p.it) 



Sp<0 



Sp>0 



dp [ p 



5p 
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6p 
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p 



p 



5p 

Pa,t = 



(2.11) 
(2.12) 



Apart from criticality, the steady state solutions {p^ = 5p/p and = 0) are approached 
exponentially, independent of the initial value /5a,t=o- The associated correlation time 
diverges as At criticality {6p = 0) the order parameter exhibits an algebraic 

decay 

P.it) = — — - (2.13) 



Pa,t = 



pt 



The dynamical behavior of the order parameter is sketched in FigureEl 

Similar to equilibrium phase transitions, it is often possible for absorbing phase 
transitions to apply an external field that is conjugated to the order parameter. Being 
a conjugated field it has to destroy the absorbing phase, it has to be independent of 
the control parameter, and the corresponding linear response function has to diverge 
at the critical point 

5 Pa 



X 



dh 



oo . 



(2.14) 



For the contact process the conjugated field causes a spontaneous creation of particles 
(0 — > A). Clearly spontaneous particle generation destroys the absorbing state and 
therefore the absorbing phase transition at all. In the above mean field scheme the 
conjugated field creates a particle at an empty lattice site with probability (1 — p^)h. 
This particle will survive the next update with probability p or will be annihilated with 
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probability 1—p, respectively. The modified differential equation [Eq. (|2.7j) ] is given by 

dtp^ = 6pp^ - ppI + p{l- p^)h. (2.15) 

The steady state order parameter is a function of both the control parameter and the 
conjugated field 

The solution with the + sign describes the order parameter as a function of the control 
parameter and of the conjugated field whereas the — sign solution yields unphysical 
results Pa < 0. Both solutions are sketched in FigureEl for various field values. As 
requested, the absorbing phase is no longer a steady state solution for h > 0. At 
criticality {5p = 0) the order parameter behaves as 

p^(Sp = 0,h) = Vh + 0{h). (2.17) 

Now we examine the order parameter behavior close to the critical point. Therefore 
we perform in Eq. (|2.15j) the limits p^ ^ 0, (5p — > 0, and h ^ with the constraints 
that p^l\fh and p^/5p are finite. The remaining leading order yields 

The Eq. (|2.9|) and Eq. (|2.17|) are recovered from this result by setting /i = and 6p = 0, 
respectively. It is straight forward to derive the dynamical behavior of the order param- 
eter p^{6p,h,t) close to the critical point (p^ ^ !)• The differential equation Eq. ()2.15|) 
yields that the steady state solution is approached from above as 

P.iSp,h,t) ~ g + ^{^y + h (l±2coe-*/€li + 0(e-2*/«ll)) , (2.19) 

where the constant Cq contains the initial conditions. The + sign corresponds to initial 
conditions with Pa,t=o > Pa,t-.oo whereas the — sign is valid for Pa,t=o < Pa.t-oo- Again, 
the steady state value is exponentially approached independent of the initial condition. 
The corresponding temporal correlation length is usually denoted as for absorbing 
phase transitions, and we find jll5j 

^,i6p,h) = , \ (2.20) 

Thus the correlation time diverges at the critical point as 

Cll(5p,/i = 0) = \6p\-^ and ^||(<5p = 0,/i) = ^ h'^/^ . (2.21) 

z p 
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Figure 7: The mean field order parameter (left) and the mean field order parameter susceptibility 
(right) of the contact process. The dashed lines mark unphysical or unstable solutions. Thick lines 
correspond to the zero-field behavior whereas the thin lines sketch the smooth non-zero field behavior. 
The susceptibility displays a finite peak for non-zero field. In the limit /i — > this peak diverges 
signalling the critical point. 



The zero-field singularity is usually associated with the temporal correlation length 
exponent z^n = 1. 

Furthermore, the derivative of with respect to the conjugated field yields the 
order parameter susceptibility 



dh 



6p\' 
2pJ 



1 -1/2 



+ h 



(2.22) 



As required, the order parameter susceptibility becomes infinite at the transition point 



oc 



oc 



1-1 



oo , 



OO . 



(2.23) 
(2.24) 



The mean field behavior of the susceptibility is displayed in Figure[7| 

It is possible to incorporate spatial variations of the order parameter in the mean 
field theory. Therefore Eq. (|2.15|) has to be modified and we assume that the order 
parameter obeys close to criticality (p^ <^ 1) 

dtp.ix,t) = Spp,ix,t) - pp,{x,tf + ph + rvVfei), (2.25) 
where Ps.{x^t) represents the local particle density and where we add the diffusive cou- 



pling with r > 0. Usually, continuous equations like Eq. (|2.25|) are obtained from 
microscopic models by a coarse graining procedure or they are phenomenologically 
motivated by a Landau expansion [301. According to Landau, a given system is de- 
scribed by an appropriate functional containing all analytic terms that are consis- 
tent with the symmetries of the microscopic problem. Consider small spatial devi- 
ations from the homogeneous steady state value p^{5p,h) [Eo. (|2.18|) ] by introducing 
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SPc^ix, t) = Pa^ix, t) — p^{5p, h). Tlie order parameter variations 5p^{x, t) obey the differ- 
ential equation 



dtdp^{x,t) = [5p - 2pp^i6p,h)] 5p^{x,t) + TV^5p^ix,t) + 0(5 pi 



(2.26) 



where nonlinear terms are neglected. Working in Fourier space and using Eq. ()2.18|) we 
find that small deviations with wavevector k decay as 



dtSp,{k, t) 
Here, the characteristic length 



5pSht) 



2pJ 



+ h 



-1/4 



(2.27) 



(2.28) 



is introduced that describes spatial correlations of the order parameter variations |llf)| 
I117j . Approaching the transition point, the spatial correlation length becomes infinite 



T\6p\ 



-1/2 



and 



C^{6p = 0,h) 




(2.29) 



Thus the mean field value of the exponent of the correlation length at zero field is 
= 1/2. Note that the temporal and spatial correlation length are related via jll5j 

c2 



(2.30) 



Solving Eq. 1)2. 27() and taking the inverse transform we obtain the order parameter 
variations 



5p^{x, t) 



Pa.O 



1 



D/2 



-t/iw p-(2i-;^o)V4rt 



(2.31) 



.47rrt 

where we have assumed a seed-like initial density variation 6p^{x, t = 0) = /5a,o^(iE — Xq)- 
Similar to the homogeneous case the lifetime of the localized density fluctuation is 
determined by the temporal correlation time . Furthermore the activity diffuses over 
regions of distance V^Tt from the origin x^. 

Fluctuations within the steady state can be investigated by adding an appropriate 
noise term r]{x, t) in Eq. (|2.25|) representing rapidly-varying degrees of freedom (see 
e.g. |115l Ill6] ). The noise has zero mean and the correlator is assumed to be given by 



( r/(x, i) r/(x', t') ) = K 5(x - x') 5(t - ' 



(2.32) 



In that case, the mean field steady state fluctuations are given by jll6j 
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For zero-field the fluctuations reduce to 

f if 5p<0 
{5pi) oc (2.34) 
[ K/2p if (^p > , 

i.e., the fluctuations do not diverge at the critical point but exhibit a finite jump. 

In summary, the mean field theory presented above describes the absorbing phase 
transition of the contact process. The steady state and dynamical order parameter 
behavior have been derived and certain critical exponents are determined. In contrast 
to the one characteristic length scale of ordinary critical equilibrium systems, absorbing 
phase transitions are characterized by two length scales, and ^y, with different critical 
exponents and i^y . This scaling behavior of directed percolation can be interpreted 
in terms of so-called strongly anisotropic scaling jll8j . where a different steady state 
scaling behavior occurs along different directions (indicated by different critical expo- 
nents, V\\/vi_ 7^ 1). Paradigms of strongly anisotropic scaling are Lifshitz points jll9j . 
e.g. in Ising models with competing interactions |12ni 1121] . 

Unfortunately, the notations of some critical exponents for absorbing phase tran- 
sitions differ from those of equilibrium systems. In particular the field dependence of 
the order parameter at the critical isotherm is written as [Jj 

p,((5p = 0,/i) 6c /i^/'", (2.35) 

in contrast to the usual equilibrium notation Moch^^^. Thus the field exponent a 
corresponds in equilibrium to P 6, i.e., the exponent a of absorbing phase transitions is 
identical to the gap exponent A = P S oi equilibrium systems. In order to avoid friction 
with the literature of absorbing phase transitions, where 6 is reserved to describe the 
dynamical scaling behavior, we will use the notation of Eq. (|2.35j) in the following. 
Eventually we summarize that the mean field theory is characterized by the exponents 
/? = 1, (T = 2, z^li = 1, = 1/2, as well as 7 = 1, where the latter one describes the 
divergence of the zero-field susceptibility (xoc l^pl""^). 

2.2 Numerical simulation methods 

Despite the simplicity of the contact process, a rigorous solution is still lacking. Like in 
equilibrium the mean field theory provides correct results only above the upper critical 
dimension D^. It fails below since different sites are strongly correlated whereas they 
are considered as independent within the mean field approach. Beside of approximation 
schemes like series expansions [1221 II2S1 II211 ESI EEl 1123 II2H1 1123 11313 HHH II321 
I133j and e-expansion within renormalization group approaches, numerical simulations 
were intensively used in order to determine the critical behavior of the contact process 
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and related models. There are two primary ways to perform numerical simulations 
of lattice models exhibiting absorbing phase transitions. The first method is to study 
the steady state behavior of certain quantities of interest, e.g., the order parameter 
and its fluctuations as a function of various parameters like the control parameter and 
the external field (e.g. |134| 11351 1136[ [T37| ). Especially the equation of state as well 
as the susceptibility can be obtained in that way [ISHl UMl EDI ED HH ESI- But 
approaching the critical point, the simulations are affected by diverging fluctuations 
and correlations causing finite-size and critical slowing down effects (see e.g. jll2j ). 
Techniques such as finite-size scaling analysis have to be applied in order to handle 
these effects. 

In contrast to steady state measurements, the second method is based on the dy- 
namical scaling behavior. Often the evolution of a prepared state close to the absorbing 
phase is investigated. For example the spreading of activity started from a single seed 
is studied in the vicinity of the critical point |144j . This method is known to be very 
efficient and the critical exponents can be determined with high accuracy (see for in- 
stance ElESlEniEZlEHlEnilinnilinil)- a drawback of this technique is that it is 
restricted to the vicinity of the absorbing state. For example the conjugated field drives 
the systems away from the absorbing phase and can therefore not be incorporated. In 
the following both methods are discussed and the associated exponents are defined. For 
the sake of concreteness corresponding data of the five-dimensional contact process are 
presented as an exemplification and are compared to the mean field results. 

2.2.1 Steady state scaling behavior 

In order to perform simulations of a system exhibiting an absorbing phase transition 
we have to specify the lattice type and the boundary conditions. In case of the contact 




Figure 8: The density of active sites as a function of time for the five-dimensional contact process 
(L = 16, 5p — 0.02, and pa,t=o ~ !)• After a transient regime, which depends on the initial configuration, 
the density of active sites fiuctuates around the steady state value (dashed line). 
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process it is customary to consider D-dimensional simple cubic lattices of linear size L 
and periodic boundary conditions. Steady state simulations usually start far away 
from the absorbing state, e.g. with a fully occupied lattice. The system is updated 
according to the microscopic rules presented on page EZI using a randomly sequential 
update scheme. Therefore all occupied sites are listed and one active site after the other 
is updated, selected at random. After a sufficient number of update steps the system 
reaches a steady state where the number of active sites fluctuates around the average 
value (see Figure|HI). It is customary to interpret one complete lattice update with 
active sites as one time step [t ^ t + 1). Thus an elementary update of one lattice site 
corresponds to a time increment t ^ t + Monitoring the density of active sites 

Pa.{t) in the steady state one obtains an estimate of the order parameter = {p^{t)) as 
well as of its fluctuations 

Ap. = ((p.(t)2) - , (2.36) 

where (...) corresponds to the temporal average 

{pM') = \;jlp.{tf- (2.37) 

t=i 

Of course reliable results are only obtained if the number of update steps significantly 
exceeds the correlation time T » ^y. This becomes notoriously difficult close to the 
critical point since — > oo. Furthermore the accuracy can be improved if an additional 
averaging over different initial configurations is performed. 

The average density of active sites as well as its fluctuations are plotted in FigurelHl 
as a function of the control parameter p for various system sizes L. As can be seen the 
order parameter tends to zero in the vicinity of p ~ 0.532. Assuming that the scaling 
behavior of p^ obeys asymptotically the power law 

Pa(p) 6c (p-p,)^ (2.38) 

the critical value is varied until a straight line in a log-log plot is obtained (see Fig- 
urelH)). Convincing results are observed for p^ = 0.53237±0. 00006 and the corresponding 
curve is shown in Figurel^l For p^ = 0.53231 and p^ = 0.53243 significant curvatures in 
the log-log plot occur. In this way it is possible to estimate the critical value p^ as well 
as its error-bars. Then a regression analysis yields the value of the order parameter 
exponent (3 = 0.991±0.015 in agreement with the mean field value (3 = 1. 
Usually the order parameter fluctuations diverge at the transition point 

Ap, 6c (p-p,)^^'. (2.39) 

But as can be seen the fluctuations are characterized by a finite jump at p^ correspond- 
ing to the value 7' = 0. This agrees well with the mean field result Eq. 1)2.34(1 . We will 
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Figure 9: The order parameter (upper left) and its fluctuations (upper right) of the five- dimensional 
contact process for various system sizes L. The data are averaged over at least T = lO'^ lattice updates. 
The lower figures show the order parameter as a function of p — Pc • The dashed line corresponds to 
the mean field behavior with (3 — 1. A straight line is obtained for Pc = 0.53237, whereas Pc ~ 0.53243 
and Pc — 0.53231 lead to significant curvatures (see lower right figure). For the sake of simplicity lines 
are plotted instead of symbols. 

see in the next chapter that low-dimensional systems are characterized by a non-zero 
exponent. 

The effects of finite system sizes deserve comment. Of course a finite system reaches 
the absorbing state with a finite probability even for p > p^. But this probability tends 
to zero for L — > oo. It turns out that sufficiently large systems will never reach the 
absorbing phase within reasonable simulation times for p > p^. The situation changes 
close to the transition point since the finite system size L prevents the correlation 
length from becoming infinite. As a result so-called finite-size effects occur, i.e., 
the corresponding singularities become rounded and shifted (see for example jll2p . A 
typical feature of finite-size effects in equilibrium systems is that a given system may 
pass within the simulations from one phase to the other. This behavior is caused by 
critical fluctuations that increases if one approaches the transition point. In case of 
absorbing phase transitions the scenario is different. Approaching the transition point 
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Figure 10: The order parameter (upper left) and its fluctuations (upper right) of the five- dimensional 
contact process for various field values (^ = 3 10~^, 8 10~^, 10~^). The order parameter and the fluctu- 
ations exhibit an analytic behavior for non-zero fleld. The order parameter susceptibility (lower left) 
displays the characteristic peak at the critical point (dashed line). The lower right flgure shows the 
order parameter as a function of the conjugated field at p = Pc for various system sizes L = 4, 8, 16 
(from right to left). The dashed line corresponds to the mean field behavior with 13/ a = 1/2. 



the correlation length ^ increases and as soon as ^ is of the order of L the system 
may pass to an absorbing state and is trapped forever. For example finite-size effects 
are expected to occur for dp < TL~^ according to the mean field result Eq. (|2.29() . In 
that case the absorbing state is reached by finite fluctuations even for p > p^. A simple 
way to handle that problem is to increase the system size before finite-size effects 
occur (see Figurel^I)- Considering overlapping data regions for different system sizes 
ensures that the obtained results are not affected by finite-size effects. Additionally it 
is suggested (see, for instance [135 ) to consider metastable or quasisteady-state values 
of the order parameter. But as pointed out in |152j this method is inefficient and can 
be seriously questioned since the metastable values are not well defined. 

A better way to analyze finite-size effects for absorbing phase transitions is to incor- 
porate the conjugated field. Due to the conjugated fleld the system cannot be trapped 
forever in the absorbing phase. Therefore steady state quantities are available for all 
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values of the control parameter. Analogous to equilibrium phase transitions, the con- 
jugated field results in a rounding of the zero-field curves, i.e., the order parameter and 
its fluctuations behave smoothly as a function of the control parameter p for finite field 
values (see FigureEU- For /i — > we recover the non-analytical behavior of and 
A/9;,, respectively. Furthermore, the order parameter susceptibility x can be obtained 
by performing the numerical derivative of the order parameter Ps.{5p, h) with respect to 
the conjugated field [Eq. ()2.14|) ]. Similar to the mean field solution the susceptibility 
displays the characteristic peak at the critical point p^. In the limit /i ^ this peak 
diverges, signalling the critical point. 

The order parameter and fluctuation data presented in Figure lTUl are obtained from 
simulations where the correlation length is small compared to the systems size. Thus 
these data do not suffer from finite-size effects. But approaching the transition point 
finite-size effects emerge. Figure^] shows the order parameter at p^ as a function of the 
conjugated field. The behavior of the infinite system [L S> ^) agrees with the mean 
field result p^cx/i^/^. Approaching the transition point {h — > 0), finite-size effects, i.e., 
deviations from the behavior of the infinite system occur. The larger the system size L 
the larger is the scaling regime where the power-law behavior is observed. In this way it 
is possible to study finite-size effects of steady state measurements in absorbing phase 
transitions. A detailed analysis of the order parameter, the fluctuations, as well as 
of a fourth order cumulant is presented in the next chapter. Similar to equilibrium 
phase transitions, we formulate scaling forms which incorporate the system size as an 
additional scaling field. 

2.2.2 Dynamical scaling behavior 

In the following we investigate the dynamical scaling behavior at zero field close to 
the critical point. First, we consider how the order parameter decays starting from 
a fully occupied lattice, i.e., starting from a so-called homogenous particle source. 
Studying the temporal evolution of the system one has to average over different runs 
Pa(i) = [Pa(t,n,„d)] with 



Here, denotes the seed number of the random number generator. Simulation results 
for various values of the control parameter are shown in Figure lTTl 

Similar to the mean field behavior [Eqs. (|2. 11112. the order parameter tends 
exponentially to the steady state value above p^ (positive curvature in the log-log plot). 
Below the transition point the absorbing state is approached exponentially (negative 
curvature). At the critical point the order parameter decays asymptotically as 





Pa(i) oc t 



(2.41) 
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As can be seen in Figure lTTl the decay exponent agrees with the mean field value a = 1. 
Again, deviations from the power-law behavior allow to estimate the critical value p^. 
The obtained result and its uncertainty is in agreement with the steady state measure- 
ment = 0.53237 ± 0.00006. 

More accurate estimates for the critical value p^ and for the critical exponents 
can be obtained by considering the activity spreading generated from a single active 
seed |144j . a so-called localized particle source. Three typical snapshots of activity 
spreading of the one-dimensional contact process are shown in Fieure lT^ Below the 
critical point the activity ceases after a certain transient regime whereas above p^ the 
activity spreading continues infinitely with a finite probability. At the transition point 
a cluster is generated by the temporal evolution that exhibits fractal properties |145j . 
Interpreting active sites, e.g., as viruses this figure illustrates the spreading of epidemics 
through a population. 

In order to analyze the scaling behavior of activity spreading it is useful to measure 
the survival probability i-'a(t), the number of active sites N^{t) as well as the average 
mean square distance of the spreading R^{t)'^ jl44j . Here, N^{t) is averaged over all 
clusters whereas R^{t)'^ is averaged only over surviving clusters. At criticality, the 
following power-law behavior are expected to occur asymptotically 

P^{t) 6c t-^ , N^{t) 6c , R^itf 6c t^A , (2.42) 
10" 
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Figure 11: The dynamical behavior of the order parameter of the five-dimensional contact process for 
L = 32. The parameters of the presented curves are 5p = 0, 5p — ±0.02 (thick lines), Sp = ±0.00006 
(thin lines) respectively. The data are averaged over at least N — 10^ different runs. Started from a 
fully occupied lattice (pa,t=o ~ 1) the order parameter saturates for p > p^ and decays exponentially 
for p < Pa- At the critical point the order parameter displays an algebraic decay. The dashed line 
corresponds to the mean field exponent a — 1 and is shown to guide the eye. 
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Figure 12: Snapshots of activity spreading in the one-dimensional contact process (time increases 
downward from t = to 10'^). The parameters of the presented clusters are 5p — and Sp — ±0.02. 

where we have used the notation of [7!. The so-called dynamical exponent z describes 
the relationship between spatial and temporal correlations and is therefore identified 
as 

z = ^. (2.43) 

According to Eq. H2.3U() the mean field value is given hy z = 2. The survival prob- 
ability P^{t) is shown in FigureEl Again, a positive curvature for t ^ oo indicates 
in a log-log plot the active phase whereas a negative curvature indicates the absorb- 
ing phase. Similar to the steady state [Eq. H2.38|l ] and dynamical [Eq. (|2.4H) ] order 
parameter behavior, the survival probability obeys power-law behavior Eq. 1)2. 42() only 
asymptotically, i.e., confluent singularities have to be taken into consideration. Often, 
local slopes are considered 

In c 

to handle this problem. The parameter c defines the distance of the local slopes and 
typical values used in previous studies are for example c = 4 jl53j , c = 5 jl491 1154j , 
c = 8 jl45j as well as c = 10 jl55j . Of course, for c — > 1 the local slopes equals the 
so-called effective exponent |156) 

dlnP^jt) 

5.«(t) = -^^. (2.45) 

The analysis of the local slope is illustrated in Figure lTTn Above criticality the local slope 
tends to zero whereas (5iocai(i) oo for 6p < 0. Both limiting cases are separated by 
the local slope at criticality and the exponent 6 can be estimated from an extrapolation 
to f ^ oo. 

The exponents 9 and z can be determined in an analogous way and the mean field 
values 6 = 1, 9 = 0, and z = 2 are obtained in case of the five dimensional contact 
process. The dynamical scaling behavior at equals a critical branching process (see 
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Figure 13: The survival probability of the five-dimensional contact process close to criticality (left 
figure). The parameters of the presented data are &p = G (straight asymptotic line) and &p = ±0.02 
(curved lines). The dashed line corresponds to the asymptotic mean field behavior (5 = 1). The right 
figure shows the local slope 5iocai [Eq. 12.4411 1 . The parameters of the presented data are &p = (thick 
line) and &p = ±0.00002 (upper and lower thin line). The data are averaged over A*' = 10^ independent 
runs and c = 8 is used to calculate the local slope. As can be seen the mean field value (5 = 1 is 
approached for 5p = Q and t —* oo, whereas the local slopes tend to or oo for 5p ^ 0. 



appendix) and the value ^ = corresponds to N^{t) ~ const. The accuracy of the 
exponents can be checked using the hyperscahng relation [14^; 

e + 26 = (2.46) 

which is valid for D < D^. A derivation of this scaling law is presented in the next 
chapter. Notice that the mean field values 6 = 1,6 = 0, and z = 2 fulfill the scaling law 
Eq. (|2.46|) for = 4, indicating that four is the upper critical dimension of the contact 
process. 

Other exponents can be derived from 5, 9, and z by additional scaling laws. For 
example, the average number of active sites per surviving cluster scales as A'^^,sur(i) = 
N^{t) / P^{t) 6c As usually, the fractal dimension Df is defined via N^^s^,.{t) 6c R{t)^^ 
yielding |145j 

Dt = z{e + 6). (2.47) 

In case of the one-dimensional contact process, the critical clusters are characterized 
by the fractal dimension Df ~ 0.75, whereas the mean field values yield Df = 2. 

In summary, activity spreading measurements provide very accurate estimates for 
the critical value as well as for the exponents 5, 9, and z of systems exhibiting absorb- 
ing phase transitions (see e.g. |ll[I13[IiaiI171IIlHl[IiniIIMlIinil)- The efficiency of 
this technique can be increased significantly by performing off-lattice simulations using 
the fact that clusters of active sites are fractal at the critical point, i.e., huge parts of 
the lattice remains empty. Storing just the coordinates of active sites in dynamically 
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generated lists makes the algorithm much more efficient compared to steady state ex- 
ponents. Another advantage of such ofF-lattice simulations is that no finite-size effects 
occur. But a drawback of activity spreading measurements is that the determination 
of the exponents 5, 6, and z via Eq. (|2.42|) is not sufficient to describe the scaling be- 
havior of absorbing phase transition. Due to the hyperscaling relation [Eq. (|2.46j) ] only 
two independent exponents can be determined. But we will see in the next chapter 
that the scaling behavior of the contact process is characterized by three independent 
exponents. Thus additional measurements have to be performed in order to complete 
the set of critical exponents. 

Additionally to the consideration of lattice models, direct numerical integrations 
of the corresponding Langevin equations are used to investigate the behavior of the 
contact process or other related models jl57| I158( I159| . But this technique runs into 
difficulties, for example negative densities occur if the absorbing phase is approached. 
To bypass this problem, a discretization of the density variable was proposed in jl57j . 
This approach has been applied to various absorbing phase transitions. The obtained 
estimates of the exponents agree with those of simulations of lattice models but are of 
less accuracy. 
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The problem of directed percolation was introduced in the mathematical literature in 
1957 by Broadbent and Hammersley to mimic various processes like epidemic spreading, 
wetting in porous media as well as wandering in mazes jl6Uj . In a lattice formulation, 
directed percolation is an anisotropic modification of isotropic percolation. Consider 
for instance the problem of bond directed percolation as shown in Figure lTH Here, 
two lattice sites are connected by a bond with probability p [Sni EZl EHj • If the bond 
probability is sufficiently large, a cluster of connected sites will propagate through the 
system. The probability that a given site belongs to a percolating cluster is the order 
parameter of the percolation transition and obeys 



for p > Pc, whereas it is zero below the critical value p^. Due to a duality symmetry 
the critical value of bond percolation on a square lattice is p^ = 1/2. 

In case of directed percolation (DP) a preferred spatial direction is introduced. For 
example, only branching of the cluster in the downward direction in FigureEl is al- 
lowed. Similar to isotropic percolation one considers the probability Ppcrc(p) that a 
path of connected sites percolates through the entire system. Again, depending on 
the probability p, a continuous phase transition takes place at a certain threshold p^ 
and the percolation probability Ppaidp) vanishes algebraically [Eg. p. II) ]. Thus both 
percolation transitions are characterized by the order parameter Pp^icip), but as ex- 
pected the critical value for directed percolation (p^ = 0.644700185(5) |133j ) is larger 
than the isotropic value p^ = 1/2. Another important difference is that isotropic per- 
colation can be mapped exactly to the equilibrium q-state Potts model jl611 1162| in 
the limit g — > 1 jl63| 116^ and the critical exponents are known exactly, /?' = 5/36 and 
1/ = 4/3 |165|I166'] . On the other hand, directed percolation is a non-equilibrium system 
and is not exactly solved so far. Furthermore, the critical exponents do not seem to be 
given by rational numbers in contrast to isotropic percolation. And finally the upper 
critical dimension of isotropic percolation is = 6 (HE] whereas = 4 jl()3| I167j in 
case of directed percolation. In other words, despite their similarities isotropic and 
directed percolation belong to different universality classes. 
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Figure 14: Isotropic bond percolation (left) and directed bond percolation (right) on a rotated square 
lattice. Only the cluster is shown that is connected to the marked lattice site (bold circle). The same 
bond configuration leads to different clusters for isotropic and directed percolation. In case of directed 
percolation the cluster propagates only along the preferred direction (arrows). Therefore, each site 
generates an individual cluster for a given bond configuration. 

Although, directed percolation is a geometrical problem according to the original 
definition, it appears in different contexts as a dynamic process. In those cases the 
preferred direction is interpreted as a representation of time instead of space. There- 
fore, directed percolation describes the propagation or spreading of a non-conserved 
quantity (particles, viruses etc.) through a one-dimensional cellular automata. Each 
row represents the state of the system at time step t and particles are associated to 
lattice sites that belong to a cluster at that time step. Obviously, an empty row is an 
absorbing state if no spontaneous particle creation takes place. In this way the geo- 
metrical percolation transition of two-dimensional directed percolation is mapped to an 
absorbing phase transition of a one-dimensional cellular automaton. Furthermore, the 
percolation order parameter [Eq. (|33)] corresponds to the ultimate survival probability 

P^Up) = lim P.{t,p) (3.2) 

t^oo 

that a cluster generated by an active single seed propagates through the entire sys- 
tem. Thus the critical behavior of directed percolation can be described by both the 
percolation probability and the steady state density order parameter 

Pp,,,(p) 6c (p-Pc)^', /Oa(p) 6c {p-PcY ■ (3.3) 

Today, directed percolation is considered as the most important universality class 
of absorbing phase transitions. Directed percolation transitions occur, for example, 
in systems describing the spreading of epidemics like the contact process, catalytic 
reactions, the percolation in porous media, etc. (see for recent reviews [7| |H]). The 
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widespread occurrence of directed percolation is reflected by the universality hypoth- 
esis of Janssen and Grassberger [231 124j : Short-range interacting systems, exhibiting 
a continuous phase transition into a unique absorbing state generically belong to the 
universality class of directed percolation, provided they have no additional symme- 
tries. Different universality classes are expected to occur in the presence of additional 
symmetries, e.g. particle-hole symmetry, or relevant disorder effects. 

Similar to equilibrium critical phenomena, the universality of directed percolation is 
understood by renormalization group treatments of an associated continuous field the- 
ory. Amazingly the field theoretic description of directed percolation is well established 
in high energy physics, termed Reggeon field theory (see e.g. jlU31ll68] ^. Reggeon field 
theory describes hadron scattering at high energies and low-momentum transfer. The 
time direction of directed percolation is related to the longitudinal momentum, the 
spatial direction equals the impact parameter and the total cross section corresponds 
to the square of the percolation probability, i.e., to the order parameter jlU3j . 

The process of directed percolation might be represented by the Langevin equa- 
tion |2n] 

dtp.ix,t) = rp,{x,t) - uplix,t) + TV^p,{x,t) + r]{x,t) (3.4) 

which describes the order parameter p^{x, t) on a mesoscopic scale and where t] denotes 
the multiplicative noise which accounts for fiuctuations of the particle density p^{x,t). 
According to the central limit theorem, r](x, t) is a Gaussian random variable with zero 
mean and whose correlator is given by 

{r]{x,t)r]{x',t')) = K p,{x,t) 6{x- x') 5{t- t') . (3.5) 

Notice, that the multiplicative noise ensures that the systems is trapped in the ab- 
sorbing state Pi^ix, t) = 0. Furthermore, higher order terms as Pa.{x, t)^, Pa.{x, t)^, ... (or 
V^Pa(^) i)) V^/3a(x, i), . . .) are irrelevant under renormalization group transformations 
as long as u > 0. Negative values of u give rise to a first order phase transition whereas 
u = is associated to a tricritical point |169[ llTDj similar to the Landau functional of 
ferromagnets [Eq. H1.34() ] with terms of order tvP |3Uj . 

Obviously, the above Langevin equation corresponds to the mean field equation of 
the contact process [Eq. (|2.25)) ] if the noise can be neglected. More precisely, the mean 
field equation describes the macroscopic behavior correctly if the noise is irrelevant 
on large scales, whereas it fails if the noise is relevant on large scales. The question 
whether the noise is relevant or irrelevant can be answered by the so-called dimensional 
analysis. Similar to the renormalization group approach in equilibrium (see II. 4|) one 
considers how the coupling constants r, u, T, k of the Langevin equation behave under 
the rescaling transformations (see e.g. [71 

x^bx, t^¥t, p. ^ b-f^/"^ p, , (3.6) 



46 



Directed percolation 



with 6 > 1. Replacing x, t, and in Eq. (|3.4j) and Eq. (|.S.5|) by the rescaled quantities 
yields the recursion relations for the couplings constants 

r ^ b'r, u ^ h^-f^/^^u, r ^ 6^-2 r, k ^ lP/''^+^-^ k. (3.7) 

Thus the mean field solution is consistent for D > A since the Langevin equation is 
at criticality (r = 0) invariant under the transformations Eq. I|3.6j) for the mean field 
values /? = 1, z^x = 1/2, z = 2. The situation is different below the upper critical 
dimension = A where the noise is relevant and the mean field solution does not 
describe the macroscopic behavior of the system. 

Below the upper critical dimension, renormalization group techniques have to be 
applied to determine the critical exponents and the scaling functions. In that case path 
integral formulations are more adequate than the Langevin equation (see e.g. jl71j ^. 
Stationary correlation functions as well as response functions can be determined by 
calculating path integrals with weight exp(— 5), where the dynamic functional S de- 
scribes the considered stochastic process. Up to higher irrelevant orders the dynamic 
functional associated to directed percolation is given by 11721 11731 1174j 



5[n,n] =/d-^dtn 



dtU — (r + V^)n — ^ —h — un^ n 



(3.8) 



where the particle density equals n{x, t) and where n{x, t) denotes the response field 
conjugated to the Langevin noise field |175j . This functional corresponds to the La- 
grangian of Reggeon field theory |168j . Rescaling the fields 

TT-ix-it) = fJ'S{x,t), n{x,t) = fi^^ s{x,t) (3-9) 

the functional S is invariant under the duality transformation (so-called rapidity rever- 
sal in Reggeon field theory) 

s{x,t) < — > -s{x,-t) (3.10) 

for fj? = 2u/k. Note that is a redundant variable from the renormalization group 
point of view 1^51 1176j . As usual, the duality transformation defines a dual stochastic 
process that might differ from the original one |177j . Close to criticality, the average 
density of particles of the dual process is connected to the survival probability |176j 
via 

PM ~ /iVr'(i)- (3.11) 
Here, Pa.{t) denotes again the probability that a cluster generated by a single seed is 
still active after t time steps. On the other hand, Pa^^'H^) describes the particle decay 
of the dual process started from a fully occupied lattice. The self-duality of directed 
percolation implies Pa{t) = pt^''\t) and therefore |144| 1176] 

P.it) ~ fi^pM- (3.12) 
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Figure 15: The asymptotic equivalence of pa{t) and Pa(t) for the five dimensional contact process. 
The survival probability (here multiplied by a factor) measures the probability that a cluster 

generated by a single seed is still active after t time steps. On the other hand Pa.{t) describes the decay 
of the particle density started from a fully occupied lattice. 



Numerical simulations confirm this result, as can be seen in FigureEl for the five 
dimensional contact process. The asymptotic equivalence ensures that both quantities 
have the same exponents |144) . including 



(3.13) 



It is worth mentioning that Eq. (|3.12|) follows from the rapidity reversal symmetry of 
the dynamical functional 5[n, n], i.e., it is a specific property of the directed percolation 
universality class. Thus, compared to general absorbing phase transitions, the number 
of independent critical exponents for directed percolation is reduced. Furthermore, the 
self-duality is restricted to the field theoretical treatment (e.g. after a course-graining 
procedure) and does not necessarily represent a physical symmetry of microscopic mod- 
els US]- 

Using standard renormalization group techniques such as e-expansion, the critical 
exponents as well as the scaling functions were obtained. For example the exponents 
P, z, and uj_ are in two loop order pi ITTsl [TTol ITW)] 
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with € = 4 — D. Furthermore, the Widom- Griffiths scahng form of the equation of state 
is given close to the transition point (to be specific for x < and |x| <C 1) by jl81j 
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Below, these universal quantities will be compared to various lattice models belonging 
to the directed percolation universality class. 



3.1 Lattice models of directed percolation 

In the following we consider various lattice models that belong to the universality class 
of directed percolation. First, we revisit the contact process that is well known in 
the mathematical literature. Second we consider the Domany-Kinzel (DK) cellular 
automaton jl82j which is very useful in order to perform numerical investigations of 
directed bond and directed site percolation. Furthermore, the Domany-Kinzel automa- 
ton exhibits a non-trivial phase diagram and allows therefore to study how non-directed 
percolation behavior arises if an additional symmetry occurs. Third we consider the 
pair contact process (POP) |134j that is characterized in contrast to the other models 
by infinitely many absorbing states. In contrast to the first two models the universal 
behavior of the pair contact process is still a matter of controversial discussions in the 
literature. In particular the scaling behavior in higher dimension and the effects of addi- 
tional particle diffusion are investigated. Furthermore we brieffy discuss the threshold 
transfer process |136j as well as the Ziff-Gulari-Barshad model [22] • The latter one 
mimics the catalysis of the carbon monoxide oxidation. 

3.1.1 Contact process 

The contact process (CP) is a continuous-time Markov process that is usually defined on 
a D-dimensional simple cubic lattice. A lattice site may be empty (n = 0) or occupied 
(n = 1) by a particle and the dynamics is characterized by spontaneously occurring 
processes, taking place with certain transition rates. In numerical simulations the 
asynchronous update is realized by a random sequential update scheme (see section 1221 
for details): A particle on a randomly selected lattice site is annihilated with rate one, 
whereas particle creation takes places on an empty neighboring site with rate \n/2D, 
i.e., 

n, = 1 -j^ = 0, (3.19) 
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ni = 



\n/2D 



ni = l 



(3.20) 



where n denotes the number of occupied neighbors of n^. Notice that the rates are de- 
fined as transition probabilities per time unit, i.e., they are not normalized probabilities 
and may be larger than one. Thus rescaling the time will change the transition rates. 
In simulations a discrete time formulation of the contact process is performed. In that 
case a particle creation takes place at a randomly chosen lattice site with probability 
p = A/(l + A) whereas particle annihilation occurs with probability 1 — p = 1/(1 + A). 
In dynamical simulations, the time increment is associated to each attempted 

elementary update step, where A''^ denotes the number of active sites (see sections l2. II 
and l2.2.T|) . It is usual to present the critical value in terms of A^. instead of p^. To avoid 
friction with the literature of absorbing phase transitions we list the corresponding 
values of A^ in Table|^ 

For the sake of completeness, we mention that the master equation of the contact 
process may be written as (see e.g. (Zj) 



Here, P{ni, . . . , n^jt) denotes the probability to find the system at time t in the config- 
uration (tt-i, n2 . . . , njv). Although a complete analytical solution of the contact process 
is still lacking, a number of rigorous results are available. For example, the transition 
is known to be continuous, upper and lower bounds on Ac exist, as well as a complete 
convergence theorem has been proven (see e.g. jl77( 11831 11H4( Hl^H] ). 

3.1.2 Domany-Kinzel automaton 

Another important 1 + 1-dimensional stochastic cellular automaton exhibiting directed 
percolation scaling behavior is the Domany-Kinzel (DK) automaton |182j . It is defined 
on a diagonal square lattice with a discrete time variable and evolves by parallel up- 
date according to the following rules (see Figure !!^ . A site at time t is occupied with 
probability p2 (pi) if both (only one) backward sites (at time t — 1) are occupied. Oth- 
erwise the site remains empty. If both backward sites are empty a spontaneous particle 
creation takes place with probability po. Similar to the contact process, the sponta- 
neous particle creation destroys the absorbing phase (empty lattice) and corresponds 
therefore to the conjugated field h. 

As pointed out above, the geometrical problem of Z) -|- 1-dimensional directed perco- 
lation corresponds to a D-dimensional dynamic process. This mapping becomes evident 
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P{ni, . . . , 1, Tii^i, . . . , riff, t) 



nN,t) 



(3.21) 
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Figure 16: The 1 + 1-dimensional Domany-Kinzel automaton. Occupied sites are marked by full circles. 
A spreading of particles is sketched in the left part. The state of the one dimensional system at time 
t + 1 is obtained by an iteration of the dynamics rules that are embodied in the conditional probabilities 
P(nt-|-i|ni t, n,._t). Thus the occupation Ut+i of a given site at time t depends on the state of the left 
(1) and right (r) backward neighbors (ni,t,?ir,t) at time t. Spontaneous particle creation corresponds to 
the conjugated field and is usually forbidden, i.e., po = h = Q. 

for the Domany-Kinzel automaton |182j where the problems of bond and site directed 
percolation are recovered for a particular choice of the conditional probabilities and 
P2- Suppose that the bonds of the square lattice of FigurelTBl are present with proba- 
bility q and that the sites of the lattice are present with probability p. This equals the 
site directed percolation (sDP) problem for (7 = 1 whereas the bond directed percolation 
(bDP) problem is recovered for p = 1. These occupation rules correspond on the other 
hand to the conditional probabilities Pi = pq and p2 = pq{2 — q) |182j . respectively. 
Thus bond directed percolation (p = 1) is recovered for 

P2 = Pi(2-pi), (3.22) 

whereas site directed percolation (q = 1) corresponds to 

P2 = Pi- (3.23) 

Analogous to the contact process (see pagel^Tj). mean field treatments present some 
insight into the critical behavior of the Domany-Kinzel automaton (see e.g. |186( I187| 



ni 71-2 




Figure 17: The mean field treatment of the 1 + 1-dimensional Domany-Kinzel automaton. In the 
simplest case one considers a lattice consisting of two sites i — 1,2. The sites are updated (n^ — » n'^) in 
parallel according to the conditional probabilities presented in Figure tT^ 
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1188(1189] ). Here, we consider an automaton consisting of two-sites i = 1,2. Depending 
on the state of the system {rii, 712} a new configuration {n[,n'^} is obtained by applying 
the dynamical rules that are embodied in the conditional probabilities po = h, pi, and 
P2. Neglecting correlations, the associated probabilities for a change of the number of 
active sites by An^ are 

p(An. = -2) = pl{l-p,f, 

p(An, = -1) = 2pIp,{1-p,) + 2p,p,{l-p,f , 

p(An, = 0) = 2pIp^^ + Ap,p,p,{l-p,) + p^Jl-hf , 

p{An,= 1) = 2p,/5,p2 ^ 2p2/i(i_/i), 

p{An^= 2) = p2/,2^ (3 24) 

where p^ denotes the density of (active) occupied sites and the density of empty sites 
is Pe = 1 — Pa- This reaction scheme leads to the differential equation 

d,p^{p„p„h) = 2{2p, - I) p^ - 2{2p,- p,)pI + 2{l- p^fh. (3.25) 

Focusing on the steady state behavior {dtP^, = 0) we find that the order parameter at 
zero field is given by 

275—1 

P.(Pi,P2,/i = 0) = V p,(pi,p2,/i = 0) = (3.26) 

2pi - P2 

The first solution is stable for p^ < 1/2 and corresponds to the inactive phase. The ac- 
tive phase is described by the second solution which is stable for pi > 1/2. Thus within 
this approximation, the phase diagram of the Domany-Kinzel automaton exhibits two 
phases separated by the critical line p^ = pi = 1/2 (see FieurelTH)) . Along this critical 
line the order parameter vanishes linearly (/3 = 1) in the active phase 

Pa(Pi,P2, h = 0) = dp + 0(5p2) (3.27) 

I-P2 

with 5p = {p — Pc)/Pc- The field dependence of the order parameter at the critical line 
is given by 

pM=Pc,P2,h) = J-^ + 0{h), (3.28) 
V 1 - P2 

yielding the field exponent a = 2. The asymptotic behavior of the equation of state 



pM,P2,h) ~ ^ + \ + a,h , (3.29) 



is obtained by performing the limits p^ 0, 6p ^ 0, and /i — > with the constraints 
that p^l\fh and p^^/Sp remain finite. The metric factors are given by = = 1/(1 — 
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Figure 18: The phase diagram of the Domany-Kinzel automaton. The left figure shows the phase 
diagram of the mean field approximation. The 1 + 1-dimensional behavior is sketched in the right figure. 
In both cases, the inactive phase (pa = 0) is separated from the active phase (0 < Pa < 1) by a line of 
second-order transitions (solid hue). Bond directed percolation corresponds to p2 = pi(2 — pi) whereas 
site directed percolation is obtained for pi = P2- For p2 = the Domany-Kinzel automaton corresponds 
to the cellular automata rule "18" of Wolfram's classification scheme |13l 11901 . Furthermore, the fully 
occupied lattice (p^ = 1) remains inactive for pi > 1/2 and P2 = 1 (thick dashed line). For P2 = 1 
the system is characterized by a particle-hole symmetry which gives rise to a different universal scaling 
behavior (marked by the full circle), so-called compact directed percolation. 

^2)- Thus the mean field behavior of the Domany-Kinzel automaton is characterized by 
the same critical exponents and despite of the metric factors a^, by the same scaling 
function as the mean field contact process [see Eq. H2.18() ] . 

It is possible to improve the mean field approximations presented above by consid- 
ering higher orders of correlations (see e.g. jl87j ). In this way, better approximations 
of the phase diagram of the 1 -|- 1-dimensional system can be obtained. But the results 
are still mean field-like, i.e., the critical exponents as well as the asymptotic scaling 
functions equal /3 = 1, o" = 2, and Eq. (|3.29|) . 

Clearly, the above derived mean field results are not valid for p2 = 1. In that case 
no particle annihilation takes place within a domain of occupied sites, i.e., the creation 
and annihilation processes are bounded to the domain walls where empty and occupied 
sites are adjacent. This corresponds to the particle-hole symmetry 

n ^ 1-n (3.30) 

which is also reflected in the mean field differential equation [Eq. (|3.25|) ] 

d,pM,p, = l,h = 0) = 4(2pi-l)p,(l-/5j. (3.31) 

Obviously, the fully occupied lattice is another absorbing state. The behavior of the 
domain wall boundaries can be mapped to a simple random walk jl82| and the domains 
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Figure 19: The order parameter and its fluctuations of the 2 + 1-dimensional generalized Domany- 
Kinzel automaton on a bcc lattice for various field values (from h = 10"'^ to h = 210""). The data 
corresponds to the site directed percolation (sDP) process. For non-zero field (lines) pa and Apa exhibit 
an analytical behavior. For zero field (symbols) the typical critical behavior of the order parameter and 
its fluctuations are obtained. The data are obtained from simulations on various system sizes L < 2048 
with periodic boundary conditions. 

of particles grow on average for > 1/2 whereas they shrink for < 1/2. Thus the 
steady state density is zero below p^ = 1/2 and = 1 above p^. At the critical 
value the order parameter p^ exhibits a jump. The associated critical exponent (3 = 
differs in all dimensions from the directed percolation values < /3dp - 1- Since domain 
branching does not take place the dynamical process for = 1 is termed compact 
directed percolation (CDP). The critical behavior equals that of the 1 + 1-dimensional 
voter model |177j and is analytically tractable due to the mapping to random walks. 
Exact results are derived for the critical exponents jl82l 1191] (/? = 0, /3' = 1, i/y = 2, 
i/x = 1, see also Table|21) as well as for certain finite-size scaling functions jl92j . In 
particular, the domain growths from a single seed exhibits a one-to-one correspondence 
to a pair of annihilating random walkers. That correspondence allows to calculate the 
percolation probability PM] [Eq. 



yielding the percolation exponent (5' = 1. In contrast to directed percolation the 
universality class of compact directed percolation is characterized by the inequality 



The number of independent critical exponents is therefore greater than for directed per- 




(3.32) 



(3.33) 
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eolation. In summary, the universality class of compact directed percolation is charac- 
terized by a continuously vanishing percolation order parameter Pp^^^ and algebraically 
diverging correlations lengths |191j . indicating a second order phase transition. But 
due to the misleading discontinuous behavior of the steady state order parameter 
the phase transition was often considered as first order. 

It is straight forward to generalize the 1 + 1-dimensional Domany-Kinzel automaton 
to higher dimensions (see e.g. |145l 11931 1142j ^. In the following, we consider cellular 
automata on a D + 1-dimensional body centered cubic (bcc) lattice where the time 
corresponds to the [0, 0, ... , 0, 1] direction. A lattice site at time t is occupied with 
probability p if at least one of its 2^ backward neighboring sites (at time t — 1) is 
occupied. Otherwise the site remains empty. This parameter choice corresponds to 
the condition probabilities Pi = P2 = ■ ■ ■ = P2D = P, i-e., site directed percolation 
is considered. In Figure^l we present simulation results for the 2 -|- 1-dimensional 
Domany-Kinzel automaton. The data show the expected critical behavior of the order 
parameter and its fluctuations. 

3.1.3 Pair contact process 

The pair contact process (POP) was introduced by Jensen |134j and is one of the sim- 
plest models with infinitely many absorbing states showing a continuous phase tran- 
sition. The process is defined on a D-dimensional cubic lattice and an asynchronous 
(random sequential) update scheme is applied. A lattice site may be either occupied 
(n = 1) or empty (n = 0). Pairs of adjacent occupied sites, linked by an active bond, 
annihilate each other with rate p otherwise an offspring is created at a neighboring 
site provided the target site is empty (see Fieure OU)) . The density of active bonds 
Pa is the order parameter of a continuous phase transition from an active state to an 
inactive absorbing state without particle pairs. Similar to the contact process and to 
the Domany-Kinzel automaton a spontaneous particle creation acts as a conjugated 
field |14Uj . The behavior of the PCP order parameter and its fluctuations are plotted 
in Figure On 

Since isolated particles remain inactive, any conflguration containing only isolated 
particles is absorbing. In case of the 1 + 1-dimensional pair contact process with L sites 
and periodic boundary conditions the number of absorbing states is asymptotically 
given by the golden mean |194j 



In the thermodynamic limit the pair contact process is characterized by infinitely many 
absorbing states. Due to that non-unique absorbing phase the universality hypothe- 




(3.34) 
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Figure 20: The 1 + 1-dimensional pair contact process |131| . Lattice sites may be either empty (open 
circles) or occupied by a particle (full circles). Pairs of occupied sites (solid bonds) are considered 
as active whereas isolated particles remain inactive. A pair is annihilated with rate p, otherwise an 
offspring is created at an empty neighboring site selected at random. 

sis of Janssen and Grassberger can not be applied. Therefore, the critical behav- 
ior of the pair contact process was intensively investigated by simulations (includ- 
ing |131l 11351 11951 11961 1197j ) . It was shown numerically that the critical scaling behav- 
ior of the 1 + 1-dimensional pair contact process is characterized by the same critical 
exponents |134| 1135j as well as by the same universal scaling functions as directed per- 
colation jl4Uj . In particular the latter result provides a convincing identification of 
the universal behavior. Thus despite the different structure of the absorbing phase the 
1 + 1-dimensional pair contact process belongs to the directed percolation universality 
class. This numerical evidence confirms a corresponding renormalization group analysis 
predicting DP universal behavior |198j in all dimensions. But the scaling behavior of 
the PCP in higher dimension is still a matter of controversial discussions. A recently 

0.10 [ — , — , — , — , — , — , — , — , — , — , — , — , — , 1 




Figure 21: The order parameter and its fluctuations (inset) of the 2 + 1-dimensional pair contact process 
on a square lattice for various values of the field (from h — 10"'' to ft = 6 10"^). The symbols mark 
zero-held behavior. The data are obtained from simulations on various system sizes L < 2048 with 
periodic boundary conditions. 
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performed renormalization group analysis conjectures that the pair contact process ex- 
hibits a dynamical percolation-like scaling behavior 119911201)] . A dynamical percolation 
cluster at criticality equals a fractal cluster of ordinary percolation on the same lattice. 
Thus, the dynamical percolation universality class |2()H \2()2\ VliYAi I2()4j differs from the 
directed percolation universality class. In particular the upper critical dimension equals 

= 6 instead of -De = 4 for DP. So far the investigations of the PCP are limited to 
the 1-1-1- and 2 -|- 1-dimensional |2U5j systems. Only recently performed numerical sim- 
ulations |14,Sj provide estimates of the scaling behavior in higher dimensions, including 
the exponents as well as the scaling functions. The results are presented in this work. 

A modification of the pair contact process, the so-called pair contact process with 
diffusion (PCPD), attracted a lot of attention in the last years (e.g. |H OUBl IWl 

una eohi mni Em ussi En Eia eei ehi eei eisi Em eihi). in contrast to the 

original PCP isolated particles are allowed to diffuse in case of the PCPD. Clearly, 
the diffusion of isolated particles changes the structure of the absorbing phase. It only 
contains the empty lattice and a state with a single endlessly wandering particle. The 
question of the universal scaling behavior of the pair contact process with diffusion 
is still a matter of controversial discussions in the literature (see |219j for a recent 
review). For example, non-universal continuously varying exponents [2111 1155] . a new 
universality class |208ll213ll214j as well as a crossover to the universality class of directed 
percolation j215j are reported. The striking variety of the observed conflicting scaling 
scenarios reflects the non-trivial character of the pair contact process with diffusion 
which remains a topic for future research. 

3.1.4 Other models 

In this section we briefly consider two lattice models exhibiting a directed percolation- 
like phase transition. The first one is the threshold transfer process (TTP) which 
was introduced in |136j . Here, lattice sites may be empty (n = 0), occupied by one 
particle (n = 1), or occupied by two particles (n = 2). Double occupied lattice sites 
are considered as active. In that case both particles may move to the left (1) and right 
(r) neighbor if possible, i.e.. 

Hi — > rii -|- 1 if ni < 2 , 

— > + I if < 2 . (3.35) 

Additionally to the particle movement, creation and annihilation processes are incorpo- 
rated. A particle is created at an empty lattice site (0 — > 1) with probability r whereas 
a particle annihilation (1 — > 0) takes place with probability 1 — r. In the absence of 
double occupied sites the dynamics is characterized by a fluctuating steady state with 
a density r of single occupied sites. The density of double occupied sites is identified 
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as the order parameter of the process, and any configuration devoid of double occu- 
pied sites is absorbing. The probability r controls the particle density, and a non-zero 
density of active sites occurs only for sufficiently large values of r. In contrast to the in- 
finitely many frozen absorbing configurations of the pair contact process, the threshold 
transfer process is characterized by fluctuating absorbing states. Nevertheless steady 
state numerical simulations of the 1 + 1-dimensional threshold transfer process yield 
critical exponents that are in agreement with the corresponding DP values |136j . So 
far no systematic analysis of the TTP in higher dimensions was performed. 

In this work we focus on the 2 + 1-dimensional model. In that case one tries to 
transfer both particles of a given active site to randomly chosen empty or single occupied 
nearest neighbor sites. Furthermore we apply an external field that is conjugated to the 
order parameter. In contrast to the models discussed above the conjugated field can not 
be implemented by a particle creation. A particle creation with rate h would affect the 
particle density, i.e., the control parameter of the phase transition. But the conjugated 
field has to be independent of the control parameter. Therefore we implement the 
conjugated field by a diffusion-like field that acts by particle movements. A particle 
on a given lattice site moves to a randomly selected neighbor with probability h, if 
n < 2. Thus the conjugated field of the TTP differs from the conjugated field of the 
Domany-Kinzel automaton, the contact process, and the pair contact process. 

Another model exhibiting a directed percolation-like absorbing phase transition 
is the Ziff-Gulari-Barshad (ZGB) model [5^]. This model mimics the heterogeneous 
catalysis of the carbon monoxide oxidation 

2C0 + O2 — > 2CO2 (3.36) 

on a catalytic material, e.g. platinum. The catalytic surface is represented by a square 
lattice where CO or O2 can be adsorbed from a gas phase with concentration y for 
carbon monoxide and 1 — y for oxygen, respectively. The concentration y is the control 
parameter of the model determining the density of the components on the catalytic 
surface. Adsorbed oxygen molecules dissociate at the catalytic surface into pairs of 
O atoms. It is assumed that the lattice sites are either empty, occupied by an CO 
molecule, or occupied by an O atom. Adjacent CO and O react instantaneously and 
the resulting CO2 molecule leaves the system. Obviously, the system is trapped in 
absorbing configurations if the lattice is completely covered by carbon monoxide or 
completely covered by oxygen. The dynamics of the system is attracted by these 
absorbing configurations for sufficiently large CO concentrations and for sufficiently 
large O2 concentrations. Numerical simulations show that catalytic activity occurs 
in the range 0.390 < y < 0.525 |146j only (see Figure l22|) . The system undergoes a 
second order phase transition to the oxygen passivated phase whereas a first order 
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Figure 22: The schematic phase diagram of the Ziff-Gulari-Barshad model for the catalytic carbon- 
monoxide oxidation. The steady state concentrations of CO and O2 on the catalytic surface are plotted 
as function of the CO concentration y in the gas phase. The system undergoes a second order phase 
transition to the oxygen passivated phase whereas a first order phase transition takes place if the 
CO passivated phase is approached. The right figure shows the density of vacant lattice sites and its 
fluctuations. 

phase transition takes place if the CO passivated phase is approached. In particular, the 
continuous phase transition belongs to the universality class of directed percolation |22n| 
I146U151] . This might be surprising at a first glance since the ZGB model is characterized 
by two distinct chemical components, CO and O. But the catalytic activity is connected 
to the density of vacant sites, i.e., to a single component order parameter j22nj . Thus 
the observed directed percolation behavior is in full agreement with the universality 
hypothesis of Janssen and Grassberger. 

But one has to stress that the ZGB model is an oversimplified representation of the 
catalytic carbon monoxide oxidation. A more realistic modeling has to incorporate for 
example mobility and desorption processes as well as a repulsive interaction between 
adsorbed oxygen molecules (see e.g. pnil221| '). These features affect the critical behavior 
and drive the system out of the directed percolation universality class. 

Eventually we just mention that depinning phenomena can be mapped to directed 
percolation. Whereas certain interface depinning phenomena are related to 1 + 1- 
dimensional DP, higher dimensional behavior is related to line pinning phenomena. 
Instead of going into the details we refer the interested readers to the reviews [Tj |H1 110] ■ 

3.2 Steady state scaling behavior 

In the following we consider the steady state scaling behavior of different models in 
various dimensions. We focus our attention on the order parameter, i.e., the equa- 
tion of state, its fiuctuations as well as on the order parameter susceptibility. First 
we consider the scaling behavior below the upper critical dimension and compare the 
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results to those of an renormalization group approach. Then we show that the scaling 
behavior above the upper critical dimension is in full agreement with the mean field 
predictions. The scaling scenario at the upper critical dimension is characterized by 
logarithmic corrections, and will be discussed in detail. Furthermore, certain universal 
amplitude combinations are considered which are related to the order parameter and 
its fluctuations. Analogous to equilibrium phase transitions, these amplitude combina- 
tions are directly related to particular values of the universal scaling functions and are 
of great experimental interest. A finite-size scaling analysis of the order parameter, its 
fluctuations as well as a fourth order cumulant is presented at the end of this section. 

3.2.1 Below the upper critical dimension 

Similar to equilibrium phase transitions we assume that sufficiently close to the critical 
point the order parameter p^j its fluctuations Ap^, and the order parameter suscepti- 
bility X can be described by generalized homogeneous functions (see Eq. (|1.28j) ) 

p,{dp, h) ~ A"'^ R{a^dp A, a,h X") , (3.37) 

a^Ap,{5p,h) ~ X'^' D{ap6pX,a^hX''), (3.38) 

axX{^,h) ~ X"' X{a,6p X,a,,hX''). (3.39) 

The fluctuation exponent fulfills the hyperscaling law [see Eq. ()A.6|) ] 

y = v^D - 2p. (3.40) 

Note that in contrast to the scaling behavior of equilibrium phase transitions the non- 
equilibrium absorbing phase transition of directed percolation is characterized by 7 7^ 
7'. Furthermore, the spatial and temporal correlation length are expected to scale as 

a^^{dp,h) ~ A-'^^ H^(a>A,a,/iA'^), (3.41) 

0,1^11 (5p,/i) ~ A""ii S||(a/pA,a,/iA'^). (3.42) 

As usual, h denotes the conjugated field and 5p denotes the distance to the critical 
point, e.g. (5p = (A — X^)/Xc for the contact process, dp = {jp — Pc)/Pc for site directed 
percolation and for the pair contact process, 6p = {r — r^)/r^ for the threshold transfer 
process, etc. Additionally to the above Hankey-Stanley scaling forms one may consider 
the corresponding Widom- Griffiths scaling forms, e.g. 

a^h{dp,p:) ~ A-" H{a,5pX,p^X^). (3.43) 

Once the non-universal metric factors are chosen in a specified way, the universal 
scaling functions -R, i?, D, X, and are the same for all systems within a given 
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Figure 23: The determination of the non-universal metric factors Op and for site directed percolation 
(sDP), the contact process (CP), and the pair contact process (PCP) in D = 2. The left figure shows 
the order parameter at zero-field {h = 0) whereas the right figure displays the order parameter at the 
critical value {5p = 0). Plotting pl^^ vs. apSp and p^^^ vs. Uhh the curves tend asymptotically to the 
function f{x) = x (dashed lines) when the transition point {5p — >■ and ft ^ 0) is approached. 

universality class. All non-universal system-dependent features (e.g. the lattice struc- 
ture, the range of interaction, as long as the interaction decreases sufficient rapidly 
as function of separation, etc.) are contained in the non-universal metric factors Op, 
dh, , • • ■ |M] • The above scaling forms are valid for D ^ D^. At the upper critical 
dimension they have to be modified by logarithmic corrections jl41j . 

Note that all universal functions are characterized by the same arguments. Identical 
scaling arguments arise naturally in equilibrium thermodynamics where in principle all 
functions can be derived from a single thermodynamic potential, e.g. the free energy (see 
11.2)1 . In case of non-equilibrium phase transitions one can argue that the functions can 
be derived from a corresponding Langevin equation. Furthermore, beyond the equation 
of state [Eq. 1)3. 37() ] , one new metric factor is introduced for each additional quantity 
being under consideration, here a^, a^, a^, and ay. But this phenomenological ansatz 
does not imply that the metric factors are independent of each other. For example only 
two metric factors are needed to describe the scaling behavior in equilibrium |222j . 
often termed two-scale-factor-universality. 

Similar to equilibrium [Eq. H1.29() ] we norm the universal scaling functions by 

E(1,0) = 1, ^(0,1) = 1, ^(0,1) = 1. (3.44) 

The first two normalizations imply 

H{l,l) = 0, -^'(O,!) = 1. (3.45) 

In that way the non- universal metric factors a^, a^, and are determined by the 
amplitudes of the power-laws 

p.(5p,/i = 0) ~ {%5pf, (3.46) 
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p^{dp = Q,h) ~ {a,^hf'% (3.47) 
Ap,(dp = 0, /i) ~ (a,/i)-^'/". (3.48) 

The determination of the metric factors and a,^ for site directed percolation, the 
contact process, and the pair contact process is shown in Figure !^ For this purpose 
the best known estimates for the critical exponents, as given in TableQ are used. The 
values of the non-universal metric factors are listed in |142| I143j . 

Taking into consideration that the susceptibility is defined as the derivative of the 
order parameter with respect to the conjugated field [Eq. (|2.14() ] we find 

X(3;,y) = dyR{x,y), = , (3.49) 

as well as the scaling law 

7 = a - (5. (3.50) 

This scaling law corresponds to the well known Widom law [Eq. (|1.9() ] of equilibrium 
phase transitions. Furthermore, comparing Eq. (|3.39)) for (5p = to the definition of the 
susceptibility 

a,x(5p,/i) ~ (a./i)-^/"X(0,l), X = d^p^ = d^ia^hf'^ (3.51) 

leads to the condition 

X(0,1) = ^. (3.52) 
a 

The above result offers a useful consistency check of the numerical estimates of the 
susceptibility. It is worth mentioning that the validity of the scaling form Eq. (|3.37|) 
implies the required singularity of the susceptibility [Eq. H2.14() ] , confirming that the 
applied external field is conjugated to the order parameter. 

Choosing a^hX" = 1 \n Eqs. (|3. 37113. 39|) we obtain the scaling forms 



p.{dp,h) ~ (a, /i)^/'^ i?(aXa, /i)-i/^ 1) , (3.53) 

a^^p^{5p,h) ~ {a^h)-^'/'' b{a,5p{a^h)-^''',l), (3.54) 

a^x{5p,h) ~ {a^h)-^'^ X{a,dp{a,^h)-^'%1). (3.55) 
Thus plotting the rescaled quantities 



p^{5p, h) {a, /i)"^/'^ , Ap^{5p, h) (a, /i)^'/- ^ xi^P, h) (a, h)-^/^ (3.56) 

as a function of the rescaled control parameter Updp (o^/i)"^'''^ the corresponding data of 
all systems belonging to the same universality class have to collapse onto the universal 
curves R{x, 1), D{x, 1), and X(x, 1). This is shown for R{x, 1) in Figure l2H where the 
rescaled order parameter is plotted for various two-dimensional models. As can be 
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Figure 24: The universal scaling function R{x, 1) of the directed percolation universality class (left). The 
data are plotted according to Eq. 13.531 1. All models considered are characterized by the same universal 
scaling function, an impressive manifestation of the robustness of the directed percolation universality 
class with respect to variations of the microscopic interactions. Neglecting the non-universal metric 
factors each model is characterized by its own scaling function (right). For all models considered 
the scaling plots contain at least four different curves corresponding to four different field values (see 
e.g. Figure lT^ and Figure l^H . The circles mark the condition -R(0, 1) — 1. 



seen the data of all considered models collapse onto the same scaling function, clearly 
supporting the assumption that all models belong to the universality class of directed 
percolation. Furthermore, the data-collapse confirms the accuracy of the numerically 
estimated values |151j of the critical exponents. 

The universal scaling function R{x, 1) is displayed in Figure 031 for D = 1, 2, 3. For 
each dimension, the data of the site directed percolation process, of the contact process 
as well as of the pair contact process collapse onto the unique scaling function. As 
expected the scaling functions vary with the spatial dimensionality below D^. 

Similar to equilibrium the powerful and versatile e-expansion provides estimates 
for almost all quantities of interest (see pagel47|). in particular for the critical expo- 
nents EHl ESI CEOl and for the equation of state [Ml HIS]- Using a paramet- 
ric representation [2241 I225j of the directed percolation phase transition Janssen et 
al. showed that the equation of state is given within the e^-approximation by the re- 
markable simple scaling function |181j 

H{e) = e (2 - 9) + 0{e^) . (3.57) 

Here the scaling behavior of the quantities p^, 5p, and h is transformed to the variables 
R and 6 through the relations 



b6p = R{l-e), 




(3.58) 
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Figure 25: The universal scaling function R(x, 1) of the directed percolation universality class in various 
dimensions. The two- and three-dimensional data are vertically shifted by a factor 8 and 64 in order to 
avoid overlaps. The circles mark the condition -R(0, 1) = 1. The gray lines correspond to an e-expansion 
of a renormalization group approach 11811 . 

The equation of state is given by 




(3.59) 



with the metric factors a and b. The whole phase diagram is described by the parameter 
range R> and 6 G [0, 2]. Compared to the Widom-Griffiths scaling form [Eq. (|3.17j) ] 
the parametric representation is not restricted to sufficiently small scaling arguments in 
the absorbing phase. The scaling function R{x, 1) can be calculated from the equation 
above, and the results are compared to the numerical data in Figure 031 As expected, 
the differences increase with increasing perturbation parameter e, i.e., with decreasing 
dimension and are especially strong in D = I |14nj . Thus higher orders than 0{e'^) are 
necessary to describe the scaling behavior of directed percolation in low dimensions. 
A more detailed discussion of the accuracy of the renormalization group data will be 
presented in 13.2.41 

The universal scaling functions of the order parameter fluctuations and the order 
parameter susceptibility are shown in FigureUHl The susceptibility is obtained by per- 
forming the numerical derivative of the order parameter with respect to the conjugated 
field. The perfect data-collapses confirm the scaling forms Eqs. (|3. 38113. 3'9|) . All scaling 
functions exhibit for D = 1,2,3 a clear maximum signalling the divergence of Ap^ and 
X at the critical point. The susceptibility data fulfills Eq. (|3.52|) . reflecting the accuracy 
of the performed analysis. In summary all considered models are characterized by the 
same critical exponents and the same steady state scaling functions R, D, and X. Thus 
the models are characterized by the same scaling behavior for D = 1,2,3. 
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Figure 26: The universal scaling functions of the fluctuations D{x, 1) (left) and of the susceptibility 
X(x, 1) (right) in various dimensions. The two- and three-dimensional fluctuation data are vertically 
shifted by a factor 8 and 64 in order to avoid overlaps. The circles mark the condition D{0, 1) = 1. In 
case of the susceptibility the two- and three-dimensional data are vertically shifted by a factor 3/2 and 
9/4 in order to avoid overlaps. The circles mark the condition X(0, 1) = P/a and reflect the accuracy 
of the performed analysis. 



3.2.2 Above the upper critical dimension 

According to the renormalization group scenario the stable fixed point of the recursion 
relations for D > is usually the trivial fixed point with classical, i.e., mean field 
universal quantities. Thus, in contrast to the situation below the critical exponents 
as well as the universal scaling functions are independent of the particular value of 
the dimension for D > D^. In most cases it is possible to derive these mean field 
exponents and even the scaling functions exactly since correlations and fluctuations can 
be neglected above D^. The mean field behavior of the contact process and the Domany- 
Kinzel automaton was considered in section 12.11 and in section 13.1.21 respectively. A 
corresponding analysis of the pair contact process can be found in |14Uj . The obtained 
universal scaling functions are given by (see e.g. |116j ) 




i?MF(x,y) = 7T + \ y + o ' (3-60) 



DyiT,{x,y) = = , (3.61) 
'y + {x/2f 

^MF{x,y) = , ^ (3.62) 

2 ^y + {x/2f 

- ,MF{x,y) = I ^ = (3.63) 
■ y + {x/2f 
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Figure 27: The universal scaling function of the order parameter R{x, 1) (left) and the fluctuations 
D{x, 1) (right) for D = 5. The five-dimensional data agree with the corresponding mean field functions 
Eq. IpTfnjl and Eq. IpTTT]! . respectively. 



^\\,MFix,y) = I = (3.64) 

\Jy + {^m 

yieldling the mean field exponents /3mf = l, ^mf = 2, 7mf = 1, 7mf = 0, i^x,mf = 1/2, 
and i^[|,MF = 1- These functions obey the scaling forms Eqs. H3.53I3.55|) . 

The scaling behavior of the fluctuations deserves comment. The exponent 7[jp = 
corresponds to a jump of the fluctuations and the scaling form Eq. (|3.3<S|) reduces to 



a/^Ap^{5p,h) ~ byij,{a^5p X^a^hX") . (3.65) 
Using again Z)mf(0, 1) = 1, the non-universal metric factor Ca is determined by 

Apa(<5p = 0, /i) 

The reduced scaling form Eq. 1)3.65(1 implies some interesting properties of the universal 
function Dmf- The value i)MF(l,0) is now given by 

n n n\ Ap,{6p,h = 0) 
= Ap^{6p = 0,h) ' 

and the mean field scaling function fulfills the symmetries 

Dmf{1,x) = £'mf(x-^/",1), (3.68) 
DmAx,1) = L)mf(1,x-"), (3.69) 

for X > 0. In particular we obtain for x ^ Dmf(1>0) = Dmf(oo, 1) and Dmf(0, 1) = 
£)mf(1) oo), respectively. This symmetry is an artefact of the zero fluctuation exponent 
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Figure 28: The universal Widorn- Griffiths scaling function H{x, 1) for various dimensions as a function 
of a:: = app~^^^ . The numerically obtained curves H{x, 1) tend to the mean field result (long dashed 
lines) Hmf{x, 1) — 1 — x with increasing dimension. 



and does not occur below the upper critical dimension. Therefore, it can be used to 
check whether a given D-dimensional system exhibits mean field behavior. 

Numerical data of the five-dimensional models are presented in Figure OTl A good 
data-collapse of the rescaled data with the universal mean field scaling functions 

RmAxA) = I + \/l + (f)' (3-70) 

5mf(x,1) = 1 + , (3.71) 

2 V'l + {x/2f 

is obtained. The observed mean field behavior for D = 5 agrees with the renormaliza- 
tion group result D, = 4 fFM\Tn^ . 

As already pointed out the scaling functions depend on the spatial dimension be- 
low D^. In particular the Widom- Griffiths scaling form [Eq. ()3.43|) ] is well suited to 
examine how the universal functions approach the mean field solution with increasing 
dimension. The mean field Widom-Griffiths scaling form is given by 

HMF{x,y) = - xy. (3.72) 

The rescaled field a^h p~"l^ is plotted in Figure OSl as a function of the rescaled control 
parameter a^bp p~'^l^ . As can be seen the scaling functions ii{x^ 1) tend to the mean 
field curve Hmf{x, 1) = 1 — x with increasing dimension. The observed agreement 
of the data of the five-dimensional models with the mean field scaling function Rmf, 
Hmf, and Dmf leads to the result that the upper critical dimension is less than five. 
This is a non-trivial result. The value of is often predicted by field theoretical 
treatments, and a reliable determination of the upper critical dimension is therefore 
of particular interest. For example, two contrary field theories conjecture for the pair 



3.2 Steady state scaling behavior 



67 



contact process = A |198j and = 6 |199| 1200) , respectively. The latter result is in 
clear contradiction to the observed mean field behavior in D = 5. In this way, an upper 
bound on can be obtained just by checking whether the numerical or experimental 
data are consistent with the usually known universal mean scaling functions. 

3.2.3 At the upper critical dimension 

The scaling behavior around the upper critical dimension is well understood within 
the renormalization group theory (see for instance |22l EHl CZI ) • For D > the stable 
fixed point of the corresponding renormalization equations is usually a trivial fixed 
point with mean field exponents. This trivial fixed point is unstable for D < D^, and 
a different stable fixed point exists with nonclassical exponents. For D = both 
fixed points are identical and marginally stable. In this case, the scaling behavior 
is characterized by mean field power-laws modified by logarithmic corrections. For 
instance the order parameter is expected to obey in leading order 

p^{6p,h = 0) 6c 6p\ln6p\^, (3.73) 
p^{6p = 0,h) 6c \^\lnh\^. (3.74) 

Greek capitals will be used in the following to denote the logarithmic correction expo- 
nents. These logarithmic correction exponents are characteristic features of the whole 
universality class similar to the usual critical exponents. Numerous theoretical, nu- 
merical, as well as experimental investigations of critical systems at the upper critical 
dimension have been performed (see for example ^7711^ 1^ 0!M OITTl 1^ O!^! 
ESSl Enni ESa OSHI EMI CHl IHa E2S1 Em Eia Eini). Logarithmic corrections 
make the data analysis notoriously difficult. Hence, most investigations focus on the 
determination of the correction exponents only, lacking the determination of the scaling 
functions at . Recent numerical jl411 1242j as well as analytical results j223j reveal 
that the concept of universal scaling functions can also be applied to the upper critical 
dimension. 

A method of analysis to determine both the universal scaling functions and the cor- 
rection exponents at the upper critical dimension was developed in |141j and successfully 
applied to site directed percolation jl42j . In this work the following phenomenological 
scaling ansatz is used for the Hankel-Stanley form of the equation of state (all terms 
in leading order) 

a,p.{6p,h) ~ A^'^^'iF llnAI"^ ^(a/p A | In A|", 0;,/i A^^^^ | In A|") . (3.75) 

According to this scaling form, the order parameter at zero field {h = 0) and at the 
critical density {5p = 0) is given by 

a^p^{6p,h = 0) ~ ap5p \lnap5p\^ -R(1,0), (3.76) 
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a,p.{5p = 0,h) ~ llnv^l^ ^(0,1) (3.77) 

with B = n + A and S = H/2 + A. Similar to D D, the normahzation ^(0, 1) = 
R{1,0) = 1 is used. Furthermore, the scaling behavior of the equation of state is given 
in leading order by 

a^p^{6p,h) ~ \/aJi \ In y/aji\^ R{x,l) (3.78) 

where x denotes the scaling argument 

X = ap6p \/ a^h ^ | In \/ a^h |* (3.79) 

with^' = n-H/2 = B-S. 

In case of directed percolation it is possible to confirm the scaling form Eq. (II175I1 
by a renormalization group analysis |223j . yielding in addition the values A = 7/12, 
n = —1/4, and H = —1/2. Thus the scaling behavior of the equation of state is 
determined by 

B = S = 1/3, ^' = 0. (3.80) 

It is worth mentioning that in contrast to the renormalization group results below 
the upper critical dimension, the logarithmic correction exponents do not rely on ap- 
proximation schemes like e- or 1/n-expansions. Within the RG theory they are exact 
results. 

Similarly to the order parameter the following scaling form is used for the fluctua- 
tions [m] 

Ap,((5p,/i) ~ A^MF|lnA|^ Z)(a/p A | In A]", a;,/i A""" | In A|") . (3.81) 

Using the mean field value = and taking into account that numerical simulations 
show jl42j that the fluctuations remain flnite at the critical point (i.e. K = 0) the 
scaling form 

a^/\p.X^,h) ~ D{x,l) (3.82) 

is obtained, where the scaling argument x is given by Eq. ()3.79|) with ^' = 0. The 
non- universal metric factor Oa is determined again by the condition 15(0, 1) = 1. 

In that way, the scaling behavior of the order parameter and its fluctuations at 
is characterized by two exponents (B = 1/3 and S = 1/3) and four unknown non- 
universal metric factors (a^, Op, a^, Oa). Following jl41j these values are determined 
by several conditions which are applied simultaneously: flrst, both the rescaled order 
parameter and the rescaled fluctuation data have to collapse to the universal functions 
R{x,l) and L'(x,l). Second, the order parameter behavior at zero fleld and at the 
critical density is asymptotically given by the simple function f{x) = x when plotting 
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Figure 29: The universal scaling functions of the order parameter (upper left) and its fluctuations 
(upper right) at the upper critical dimension = 4. The logarithmic correction exponents are given 
by B = S = 1/3 |228| and 'I' = 0. For both models considered, the scaling plots contain at least three 
different curves corresponding to three different field values. The circles mark the condition -R(0, 1) = 1 
and -D(0, 1) = 1, respectively. The lower figures show the order parameter for zero-field (left) and at the 
critical control parameter (right), respectively. The order parameter is rescaled according to Eq. 13.7611 
and Eq. 13.771 . Approaching the transition point (/i — > and 5p 0) the data tend to the function 
j(x) = X (dashed lines) as required (see text). 



[aaPa(<5p, 0)/ap(5p]^/^ vs. |lnap(5p| and \a^pj^^ h) / a^hY/^ vs. | In y/a^h |, respectively. 
Applying this analysis to data of both the contact process and site directed percolation, 
convincing results are obtained for B = S = l/3, ^ = 0, and for the values of the non- 
universal metric factors listed in |1421 1143] . The corresponding plots are presented in 
Figure|2ni and show that the concept of universal scaling functions can be applied to 
the upper critical dimension. 

Notice that no data-collapse is obtained if logarithmic corrections are neglected, 
i.e., for B = S = 0. Thus, the leading logarithmic corrections have to be taken into 
account in order to study steady state scaling functions. It is therefore remarkable that 
recently performed off-lattice simulations of the dynamical scaling behavior at D^. = 4 
reveals that logarithmic corrections of higher orders (e.g. 0(\n\nt)) are necessary to 
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describe the numerical data |244j . Although the steady state results presented here 
are quite convincing, we expect that even better results are obtained by incorporation 
higher order corrections. But in contrast to the dynamical scaling behavior, no analyt- 
ical results of higher logarithmic corrections for the steady state scaling behavior are 
available so far. 



3.2.4 Universal amplitude combinations 

Additionally to the critical exponents and universal scaling functions, a universality 
class is also characterized by various universal amplitude combinations [681 1245j . These 
amplitude combinations emerge from the universality of the scaling functions since uni- 
versal amplitude combinations correspond to particular values of the scaling functions. 
It is known from critical equilibrium systems that universal amplitude combinations 
vary usually more widely than the critical exponents. Thus these amplitude combina- 
tions are very useful in order to identify the universality class of a phase transition or 
to provide additional and often convincing evidence for the universal behavior. Fur- 
thermore, the measurement of amplitude combinations in experiments or simulations 
yields a reliable test for theoretical predictions. In particular, approximation schemes 
of the renormalization group, such as e- or 1/n-expansions, are widely used to obtain 
explicit and systematic estimates of the amplitude combinations (see section ll.4|l . 

Usually numerical investigations focus on amplitude combinations arising from 
finite-size scaling analyses (see e.g. |245j ). A well known example is the value of Binder's 
fourth order cumulant at criticality (see e.g. |112j ). Instead of those finite-size proper- 
ties we focus our attention on bulk critical behavior since bulk amplitude combinations 
are of great experimental interest. Furthermore, they can be compared to the renor- 
malization group results of Janssen et al. I181j . 

The amplitude combination considered first is related to the susceptibility behavior 
below and above the transition. In analogy to equilibrium (see the susceptibility 
diverges as 

X(<5p > 0, /i = 0) ~ a,,+ <5p-^ (3.83) 
X(<5p < 0, /i = 0) ~ a,,.{-6pr\ (3.84) 

if the critical point is approached from above and below, respectively. Using the scaling 
form Eq. (|3.39|) the susceptibility ratio 



x{5p>^,h) _ X( ap5p X, tthh X'^) 
xiSp < 0,h) ~ X{-a,5pX,a,hX'') 



B±ll^ (3.85) 
X(-l,x) 



is clearly a universal quantity for all values of the scaling argument x = a,^h\ap5p\ " . 
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Figure 30: The universal scaling function X(l, a;)/X(— 1, a;) for various dimensions and models. Site 
directed percolation data are marked by closed symbols. Open symbols correspond to data of the pair 
contact process for D = 2 and to data of the contact process for D = 3 and _D = 5, respectively. The 
dashed lines correspond to an e-expansion of a renormalization group approach |181| . The universal 
amplitude X(l, 0)/X(— 1, 0) is obtained from the extrapolation ahh\ap Spl^"^ 0. 



In particular it equals the ratio a^_+/a^__ for vanishing field 



ix.+ V X(+l,x) X(+1,0) 



lim -^^ ^ = -^^ ^ . (3.86) 

h^a X(-l,x) X(-1,0) 



In case of the mean field approximation Eq. (|3.62j) leads to 

XMF(+l,a^) 



Xmf( 



1 (3.87) 



for all X. The susceptibility ratio X(+l, x)/X(— 1, a:) is shown in Figure l3fll for various 
models and different dimensions. The five dimensional data agree well with the mean 
field prediction. For D < the ratio reflects the crossover from mean field to non- 
mean field behavior. Far away from the transition point, the critical fluctuations are 
suppressed and the behavior of the system is well described by the mean field solution. 
Approaching criticality, the critical fluctuations increase and a crossover to the D- 
dimensional behavior takes place. 

A parametric representation of the susceptibility and therefore of the susceptibility 
ratio X(-|-l, x)/X(— 1, x) is derived in jl81j . It is instructive to compare these results to 
the numerical data (see Fig.l30|) since the theoretical curve reflects the accuracy of the 
RG estimations of all involved quantities, the exponent, the scaling functions, as well 
as the non-universal metric factors |142j . The perfect agreement between the numerical 
data and the RG curve for D = 3 indicates that all quantities are approximated well. 
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In the two-dimensional case a horizontal shift is observed between the numerical data 
and the RG-estimates. Thus the RG-approach yields good estimates for the exponents 
and the scaling functions, but the metric factors are of significantly less quality. For 
D = 1 the e^-approximation does not provide appropriate estimates of the DP scaling 
behavior similar to the equation of state (see Figurel25|) . 

The consideration of the zero-field universal amplitude ratio X(+l, 0)/X(— 1, 0) of- 
fers a more quantitative check of the renormalization group results. Numerically this 
ratio is obtained by an extrapolation of the susceptibility ratio X(+l, x)/X(— 1, to 
X ^ 0. The estimated values are listed in Tabled Within the e^-expansion this ampli- 
tude ratio is given by 



X(+1,0) ^ ^ _ ^ 
X(-1,0) 3 



V 288 144 3 



(3.^ 



Despite of the negative and therefore unphysical results for D = 1 the RG results agree 
well with the numerical data (see also Fig. l3U() . For example, the two-dimensional values 
differ only by 3%. That has to be compared to the difference of the critical exponents. 
For example the RG estimate of the order parameter exponent [Eq. (|3.14j) ] differs for 
D = 2 from the best known numerical value by 6% (see Table^). The similar accuracy 
of the critical exponent and of the amplitude combination is contrary to an observation 
from critical equilibrium systems where the exponents are usually calculated more 
accurately by e-expansions than universal amplitude combinations (see section IL^ for 
an example, as well as ^ possible explanation is that the e^-approximation yields 

for e = 1,2 much better results for directed percolation than e.g. for the Ising model. 
This is also reflected by the e- approximations of the equation of state. The various forms 
of the equation of state are less complicated for directed percolation. In particular the 
parametric representation of the equation of state [Eq. 1)3. 57() ] is remarkably simple. 

Analogous to the susceptibility the universal amplitude ratio of the fluctuations is 
given by 

Ap^{5p > 0, h) ^ D{+l,x) 

Ap,{6p<0,h) D{-l,x) ^ ■ ' 

with X = tthhlttpSpl^'^ . In case of absorbing phase transitions this ratio diverges for 
vanishing field. For 6p < the order parameter fluctuations are zero for h = (absorb- 
ing state) whereas the fluctuations remain finite above the transition {6p > 0). Thus 
absorbing phase transitions are generally characterized by 

^ oo, (3,90) 

D(-1,0) 

In Figure l3T] the fluctuation ratio is plotted as a function of the scaling variable x 
for various dimensions and models. The fluctuation ratios diverge for x ^ in all 
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Figure 31: The universal scaling function D(l, a::)/_D(— 1, a;) for various dimensions. Site directed per- 
colation data are marked by closed symbols whereas open symbols correspond to data of the contact 
process. The dashed line corresponds to the mean field scaling behavior. 



dimensions. The one-dimensional system exhibits a particular behavior characterized 
by the minimum of the corresponding curve. The origin of this behavior is that for 
D = 2,3 the universal scaling function D(x, 1) exhibits a maximum for x > 0, whereas 
for L) = 1 it is located at x < (see FigureEHl)- In case of the five-dimensional data a 
perfect agreement is observed with the mean field behavior 

Dmf{+1,x) ^ 1 + yi + 4x _^ l + 2x 
DMF{-i,x) - 1 + VI + 4x ^0 2x ' 

Finally we consider the universal amplitude combination Ry. (see Eq. (|1.49|) ) which 
can be defined in the notation of absorbing phase transitions by 

= a„+ a-' a^-^. (3.92) 

The scaling forms Eq. p.37|l and Eq. (|3.39|) lead to 

R^ = X(1,0) (3.93) 

which is obviously a universal quantity. In Figure[S21 the scaling function X(l,a;) is 
plotted as a function of x = a^h {a^ 5p)~'' for various dimensions below D^. The 
corresponding data saturates for x ^ 0, and the obtained estimates are listed in Tabled 
Compared to the universal amplitude combination X(l, 0) /X(— 1, 0) the estimates of R^ 
are of less accuracy since the data are affected by the uncertainty of the exponent 7 
and the uncertainties of the metric factors a^, jl42j . These uncertainties increase 
the error-bars significantly. But nevertheless the two- and three-dimensional data agree 
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Figure 32: The universal scaling function X(l, i) for various dimensions. The dashed lines correspond 
to an e-expansion of a renormalization group approach |181| . The universal amplitude is obtained 
from the extrapolation ahh{ap 5p)^'^ — > 0. 



quite well with the RG-results of jl81j whereas unphysical results are obtained again 
for D = 1 from the e^-expansion. 

In this way the analysis of combinations of the universal scaling functions provides 
a detailed and instructive check of the accuracy of corresponding renormalization group 
results. Compared to critical equilibrium systems the e^-expansion provides astonish- 
ingly accurate estimates of certain amplitude combinations for e = 1 and e = 2. The 
e^-approximation fails only in low dimensions (e = 3), i.e., here, higher orders than 
©(e^) are necessary to describe the scaling behavior of directed percolation. 



3.2.5 Finite-size scaling behavior 

So far simulation data are considered where the correlation length ^ is small compared 
to the system size L. Thus, the data presented above do not suffer from finite-size 
effects, such as rounding and shifting of the anomalies. As pointed out in section 
12.2.11 an appropriate method to study finite-size effects of absorbing phase transitions 
within the steady state is to examine the field-dependence of the quantities of interest 
at the critical value 6p = for various system sizes (see FigureiriT]) . Similar to critical 
equilibrium scaling theory (see [2461 ^ review), it is assumed that the system size 
enters the scaling forms [Eqs. (|3. 37113. 3^]) ] as an additional scaling field 

p,{dp,h,L) ~ A^^ i?pb,,o(a>A,aft/i A'^,a^L A''^^), (3.94) 
a^Ap^{6p,h,L) ~ A^' Z)pbc,D(ap(5p A,a,,/iA'^,ai,L A^""^) . (3.95) 
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The choice of the corresponding scahng power i/^ is motivated by the phenomenolog- 
ical finite-size scahng theory which rests on the assumption that finite-size effects are 
controlled within the scaling regime by the ratio L/(j_ [2471 1248] . Approaching the tran- 
sition point, finite-size effects are expected to occur when O(^) = L. If both lengths, 
L and are significantly larger than all other length scales of a given system (e.g. the 
lattice spacing) the above finite-size scaling forms are universal. Thus, they may be 
used to identify a system's universality class, additionally to the so far considered scal- 
ing forms associated to the thermodynamic limit L — > oo. For example, it was shown 
by the use of universal finite-size scaling functions that the liquid-gas phase transition 
of a Lennard- Jones system belongs to the Ising universality class [2491 125UT 12511 . 

The index pbc, □ in Eq. ()3.94|) and Eq. ()3.95|) indicates that the universal finite-size 
scaling functions depend on the particular choice of the boundary conditions as well as 
on the system shape (see e.g. |253( 1254] I255| TT2T\ 1256] ) . But different lattice structures 
are contained in the metric factor a^. Throughout this work we focus on hypercubic 
lattices of size (aspect ratio 1) with periodic boundary conditions (pbc). Of course, 
the universal scaling functions Eqs. (|3. 37113. 31^ are recovered in the thermodynamic 
limit, e.g. 

Rpbc,n{x,y,oo) = R{x,y). (3.96) 

Additionally to the order parameter and its fluctuations we consider the fourth- 
order cumulant Q which is defined as (see e.g. |257( 1112] ) 

Other combinations of order parameter moments can be constructed and were inves- 
tigated for several systems exhibiting absorbing phase transitions |1961 12U51 12111 1245] . 
Here, we restrict our attention to Q since the cumulant behavior is sufficient for our in- 
vestigation |152j . For non- vanishing order-parameter the cumulant tends to Q = 2/3 in 
the thermodynamic limit. In case of a zero order parameter the cumulant vanishes if the 
order parameter is characterized by a Gaussian distribution symmetrically distributed 
around zero. The latter case is observed in equilibrium systems, e.g. the Ising model 
for T > T^. But for absorbing phase transitions the order parameter is non- negative 
per definition. Thus the order parameter is characterized by a non-trivial distribution 
and the above scenario does not apply. 

Nevertheless, it is expected that the cumulant obeys the scaling form 

Q{5p,h,L) ~ Qpbc,n(ap5p A,afe/i A'^,aiL A"*"^) . (3.98) 

Notice that no metric factor aq is introduced since the cumulant is already dimension- 
less. Choosing a^LX^'^^ = 1 we find for zero field 

Q(0,0,L) = Q{6p,0,L) 
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Sp—O 

= QpbcD (0,0,1) (3.99) 

which is universal despite of the boundary condition and the system shape. The univer- 
sal value Qpbc,n(0, 0, 1) corresponds to an intersection point if one plots Q as a function 
of the control parameter p for various system sizes L. Thus it is possible to determine 
the critical value from the common intersection point. This cumulant intersection 
method is very useful and was applied in numerous works (see for instance jll2j and 
references therein). But as pointed out in section r2. 2.11 it is a characteristic feature of 
absorbing phase transitions that the steady state order parameter moments (p^) vanish 
as soon as 0{(,±) = L even in the active phase. Thus the powerful cumulant intersection 
method can not be applied for absorbing phase transitions. 

In order to bypass this problem we consider the cumulant as a function of the 
conjugated field for 6p = 0. It is convenient to norm the universal scaling function 
Qpbc.D by the condition 

(0,1,1) = 0. (3.100) 

Since the metric factor is already known, the condition above can be used to de- 
termine the metric factor a^. Taking into account that the spatial correlation length 
scales at criticality as [Eq. ()3.41|) ] 

~ (a.h)-'^/'', (3.101) 

Eq. 1)3. 100(1 implies that the universal function Qpbc,D is positive if a^L > ^ and 
negative if a^L < a^^±. Note that in case of equilibrium phase transitions Eq. (|3.1()()|) 
is useless since the cumulant is usually positive. 

As mentioned above we focus on the scaling behavior at p = p^ and investigate the 
universal scaling functions i?pbc,n(0,x, 1), Dp^^ a{0, x, 1), <5pbc,D(0,x, 1), by plotting the 
rescaled quantities 

h, L) {a^ Lfl^^ , A/>,(0, h, L) {a^ L)'^'!-^ (3.102) 

and Q{0,h,L) as a function of the rescaled field x = a,j/i(a£,L)°'/'^-^ . Obviously, only 
a horizontal rescaling is needed in order to produce a data-collapse of the cumulant 
data. Thus the determination of the critical exponent i/^ via the fourth order cumu- 
lant is more accurate than a corresponding analysis of the order parameter or of the 
fluctuations. Therefore cumulants are widely used in order to estimate the correlation 
length exponent. The finite-size scaling analysis of site directed percolation, of the 
contact process, and of the pair contact process is presented in Figurel33l Using the 
best known values for the critical exponent (see Table^) perfect data-collapses are 
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Figure 33: The universal finite-size scaling analysis for the directed percolation universality class in 
various dimensions. The filled circles mark the condition Eq. 13.10011 . The horizontal lines correspond 
to the cumulant limit 2/3. The vertical lines separate the finite-size regime with a^L < a±^± (left) 
from the regime with UlL > a±^± (right). The data are obtained from simulations of systems sizes 
L = 64, 128, 256, 512 for D = 1, L = 64, 128, 256 for D = 2, and L = 32, 64, 128 for D = 3. 



observed. Note that different lattice structures do not affect the universal scaling func- 
tions. Whereas bcc lattice types are used for site directed percolation, simples cubic 
lattices are used for the contact process and pair contact process. As can be seen the 
order parameter and its fluctuations obey the expected algebraic field dependence for 
X ^ 1, i.e., as long as a^£,± ^ a^L. The fluctuation curves are characterized by a clear 
maximum signalling its singular behavior in the thermodynamic limit. Approaching 
the transition point the fourth-order cumulant tends to minus infinity for all models in 
all dimensions 

Qpbc,n(0, 3;, 1) — > — oo for x — > . (3.103) 

This behavior is caused by the vanishing steady state fluctuations, and it is assumed 
that the divergent fourth-order cumulant is a characteristic feature of all absorb- 
ing phase transitions |152j . A ratio that remains finite at criticality is obtained via 
U = {{pDipD - {p.){pl?)l{{p.){pt) - {P.){PI?) HSg. The value of [/ for x ^ char- 
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acterizes the universality class. Numerical investigations yield Ud=i = 0.833±0.011, 
Ua=2 = 0.704±0.013, and Ud=3 = 0.61±0.02. 

So far we considered the finite-size scaling behavior below the upper critical dimen- 
sion . The situation becomes more complicated above where the mean field theory 
applies. Naively, one would expect that the finite-size scaling behavior is described by 
the mean field exponents (Smf, ctmf, and u^^mf and e.g. the order parameter obeys 

p,{6p, h, L) - A-z^MF R^^^^(^a^5p A, a^h X^^^a^^L X-'^^mf) . (3.104) 

But as can be seen in Figure|^ for the order parameter cumulant no data-collapse 
occurs for z^x^mf = 1/2 reflecting the failure of phenomenological finite scaling within the 
mean fleld regime (see i78i for a readable review) . This breakdown of common finite-size 
scaling is well established in equilibrium although details are discussed controversially 
in the literature [2M1 EEOi ESI EMI EMI EMI EMI EMI EMI EHj- As pointed out by 
Brezin on the basis of exact calculations in the limit n ^ 00, the usual finite-size 
scaling assumption breaks for D> due to the existence of a dangerous irrelevant 
variable [7^]. In that case, the above scaling forms are replaced by j260j 

p^{5p,h,L) ~ A-'^MP i^^^^^(a^^pA,a,/^A"M^a^L A-"i.MF). (3.105) 

The modified scaling form refiects that finite-size scaling above is no longer con- 
trolled by the ratio L/£^qc L\5pY^'^^ but by the ratio L/Zq. In equilibrium, the so-called 
thermodynamic length Iq scales as |26()j 

/o cx |Jpr(2/3MF+7MF)/i?. (3.106) 

Below the upper critical dimension, the hyperscaling law 2/3 + 7 = 2 — a = Vj^D holds 
and Iq coincides with the correlation length ^ . It is instructive to express the divergence 
of the thermodynamic length in terms of the value of the upper critical dimension D^. 
According to the Ginzburg criterion (see chapter ISJ , its value is given in equilibrium 

by = (2/3mf + Tmf )/z^x,mf leading to 



6c |(5pr^^>MPi3c/i5 . (3.107) 



Thus, the scaling power of the finite-size scaling forms is given by 

<MF = z^±,mf^- (3.108) 

For example, the short range interaction Ising model is characterized by z^^^mf = 1/2 and 
D^ = A yielding i^J^mf — 2/D. This value is derived by a renormalization group analy- 
sis jm and is confirmed by numerical simulations j26()| 1262] . Furthermore, Eq. (|3.1()8|1 
also describes the modified finite-size scaling behavior of long-range interacting sys- 
tems [2591 126f| . In that case, the interactions decay as r"'^+'^' [a > 2 corresponds to 
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Figure 34: The universal finite-size scaling analysis for the directed percolation universality class above 
the upper critical dimension. Using the standard finite-size scaling form with /3mf = 1, ctmf = 2, and 
v±.M¥- = 1/2 no data-collapse occurs (left). Instead the finite-size scaling behavior is described by a 
modified scaling form using /3mf = 1, o"mf = 2, and i^J^^mf = 2/Z). The filled circles mark the condition 
Eq. I|3.1()()|l . The vertical lines correspond to the cumulant limit 2/3. The data are obtained from 
simulations of system sizes L = 8, 16, 32 for site directed percolation and L = 4, 8, 16 for the contact 
process. As can be seen in the right figure, the numerical data agree perfectly with the analytical result 
Q = 2/3 (1 - l/x) 1^1^ . 



short-range interacting systems), yielding the upper critical dimension = 2a as well 

as z^^,MF = 1/cT EMiiizniiiniiizs]- 

Analogous to the equilibrium scenario considered above, directed percolation ex- 
hibits a dangerous irrelevant variable within the mean field regime (the coupling con- 
stant u in Eq. p.8() ^. explaining the breakdown of usual finite-size scaling above D^. 
Following Eq. (|3.108() . the modified scaling form is characterized by the exponent 

<MF = ^- (3.109) 

since = 4 and j^x,mf = 1/2. As can be seen in FigurelHH a perfect universal data- 
collapse is obtained if the data of the five-dimensional models are scaled with z^j* = 
2/5. Preliminary results show that corresponding data of two-dimensional absorbing 
phase transitions with infinite particle hopping collapse onto the same universal scaling 
curves if v*^^ = l is used. This indicates that Eq. H3.1U9|) holds independently of the 
range of interactions. Furthermore, the exponent v* is D-dependent and therefore 
non-universal, whereas the finite-size scaling functions are universal. 

A recently performed field theoretical approach of the mean field scaling behavior 
of directed percolation confirms the analysis above. Similar to equilibrium |791 178j . a 
zero-momentum mode treatment allows to consider finite-size effects of periodic sys- 
tems |273j . Within the mean field regime it is possible to calculate analytically the 
critical exponent <7/v*^i^ = D [confirming Eg. (|3.1U9|) ] as well as the scaling func- 
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tions [2681 1258] 



(0,x,l) 



V2 




(3.110) 





) 



(3.111) 



As can be seen in Figure lH^ the numerical data agree perfectly with the analytical 
results. 

In conclusion, universal finite-scaling functions can be used additionally to the scal- 
ing behavior of the equation of state, the susceptibility, etc. to identify a system's 
universality class. Analogous to equilibrium, simple finite-size scaling is valid below 
but has to be modified within the mean field regime. 

3.3 Dynamical scaling behavior 

In this section we discuss the dynamical scaling behavior close to the transition point. 
First we consider how the order parameter decays starting from a fully occupied lattice, 
i.e., starting from a homogenous particle source. Second we investigate the activity 
spreading generated from a single active seed. Interpreting the initial particle density 
as another scaling field both scaling regimes are connected by a crossover scenario. 

3.3.1 Homogeneous particle source 

In section l2. 2. 21 we have seen that away from criticality the order parameter approaches 
its steady state value exponentially. At the critical point the order parameter decays 
algebraically 



where denotes the associated non-universal metric factor. The dynamical behavior 
of the order parameter can be incorporated in the scaling forms by an additional scaling 



field 

p^{6p,h,L,t) ~ X^'^ Rp^,^^a,{uii{apSp X,ahh X'' ,aLL X^''^,att X'''^\\) . (3.113) 



The scaling power of t has to equal the scaling power of the correlation time .^n, and 
the index full indicates the initial conditions. Setting Utt A~'^ll = 1 leads in the thermo- 
dynamic limit at criticality to 



(a* t) 



—a 



(3.112) 




(3.114) 
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Thus Eq. (|3.112|) is recovered for i?pbc,D,fuii(0, 0, oo, 1) = 1 and the exponents a, /?, and 
z^ll fulfill the scaling law 

a = ^. (3.115) 

A finite system size limits the power-law behavior above. The corresponding finite-size 
scaling form is obtained by choosing ai^LX~'^^ = 1 

pS5p = 0,h = 0,L,t) ~ (ar^L)-^/"^ i?pbc,D,fuii(0, 0, 1, a,t {a^L)-') , (3.116) 

where z = V\\lv^ denotes the dynamical exponent. Finite-size effects have to be taken 
into account for 

0(t) = tpss with tpss = (a^ ■ (3-117) 

For t <C tpss the scaling function obeys the power-law i?pbc.D,fuii(0, 0, 1, x) ~ 
whereas i?pbc,n.fuii(0, 0, 1, re) decays exponentially for x 1, i.e., i ^ tpss- This is 
shown in Figure l35l where the raw data as well as the rescaled data are presented for 
D = 1 and various system sizes. Again the data of all considered models collapse onto 
the single universal curve -Rpbc,n,fuii(0, 0, l,x). It is worth mentioning that the models 
are simulated with different update schemes. Whereas a parallel or synchronous up- 
date scheme is applied for site directed percolation, the contact process and the pair 
contact process are simulated with an asynchronous update, implemented by a random 
sequential update. The observed data-collapse confirms that different update schemes 
are contained in the non-universal metric factors Ot and do not affect the universal 
scaling functions. 

As already mentioned, the rapidity reversal of directed percolation implies the 
asymptotic equivalence of the survival probability and the order parameter decay at 
criticality [Eq. (ITT^ ] 

P^{t) ~ fi^pM- (3.118) 

The extension of the above dynamical scaling analysis to higher dimensions is redun- 
dant because a detailed analysis of the survival probability will be presented in the next 
section. But a comment about the rapidity reversal symmetry for site directed perco- 
lation is worth making. The order parameter p^{t) as well as the survival probability 
P^{t) for site directed percolation are plotted in FigurelHKl The survival probability 
is vertically shifted in order to check Eq. (|3.118j) . As can be seen both quantities are 
identical for all t. Thus the asymptotic equivalence (see the contact process data in 
FigurelTC)) holds in case of site directed percolation for all values of t, i.e., 

P^{t) = p^pM- (3.119) 

The factor p depends on certain non-universal parameters and can be determined easily. 
Starting from a fully occupied lattice, the density of active sites after one parallel update 



82 



Directed percolation 




Figure 35: The order parameter decay at the critical point for D — 1. The left figure shows the unsealed 
data of one-dimensional site directed percolation for L — 64, 128, 256 (from left to right). Additionally 
to the order parameter pg.(t) the survival probability P^{t) is presented in order to display the rapidity 
reversal symmetry {fi = y'2 — p^). The rescaled data, i.e., the universal scaling function 7?(0, 0, 1,2;) 
[see Eq. Ij^.H?^ ] is presented in the right figure for site directed percolation, the contact process and 
the pair contact process. The long-dashed line corresponds to the power-law behavior of the infinite 
system [Eq. Il3.ll2ll 1. 

step is given by p^{t = 1) = p, independently of the dimension and the lattice structure. 
On the other hand, the probability that the activity generated by a single seed survived 
one update step is P^{t = 1) = 1 — (1 — p)^fnn where Z{^^ denotes the number of forward 
next neighbors on the lattice. In case of 1 -|- 1-dimensional systems on a square lattice 
(•Zfnn = 2) we find at criticality p = -y/2 — . For the sake of completeness we just 
mention that Eq. (|.S.119j) holds with // = 1 for bond directed percolation 

3.3.2 Localized particle source 

As pointed out in section 12.2.21 measurements of activity spreading generated from a 
single seed have been widely applied in the last two decades. In particular they provide 
very accurate estimates of the critical value p^ as well as of the exponents 5, 9, and z. 
Here, we will focus on the scaling functions of the survival probability and the 
average number of active sites N^. At criticality both quantities obey the power-laws 
[Eq. ((2321)] 

apP, ~ {a,t)-\ a^N, - (a,^)^ (3.120) 

where 9 is often termed the critical initial slip exponent j274j . In this notation the 
rapidity reversal [Eq. (|3.12|) ] implies the scaling law a = 5. Sufficiently close to the 
critical point, and are expected to obey the scaling forms 

apP,{6p,h,L,t) ~ A-'^'^il Ppbc,D.scod(a/pA,a,/iA",a^L A-'^^,a,t A-'^il),(3.121) 
a^N,{5p,h,L,t) ~ A^^^ll iVpbc,D,sccd(ap5pA,a,/iA'^,aiLA-'^^,a,tA-''ll). (3.122) 
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Figure 36: The universal scaling functions Ppbc,n.aocd and iVpbc,n.accd of activity spreading for various 
dimensions. In case of the pair contact process the simulations are started from a natural configuration 
of inactive particles. System sizes L — 64, 128, 256, 512 are considered for D — 1, L = 64, 128, 256, 512 
for D = 2, and L = 16, 32, 64, 128 for D = 3. The dashed lines corresponds to the power-law behavior 
of the infinite system and a;*, respectively. 



Choosing attX~'^^^ = 1 the power-laws [Eq. (|3.120|) ] are recovered for Ppbc.D(0, 0, oo, 1) = 
1 as well as A^pbc,n (0, 0, cxd, 1) = 1. Again the finite-size scaling forms are obtained by 
setting a^L A^'^^ = 1 

apPMO,L,t) ~ (a^L)-''^ Ppbc,D,sced(0,0,l,a,t (a^L)-"), (3.123) 

a^NMO,L,t) ~ (a^L)^" iVpbc,D,sccd(0,0,l,a,t (a^L)-"). (3.124) 

Similar to the order parameter behavior the scaling functions Ppbc,n,scod(0, 0, 1, x) and 
■^pbc:,D,socd(0, 0, l,x) decay exponentially for t ^ tpss whereas they exhibit an algebraic 
behavior for t <^ tpss- 

Performing activity spreading simulations of site directed percolation and of the 
contact process the initial seed is implemented by a single particle on an empty lattice. 
For absorbing phase transitions with non-trivial absorbing states, like the pair contact 
process, the scaling behavior depends upon the nature of the initial configuration jl35j . 
In that case spreading activity simulations have to be performed at the so-called natural 
density of inactive particles jl351 11361 1152j . Starting with a random configuration, an 
absorbing state at criticality is prepared by the natural dynamics. Then an active seed 
is created by adding or moving a particle. Thus the numerical effort is significantly 
increased for systems exhibiting non-trivial absorbing states and only small system sizes 
are available by simulations. But nevertheless convincing data-collapses, including 
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the pair contact process data, are obtained and the corresponding universal scaling 
functions are presented in Figure|3Sl The values of the exponents used are listed in 
Tabled Thus activity spreading from a localized seed is characterized by the same 
universal scaling functions -Ppbc,n,sood and A'^pbc,D,sood for all considered models. 

The spreading exponents 6 and 9 are related to the critical exponents (5, and 
For example, the percolation probability Pp^vc is related to the ultimate survival 
probability Eq. 

Ppcrc{5p) = lim lim P^{5p,h = 0,L,t) 



~ A ^""ll Ppbc,n,scod(ap(5p A, 0,00,00) 

apOp\=l 



~ ttp^ (ap(5p)''''li Ppbc,n,scod(l,0, 00,00). (3.125) 

Comparing this result to Eq. 1)3. 3|) reveals the scaling law 

0' 

5 = — (3.126) 

^11 

In case of directed percolation, a = (5 is implied by the rapidity reversal [Eq. 1)3. 11(1 ] 
leading via Eq. H3.115|) to 

/3 = (3.127) 

as well as (5 = 

The critical initial slip exponent 9 is connected to the other exponents by the 
hyperscaling relation |136j 

3 0' D 

e + — + — = — . (3.128) 

This scaling law is fulfilled for directed percolation as well as other universality classes 
like compact directed percolation |275j . branching annihilation random walks with an 
even number of offsprings |154j , as well as absorbing phase transitions with a conserved 
field |152j . Using Eq. (|3.115|1 and Eq. (|3.126j) the hyperscaling law is often written as 

e + a + 6 = . (3.129) 
In case of directed percolation it reduces due to the rapidity reversal (q = 5) to 

[Eq. frm ] 

9 + 26 = 9 + — = —. (3.130) 

z 

Usually the hyperscaling law Eq. (|3.128() is derived via the pair connectness function 
which is defined as the probability that two sites are connected by a path of activity (see 
e.g. [2 Hl276|ll36|lll4| [7|). Here, we present an alternative and instructive way |277|ll52j 
which allows to connect the universal scaling functions Rf^n and A^^^ocd- Therefore we 
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assume that the initial particle density p^{t = 0) = p^.o appears as an additional scaling 
field 

p^{5p,h,L,t,p^^o) ~ (3.131) 

and the scaling function behaves asymptotically as (see and references therein) 

~ I X for X <^ 1 

i?p,,,o(0,0,oo,l,x) ~ <^ ^ ^ (3.132) 

[ const lor x ^ 1 . 

For example, within the mean field theory [see Eq. (|2.13|) ] the universal scaling function 
is given by 

-Rpbc,D(0,0, oo, l,x) = ^ _^ ^ . (3.133) 



Obviously the scaling function i?pbc,n is related to i?pbc,n 

full via 



-Rpbc,n,full — R'phc,a 



and due to = L^p^ related to A'^pbcn.scod by 



AT T ^ 

-^^pbc,n,sccd -'i'pbc,D 



Pa,0 = l 



Pa,0 = l/i^ 



(3.134) 



(3.135) 



The scaling power of the initial density is explained by the following crossover sce- 
nario: Starting at criticality from a low density of active sites (e.g. several seeds) the 
number of active sites increases as 6c Pa,o^^ until it reaches a maximum and crosses 
over to the expected asymptotic decay p^oit~°^. The crossover time is determined by 

O(ao/>a,o (flttco)^/""^) = 1, (3.136) 

corresponding to a merging of the survived (and former separated) clusters of activ- 
ity [7j. Despite of metric factors the crossover takes place at t^o 6c p^}J^^^^ and the 
scaling form Eq. (|3.13H) has to recover 6c for t <^ t^o whereas p^ 6c t~'^ has to be 
obtained for t ^ t^^. 

Focusing on criticality we find in the thermodynamic limit for Utt A~'^ll = 1 



,(0,0,oo,t,p.,o) ~ (a,t)-^/'^ii iipbe,a(0,0,oo,l,(t/teo)^/"~^). (3.137) 



Starting with a fully occupied lattice, i.e., with ao(atico)^^^ ^ = ^, the order parameter 
decays as 

p.(0,0,oo,i,l) ~ (a,i)-^/"ll i?pb,,n(0, 0,00,1, ao (a,i)^/""^) 
PS (att)~" ^pbc,n(0,0,cx), l,cx)) , 
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since -D/z — (5 > in all dimensions and where the scaling law a = (ijvw is used. 
Comparing this results to Eq. ()3.112|) leads to i?pbc,n (0, 0, oo, 1, oo) = 1. On the other 
hand, single seed initial configurations p^.o = are characterized by crossover times 

p,(0,0,oo,t ~ {a,t)-^l^\\ i?pbc,D(0,0,oo,l,(Vt,J^/^-^) 

~ (a,t)-^/^li+^/"-'^ aoL-^. (3.138) 

Thus a(,N^ = p^L^ ~ (fltO^ is recovered for t <^ t^^, confirming the hyperscaling 
relation Eq. (|3.128|) . Furthermore the metric factor of the initial configuration is given 
by ao = l/ojv. 

In this way the asymptotic order parameter decay [Eq. ()3.112p ] and the initial par- 
ticle growth [Eq. (|3.12()|) ] emerges as different scaling regimes of the generalized scaling 
function [Eq. (|3.131|) ]. The full crossover can be observed in simulations if the particular 
value of the initial density of active sites leads to 1 ^ t^^ <^ tpss (see jl52| for a detailed 
discussion). For example the scaling regime 6c t~°^ can not be observed in simulations 
starting from single seed since tpss < ^co, i-e., finite-size effects take place before the 
algebraic particle decay starts. On the other hand, too large initial densities p^.o ~ 1 
cause too small crossover times O{tco) = 1 and the short time scaling regime {N^oat^ 
for 1 <C t <C tco) does not occur. But intermediate initial densities allow to investigate 
the crossover behavior. An alternative way to observe the full crossover is to perform 
simulations with various initial conditions. According to Eq. (|3.136|) and Eq. (|3.137|) all 
curves corresponding to different initial configurations collapse onto the universal func- 
tion i?pbc.D(0, 0, oo, l,x) by plotting p^ (att)^/''li as a function of x = a^^^p^^ [att)^/'^"^ . 
This is shown in Figure l37l for 1 -|- 1-dimensional site directed percolation and the con- 
tact process. Note that the assumed asymptotic behavior [Eq. (|3.132j) ] of the universal 
scaling function i?pbc,D(0, 0, cxd, 1, x) is perfectly recovered confirming the scaling form 
Eq. 1)3. 131(1 . Furthermore, the data-collapse reveals that the crossover between two scal- 
ing regimes of a given universality class can be described in terms of universal scaling 
functions. We revisit crossover phenomena in a more generalized context in section 
There we address the question if the crossover between two different universality classes 
can be described in terms of universal functions. 

The hyperscaling law [Eq. (|3. 128(1 ] still implies more. Like all hyperscaling laws it 
is fulfilled for the mean field exponents at the upper critical dimension 

^MP + ^ + <5mp = — . (3.139) 

Within the mean field theory, the activity spreading of most absorbing phase transitions 
is described by a branching process. The mean field values ^mf = and (5mf = 1 lead 
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Figure 37: The universal scaling function 7?pbc,n (0, 0, oo, \,x) describes the crossover from initial par- 
ticle growth N^Sit^ to the asymptotic order parameter decay paOcf~". The data are obtained from 
simulations of site directed percolation and the contact process. Different initial particle densities 
Pa,o = n/L with n — 1, 2, 4, . . . , 512 and system sizes L = 1024,2048 are considered. Data that are 
affected by finite-size effects {0{t) > tpss) are skipped for a better visualization. The long-dashed lines 
correspond to the asymptotic scaling behavior of the universal scaling function Eq. H3.182^ . Neglecting 
the metric factors each model would be characterized by its own scaling function. In order to scrutinize 
the data which covers several decades it is usual to examine the crossover via the effective exponent (left) 
dinx In 7?(0, 0, oo, 1, a::). For the sake of completeness, the effective exponent of the 1 -1- 1-dimensional 
models is compared to the mean field behavior [dot-dashed line, see Eq. 13.13311 1 . Extended simulations 
reveal that the crossover scaling function depends on the dimension and approaches the mean field 
curve with increasing dimension (not shown). 



to 

= z^Jl + ^] . (3.140) 

\ ^11, MF / 

We have already seen that the exponents /3mf and i^y^MP can be derived within a simple 
mean field approach basing on the transition rates. Therefore, the above scaling law 
allows to determine the upper critical dimension. Furthermore a lower bound on is 
obtained by taking into account that the order parameter exponent is positive definite, 
yielding 

> Zmf . (3.141) 
In case of diffusion- like mean field processes the latter result reduces to D^>2. 



3.4 Conclusions and comments 

The scaling behavior of directed percolation is recognized as the paradigm of non- 
equilibrium phase transitions from an active to an absorbing state. The widespread 
occurrence of such models is reflected by the universality hypothesis of Janssen and 
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Grassberger. Despite its ubiquity in theoretical models no experiments observing di- 
rected percolation critical behavior are known so far. The lacking experimental real- 
ization is explained by the influence of quenched disorder or the impossibility to realize 
a perfect non-fluctuating absorbing state |2j. For a summary of experimental efforts, 
including catalytic reactions, flowing granular matter |2781 1279] . the onset of turbu- 
lence jlU2| , we refer the interested reader to [3 l28Uj . But due to its simplicity and 
robustness it is still believed that experimental counterparts to directed percolation ex- 
ists. Thus an experimental realization of directed percolation remains an outstanding 
problem of top priority |281j . 

Additionally to the experimental situation several theoretical aspects are also still 
open. For example directed percolation is not solved analytically in 1 + 1-dimensions. 
Thus, only approximative estimates (see e.g. |l282j ) exist of the still unknown critical 
exponents. Furthermore, the universality hypothesis still awaits a rigorous proof. It is 
worth mentioning that the universality class of directed percolation is even larger than 
expected from the hypothesis. In particular, the pair contact process and the threshold 
transfer process belong to the directed percolation universality class although they 
are not characterized by a unique absorbing state. Thus the hypothesis defines only 
sufficient conditions but fails to describe the DP universality class in full generality 
(see (Tj for detailed discussion). A refinement of the universality hypothesis is part of 
the larger task to provide a classification scheme for non-equilibrium phase transition. 
Compared to equilibrium a full classification of non-equilibrium phase transitions is 
still open. 
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Directed percolation is considered as a paradigm for non-equilibrium phase transitions 
into absorbing states. Different universality classes are expected to occur in the presence 
of additional symmetries. For example, a particle-hole symmetry is associated to the 
universality class of compact directed percolation (see 13.1.2(1 . Branching annihilating 
random walks with an even number of offsprings obey parity conservation (PC) and 
define another universality class |283[ I284j . In the following we consider a class of 
absorbing phase transitions which is characterized by particle conservation. According 
to Rossi et al. the additional conservation law leads to the universality class of stochastic 
absorbing phase transitions with a conserved field |285j . Similar to the universality 
hypothesis of Janssen and Grassberger the authors conjecture that all stochastic models 
with an infinite number of absorbing states in which the order parameter is coupled 
to a non-diffusive conserved field define a unique universality class. This universality 
class is of particular interest since it is related to the concept of self-organized criticality 
(SOC) jlL)51llU6lllU7lllL)^ . Compared to directed percolation this universality class is 
not well established so far. Besides field theories |2861 12871 1199j . series expansions j288j 
and path integral representations |289j most quantitative results are obtained from 
simulations. In particular, a systematic e-expansion is still lacking. In the following 
we review the numerical results and focus again on the universal scaling behavior. 
Several lattice models are introduced and the corresponding steady state and dynamical 
scaling behavior is investigated. The determination of the value of the of upper critical 
dimensions as well as the relation to self-organized criticality are discussed. 



4.1 Lattice models of absorbing phase transitions with a 
conserved field 

In the following we consider various lattice models obeying particle conservation. The 
first one is a modification of Manna's stochastic sandpile automaton |29nj . This model 
was intensively investigated in the context of self-organized criticality (see e.g. |2911 
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l292l [2931 Ens EHEI EHU EMI Enni)- The second model is a modification of the 
threshold transfer process that obeys particle conservation. This model is well suited 
for analytical treatments at the mean field level. Furthermore it allows to investigate 
the crossover between the universality class of directed percolation and the universality 
class of stochastic absorbing phase transitions with a conserved field. The third model is 
the so-called conserved lattice gas. Compared to the bosonic Manna model, no multiple 
particle occupation is allowed in the conserved lattice gas. But more than the other 
models, the Manna model can be considered as the paradigm of the universality class of 
stochastic absorbing phase transitions with a conserved field. Since universality classes 
are often labeled by the simplest model belonging to them it is termed the Manna 
universality class in the following. 

Additionally to these lattice models a conserved reaction-diffusion model is dis- 
cussed in the literature |30ni \277\ I287| I.SOlj . The model is defined by the reaction 
scheme 

A — > B, A + B — > 2A, (4.1) 

and where ^-particles diffuse with diffusion rate and -B-particles with rate Dg, 
respectively. The reaction scheme describes systems where a conserved field is coupled 
to a non-conserved order parameter (^-particles). Obviously, the total number of 
particles = N^ + Ng is conserved and a continuous phase transition occurs by varying 
A^. In the limit Dg the reaction-diffusion system belongs to the universality class 
of interest, i.e., of absorbing phase transitions where the order parameter is coupled to 
a non-diffusive field |285| I.SOlj . Field theories are derived from the reaction-diffusion 
model but the corresponding renormalization group approaches run into difficulties and 
the results are controversially discussed in the literature |S1 12871 11991 1176j . 

4.1.1 Manna model 

The Manna model was introduced in j290j as a stochastic sandpile model in which 
integer values represent local energies, number of sand-grains or particles. The Manna 
model is a bosonic lattice model, i.e., it allows unlimited particle occupation of lattice 
sites (n = 0,l,2,...). Lattice sites are considered as inactive if the particle occupation 
is below a certain threshold n < N^. For n> the lattice site is active and the particles 
are moved to nearest neighbors, selected at random for each particle, i.e., 

n — > for all sites with n> . (4.2) 

The dynamics is sketched in Figure l38l Additionally to this original Manna model |290j 
a modified version is considered in the literature where the occupation number of active 
lattice sites is not reduced to zero but 

n — > n — for all sites with n> , (4.3) 



4.1 Lattice models of absorbing phase transitions with a conserved field 



91 



O , O 



0,0,0 



o , o 



o 



Figure 38: Sketch of the dynamics of the 1 + 1-dimensional Manna model with = 2 according to 
the original rules {n ^ for n> N^). Filled circles mark active particles whereas non-active particles 
are marked by open circles. The arrows denote possible particle movements in the next update step 
according to a particular representation of the stochastic rules. 



and particles are randomly transferred to the nearest neighbors. In the latter case 
the Manna model is analytically tractable if the active particles are distributed in a 
particular way j292j . In both versions the particle transfer conserves the particle num- 
ber, i.e., no particle creation or annihilation processes take place. Applying periodic 
boundary conditions the total number of particles N is conserved. The particle den- 
sity p = N/L^ is the control parameter of a continuous phase transition from an 
inactive phase to an active phase. The active phase is characterized by a non-zero 
density of lattice sites occupied by at least particles. The density p^ corresponds 
to the order parameter |3U2j . A configuration is absorbing if it contains only empty 
and single occupied lattice sites. Thus the Manna model is characterized by an infinite 
number of absorbing configurations in the thermodynamic limit. Similar to the univer- 
sality class of directed percolation it is possible to apply an external field h which is 
conjugated to the order parameter. But considering absorbing phase transitions with a 
conserved field one has to take care that the external field does not change the particle 
number. A possible realization of the external field was developed in [l.SSj where a 




Figure 39: The order parameter behavior of the 2 + 1-dimensional Manna model and of the conserved 
threshold transfer process (CTTP) for various field values (from h = 5 10"^ to h = 10""). For non-zero 
field (lines) pa exhibits an analytical behavior. For zero field (symbols) the order parameter vanishes 
algebraically by approaching the transition point. The data are obtained from simulations of various 
system sizes L < 2048 with periodic boundary conditions. 
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diffusion-hke field is implemented. The field triggers movements of inactive particles 
which may be activated in this way. The corresponding order parameter behavior of 
the two-dimensional Manna model is shown in Fieure lHUl 

The driven-dissipative version of the Manna model was investigated intensively in 
the context of self-organized criticality (see e.g. 1^ 1^ 1^ 'l . 

Following Ben-Hur and Biham |291j the universality class of self-organized critical sys- 
tems is determined by the way how active particles are distributed to the next neighbors 
(deterministic, stochastic, directed, undirected, etc.). At first glance, the Manna model 
seems to be characterized by a stochastic particle transfer. But this is correct only for 
moderate, i.e., not to large values of the threshold N^. In the limit ^ oo the active 
particles are equally distributed to the nearest neighbors, i.e., the stochastic character 
of the particle transfer is lost. In that case the Manna model exhibits a different scaling 
behavior. The crossover between the two universality classes was observed in numerical 
simulations and confirms the conjecture of |291j . Here, we focus our attention on the 
scaling behavior of the original Manna model [Eq. 1)4. 2(1 ] for A^^. = 2 in various dimen- 
sions. The modified version [Eq. ()4.3|) ] is considered exemplarily for D = 2 in order to 
show the identical scaling behavior. 

4.1.2 Conserved threshold transfer process 

The conserved threshold transfer process (CTTP) |285j is a modification of the thresh- 
old transfer process (see section I3.1.4|) . Again lattice sites may be empty, occupied by 
one particle, or occupied by two particles. Empty and single occupied lattice sites are 
considered as inactive, whereas double occupied lattice sites are considered as active. 
In the latter case one tries to transfer both particles of a given active site to randomly 
chosen empty or single occupied sites (see Figure HO]) . Thus in contrast to the threshold 
transfer process no creation and annihilation processes of isolated particles take place 
and the particle number is conserved. Similar to the Manna model the conjugated 
field is implemented by a diffusion- like field jl.'fflj . The order parameter behavior is 
shown in FigureOHfor various external fields. Numerical simulations turn out that the 
conserved threshold transfer process exhibits the same scaling behavior as the Manna 
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Figure 40: Sketch of the dynamics of the 1+ 1-dimensional conserved transfer threshold process (CTTP). 
Filled circles mark active particles whereas non-active particles are marked by open circles. The arrows 
denote possible particle movements in the next update step according to a particular representation of 
the stochastic rules. Note that the relaxation event presented on the right site is deterministic. 
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model for D >2 |3n3j . The situation is more complicated in one-dimensional systems 
where a splitting of the universality class occurs jl52j . The dynamical reduction causes 
a change of the stochastic character of the dynamics of the CTTP. Whereas all particle 
movements are stochastic in the Manna model, roughly 40% of the relaxation events of 
the one-dimensional CTTP are deterministic (see Figure EO]) . Furthermore, a pertur- 
bation triggered by the external field inside a cluster of inactive sites performs a simple 
random walk. This pathologic behavior differs from that of the Manna model jl52j . 
Therefore we focus in the following on the conserved threshold transfer process in D > 2. 

The CTTP can be considered as a modification of the Manna model exhibiting a 
restriction of particle occupation. This occupation restriction simplifies an analytical 
treatment significantly. For example, the CTTP can be analyzed easily within a mean 
field approximation whereas the multiple particle occupation of the Manna model leads 
to a system of coupled equations |3n4j . Therefore, we focus on the CTTP in order to 
derive the mean field behavior of the universality class j3U5j . Considering uncorrelated 
particle behavior, the associated probabilities for a change of the number of active sites 
by An^ are 

p(An, = -1) = p^{pI + 2p,p^) , 

p{An,= 0) = p^{2p,pi + 2pip^ + pl) + pih{p, + p^), 

p(An, = 1) = p^p^ + p^h. (4.4) 

Here, p^ denotes the order parameter, p-, corresponds to the density of single occupied 
sites and p^ is the density of empty lattice site. This reaction scheme leads to the 
differential equation 

d,p^{p,h) = + 2p- 4p, + pI) + h{p - 2pj2^ (4.5) 

where the normalization Pc + Pi + Pa. = 1 as well as the particle conservation p = pi + 2p^ 
are used |3n5j . For zero field, the steady state condition (?t = leads to 

=0 V -I + 2p -4p, + pI = 0. (4.6) 

The first equation corresponds to the absorbing state and the second equation yields 
the order parameter as a function of the particle density |3U51 l3U6j 

= 2 - V5^^ for p>l/2. (4.7) 

The order parameter vanishes in leading order as p^ = 6p^ /4 with the reduced control 
parameter 5p = {p—pc)/p^ and (3=1. For non-zero conjugated field the order parameter 
scales at the critical density p^ = 1/2 as p^ = /i^/'^/4 with a = 2. Close to the critical 
point, the asymptotic behavior of the equation of state is given by [HUH] 
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We focus now on the dynamical behavior. The zero field behavior is determined by 
the differential equation 

dtp. = pA¥ - ^P. + pD- (4.9) 

For sufficiently small order parameter higher order terms 0{p^) can be neglected, yield- 
ing for 6p ^ 

pM = , . — (^-^o) 

4 - (4 - ()/o//o,,o) exp [-bp t) 
Asymptotically ^ oo) the order parameter behaves as 

PaW ~ -bpU-^ ' e*S (4.11) 



*<0 \ Pa.O 

*>o 4 16 V /Oa,0/ 

Apart from criticality, the steady state solutions (p^ = (5/0/4 and p. = 0) are approached 
exponentially, independent of the initial condition p. ^. The corresponding correlation 
time is given by = with = 1. At the critical density {5p = 0) the order 

parameter exhibits an algebraic decay, i.e., 

1 1 

PZo + 4t 4t 



pM = I , t — • (4-13) 



Thus the mean field behavior of the conserved transfer threshold process agrees 
with the mean field solution of directed percolation, i.e., the same critical exponents 
(/? = 1, (T = 2, z^ii = 1), and the same universal scaling functions are obtained. The 
corresponding non-universal metric factors are given by Op = 1/4, = 1/16, as well as 
at = 4. Although the scaling behavior of the CTTP and of directed percolation differ 
in low dimensions |139j they coincide on the mean field level. We will address this point 
in detail below. 



4.1.3 Conserved lattice gas 

The conserved lattice gas (CLG) |285j is a stochastic variant of a model introduced by 
Jensen |307) . In both cases, lattice sites may be empty or occupied by a particle. The 
restriction to single particle occupation differs from the bosonic Manna model. Moti- 
vated by experiments on flux flow in type-II superconductors a repulsive interaction is 
assumed |3r)7j . i.e., a particle is considered as active if at least one of its neighboring 
sites on the lattice is occupied by another particle (see FigureHTj) . If all neighboring 
sites are empty, the particle remains inactive. Active particles are moved in the next 
update step to one of their empty nearest neighbor sites, selected at random. Obvi- 
ously, the dynamics obey particle conservation and all configurations consisting only of 
isolated particles (without occupied nearest neighbors) are absorbing. Increasing the 
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Figure 41: Sketch of the dynamics of the 1 + 1-dimensional conserved lattice gas (CLG). Neighboring 
particles are considered as active (filled circles) whereas isolated particles remains inactive (open circles). 
The arrows denote how the active particles are moved in the next update step. 

particle density, a continuous phase transition from a non-active to an active phase 
takes place |285(I137] . Considering simple cubic lattices an upper bound on the critical 
density is = 1/2. It is known that the one-dimensional system exhibits a trivial 
behavior with = 1/2 (see also Ref. 22 in |285| ) whereas non-trivial behavior with 
< 1/2 occurs for D>2. It is worth mentioning that the one-dimensional CLG is 
characterized by a deterministic dynamics, i.e., the one-dimensional CLG does not be- 
long to the universality class of stochastic absorbing phase transitions with a conserved 
field. Therefore, we will consider the CLG only in higher dimensions in this work. As 
for the Manna model and the CTTP, the conjugated field is implemented by particle 
diffusion E 



4.2 Steady state scaling behavior 

In the following section we examine the steady state scaling behavior of three different 
models in various dimensions. We focus our attention on the scaling behavior of the 
order parameter, its fluctuations, as well as on the scaling behavior of the susceptibility. 
Again, the Hankey-Stanley scaling forms Eqs. (|H.;-i7IH.H9j) normed by Eq. (|H.44j) are used 
in order to show that the three models of interest belong to the same universality class. 
In all cases the conjugated field h is implemented by particle diffusion and 5p denotes 
again the distance to the critical point, i.e., 5p = {p — p^)/ Pc- A comparison of the 
scaling behavior of high dimensional systems with the mean field functions derived 
above leads to a reliable determination of the upper critical dimension. Additionally, 
we consider briefly universal amplitude combinations as well as steady state finite-size 
scaling. 

4.2.1 Equation of state and fluctuations 

In |285j . Rossi et al. conjecture that all stochastic models with an infinite number of 
absorbing states and an order parameter coupling to a non-diffusive conserved field 
define a unique universality class that is different from directed percolation. At the 
beginning we verify that the Manna model, the CTTP and the CLG belong to the 
same universality class |285| . The scaling analysis of the equation of state and of the 
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fluctuations are plotted in Figure H2l The corresponding exponents are listed in TablelHl 
The data of all models collapse onto the same scaling functions R{x, 1) and D{x, 1) 
despite the fact that bosonic and non-bosonic models on different lattice types using 
different update schemes are considered. This result is a strong manifestation of the 
robustness of the Manna universality class with respect to variations of the microscopic 
interactions. Furthermore, the data-collapse confirms the accuracy of the numerically 
estimated values of the critical exponents (see TablelHl). 

After verifying that the Manna model, the CLG as well as the CTTP share the 
same scaling behavior we show that the universality class differs from that of directed 
percolation. Therefore we compare in Figure BHl the order parameter exponents P as well 
as the scaling functions R{x, 1) for the DP class and the Manna class. As can be seen 
both classes are characterized by significantly different scaling functions and critical 
exponents. Thus as conjectured in |2S5j, stochastic absorbing phase transitions with 
a conserved field define a unique universality class different from directed percolation. 
Evidently, the difference between the scaling functions and the exponents reduces with 
increasing dimension. Since both classes are characterized by the same mean field 
equation R{x, y), the question of the value of the upper critical dimension of the Manna 
universality class arises immediately. Unfortunately, no convincing answer is provided 
by field theoretical approaches so far. The contrary values D^ = A |287j and D^ = Q jl99) 
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Figure 42: The universal scaling functions R{x, 1) and D{x, 1) of the Manna universality class for 
D = 2. The data are plotted according to Eq. (EISJ. The CTTP data are obtained from simulations 
on a square (sq) and honeycomb (he) lattice. In case of the Manna model, data of the original model 
[Eq. I4.2ll 1 using random sequential (rs) and parallel (p) update as well as data of the modified model 
(n — > 0) [Eq. H4..'i^ ] are presented. All considered models are characterized by the same universal scaling 
functions demonstrating the robustness of the Manna universality class with respect to variations of 
the microscopic interactions. For all considered models the scaling plots contain at least four diff'erent 
curves corresponding to four different field values (see e.g. Figure l^^ . The circles mark the conditions 
^(0, 1) = 1 and D{0, 1) = 1, respectively. 
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Figure 43: The universal scaling functions R{x, 1) of the Manna and directed percolation (DP) uni- 
versality class for D = 1,2,3 (left). Data of the CLG {D = 2,3), the CTTP {D = 2,3) as well as of 
the Manna model {D — 1, 2, 3 with parallel and random sequential update) are shown. The two- and 
three-dimensional data are vertically shifted by a factor 8 and 64 in order to avoid overlaps. The circles 
mark the condition -R(0, 1) = 1. At least four different field values are plotted for each model and di- 
mension. As can be seen the universal scaling functions of the Manna and DP class differs significantly, 
reflecting the different universality. The order parameter exponents /3 of both classes are presented in 
the right figure. In all cases the error-bars are smaller than the symbol sizes and no overlap within the 
error-bars occurs (see Tableland Table^J. Additionally to the exponent values the renormalization 
group estimates of the DP class [Eq. H3.14^ ] are shown. 



are conjectured from field theories. In order to obtain numerically a reliable estimate 
of the same method as described in section 15.2.31 is applied, i.e., the data of high- 
dimensional models are compared to the analytically known mean field equations. The 
result is presented in Figure lHl where the order parameter as well as fluctuation data of 
five- and six-dimensional models are shown. Furthermore, the data of a modified version 
of the two-dimensional CLG with random neighbor hopping of active particles are 
presented There, the unrestricted particle hopping breaks long-range correlations 

and the system exhibits mean field scaling behavior. In all cases the numerical data 
are in a perfect agreement with the mean field scaling functions R{x,l) and D{x,l), 
respectively. Thus the data-collapse leads to the upper bound on the critical dimension 
< 5, clearly excluding = 6 as predicted in |199j . 

The observed mean field scaling behavior for D>5 and the observed non-mean 
field scaling behavior for D <3 restricts the value of the upper critical dimension to 
3 < < 5. Neglecting non- integer values (which may occur for example in systems 
with long-range interactions) we conclude = 4. This value is also supported by the 
Ginzburg criterion (see chapter 
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Figure 44: The universal scaling function of the order parameter R{x, 1) (left) and the fluctuations 
D(x, 1) (right) for D = 5, 6. Additionally to the five- and six-dimensional models we present data of a 
two-dimensional CLG model with random neighbors (rn) At least four different field values are 

plotted for each model. All presented data agree perfectly with the corresponding mean field functions 
Eq. J^THI and Eq. igilTJ, respectively. 



Furthermore, ZJ^. = 4 is confirmed by direct numerical investigations |141j . Again, the 
scaling behavior at D = 4 is characterized by logarithmic corrections (see Figure ll^ . 
The performed numerical analysis reveals the estimates for the logarithmic correction 
exponents B = 0.20±0.06, S = 0.35±0.05, and ^' = B - S = -0.15±0.08 OT- As 
expected these values differ significantly from the corresponding values of directed per- 
colation B = S = 1/3 |223j . Since the logarithmic corrections make the analysis of 
the numerical data notoriously difficult it is desirable to calculate the exponents B 
and S within a renormalization group approach. Unfortunately, these calculations are 
unavailable so far. 

In summary, the investigations of the equation of state reveals that all considered 
models belong to the same universality class which differs from the directed percolation 
behavior for D < D^. The upper critical dimension is four, which is also supported by 
simulations of the Manna model in the context of self-organized criticality |3n8| 1309) . 
It is worth mentioning that the directed percolation and the Manna universality class 
share the same mean field scaling behavior above = 4. A similar scenario is known 
from equilibrium phase transitions, where the Ising model, the XY-model as well as the 
Heisenberg model exhibit different scaling behavior below but share the same mean 
field behavior above D^. Another example is known from restricted random walks, like 
self-avoiding random walks |31()j or loop-erased random walks [31 H 13 121 1313j , which 
differ below D^ = A whereas they exhibit the scaling behavior of a simple random walk 
above flTTlj . 



4.2 Steady state scaling behavior 



99 




Figure 45: The universal scaling functions of the order parameter (upper left) and its fluctuations 
(upper right) at the upper critical dimension = 4. The logarithmic correction exponents are given 
by B = 0.20, E = 0.35, and '1/ = —0.15 305,, . The lower flgures show the order parameter for zero 
field (left) and at the critical control parameter (right), respectively. The order parameter is rescaled 
according to Eq. 13.761 and Eq. 13.7711 . Approaching the transition point (/i — > and &p —> 0) the data 
tend to the function f{x) — x (dashed lines) as required (see text). The circles mark the condition 
,R(0, 1) = 1 and D{0, 1) = 1, respectively. 

4.2.2 Universal amplitude combinations 

Universal amplitude combinations are related to particular values of the universal scal- 
ing functions. Therefore we expect that the amplitude combinations of the Manna class 
differ from the corresponding values of the DP class, as the scaling functions R and D 
do. Here, we focus on the susceptibility ratio X{+1, x) /X{—1, x) [see Eq. (|:i85|) ]. The 
susceptibility is obtained by performing the numerical derivative of the order parameter 
with respect to the conjugated field. The data of the Manna model, the conserved lattice 
gas, and of the conserved threshold transfer process are shown in Figure ll^ for various 
dimensions. For each dimension the data of the different models collapse onto a single 
universal curve. Similar to directed percolation, the systems are well described by the 
mean field solution far away from the transition point because the critical fluctuations 
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Figure 46: The universal scaling function X(l, a;)/X(— 1, x) for various dimensions and models. The 
dashed lines mark the corresponding behavior X(l, i)/X(— 1, x) for the directed percolation universality 
class. The universal amplitude X(l, 0)/X(— 1, 0) is obtained from the extrapolation ahh\ap Sp\~'^ — > 
leading to 0.41 ± 0.05 {D = 2) and 0.80 ± 0.05 (D = 3). 



are suppressed. Approaching criticahty {h — > 0) the critical fluctuations increase and a 
crossover to the L>-dimensional behavior takes place. Below the upper critical dimen- 
sion the data of the susceptibility ratio differ from the DP behavior. In particular, the 
zero field universal amplitude ratio X(-|-l, 0)/X(— 1, 0) reflects the different universality 
classes. Numerically this ratio is obtained by an extrapolation of the susceptibility 
ratio X(-|-l, j;)/X(— 1, x) to x — > 0. The estimated values for D = 2,3 are listed in 
TableEl In order to estimate the value of the one-dimensional Manna universality class 
additional simulations are required since no saturation occurs so far. 



4.2.3 Finite-size scaling behavior 

So far we considered systems of the Manna universality class where the correlation 
length ^ is small compared to the system size L. Additionally to these simulations 
of infinite systems, finite-size effects are analysed in the literature |152) . Analogous to 
the method described in section [3.2.51 the field dependence of the order parameter p^, 
of its fluctuations Ap^ as well as of the fourth order cumulant Q at the critical point 
are examined. Beside of the provided universal scaling functions, the finite-size scaling 
analysis of the cumulant provides accurate estimates of the correlation length expo- 
nent The results are listed in Table|31 The finite-size scaling analysis is shown in 
Figure ETI for various models and D < D^. As can be seen, the order parameter and its 
fluctuations obey the expected algebraic fleld dependence for <^ a^L. The fluctu- 
ation curves are characterized by a clear maximum, signalling its singular behavior in 
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Figure 47: The universal finite-size scaling analysis of the Manna universality class in various dimen- 
sions. In case of the one-dimensional models, the expected splitting of the universality class is observed 
clearly. For D > 2 all models considered belong to the same universality class. The filled circles mark 
the condition Eq. l3.10Ut . The vertical lines correspond to the cumulant limit 2/3. The data are ob- 
tained from simulations of systems sizes L = 2048, 4096, 8192 for D = 1, L = 64, 128, 256 for D = 2, 
and L = 16, 32, 64 for _D = 3. 



the thermodynamic limit. Approaching the transition point, the fourth-order cumulant 
tends to infinity for all models in all dimensions. Good data-collapses are observed in 
case of the two- and three-dimensional models. In case of the one-dimensional models 
the expected splitting of the universality class, i.e., the anomalous scaling behavior of 
the 1 -|- 1-dimensional CTTP is seen clearly. Thus additionally to the universal scaling 
functions of infinite systems, universal finite-scaling functions are very useful in order 
to identify a system's universality class. 



4.3 Dynamical scaling behavior and hyperscaling relations 

In the following we consider the dynamical scaling behavior of the Manna model, the 
conserved lattice gas, and the conserved threshold transfer process. Similar to directed 
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percolation we investigate both homogeneous and locahzed particle sources. The ob- 
tained estimates of the critical exponents allow to check certain hyperscahng relations. 
It turns out that hyperscahng holds, although the Manna universality class is charac- 
terized by a scaling anomaly which is not understood so far. 

4.3.1 Homogenous and localized particle source 

Starting from a homogeneous particle source the order parameter is expected to obey 
an algebraic decay at criticality described by the exponent a, analogous to directed 
percolation [Eg. (|H.112|) ]. Performing simulations, = p^L^ particles are randomly 
distributed on the lattice. In this way an initial course grained homogeneous particle 
density is created which differs significantly from the corresponding initial density of 
directed percolation with p^ Q = 1. For example, a random distribution of particles leads 
for the two-dimensional CTTP to an initial density Pa,o ~ 0.1703. Of course, different 
values are obtained for the Manna model and the conserved lattice gas. Nevertheless, 
asymptotic universal scaling behavior is observed for the order parameter decay. The 
corresponding scaling plots [see Eq. 116(1 ] are shown in Figure llSl The obtained 
estimates of the exponents a, z are listed in TableOl The values of the metric factors 
can be found in |152j . 

The spreading exponents 6 and 9 are obtained by simulations starting from a local- 
ized particle source, i.e., from a single seed. Due to the non-trivial absorbing state for 
the models considered the simulations have to be performed at the natural density |135j . 
Therefore, an absorbing state at criticality is prepared by the dynamics |152j . The 
scaling analysis of the survival probability as well as of the average number of active 
particles are presented in Figure EHl The obtained values of the dynamical exponent z 
agree with those obtained from the order parameter decay measurements. The results 
are listed in TableOl It is worth mentioning that the Manna universality class is char- 
acterized by the inequality a ^ 5. For example, the two-dimensional systems yield the 
estimates a = 0.419 ± 0.015 and 5 = 0.510 ± 0.020 HHS]. This indicates the expected 
violation of the rapidity reversal [that implies Eq. (|3.12|) ]. This point is revisited in 
detail in the following section. 

4.3.2 Hyperscaling relations 

In this section we check several hyperscaling laws. The validity of the scaling law is 
immediately related to the number of independent exponents. Compared to directed 
percolation the number of independent exponents is expected to be larger than three 
due to the violation of the rapidity reversal. At the beginning we consider the scaling 
relation 7' = Dv^ — 2/3 [Eq. (|3.4U() ] . The corresponding data are plotted in Figure UHl 
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Figure 48: The universal scaling functions -Rpbc.n.tuu, -Ppbc.n.seod, and A^pbc.n.aood of the Manna univer- 
sality class. In case of activity spreading, the simulations are started from a natural configuration of 
inactive particles. System sizes L — 512, 1024, . . . , 8192 are considered for D — 1, L = 64, 128, 256, 512 
for D — 2, and L — 16, 32, 64, 128 for D = 3. The dashed lines correspond to the power-law behavior 
of the infinite system x~°', , and , respectively. 



As can be seen, the above scaling law is fulfilled within the error bars. 

The next hyperscaling law connects the field exponent a with the correlation ex- 
ponents. Taking into consideration that a weak field may trigger spreading events the 
field exponent is given by (see e.g. ^7, or Eq. (|A.7|) 1 

a = Dv^ + v^^ -^\\^- (4-15) 

Again, this scaling law is fulfilled within the error bars (see FigureEO]). 

In case of directed percolation the rapidity reversal [Eq. 1)3. 12() ] implies a = 6 and 
therefore (3 = (3' [Eq. ()3. 127(1 ]. In a more general context one expects that both expo- 
nents differ [Jj, for example in our case where a conserved field couples to the order 
parameter. The number of independent critical exponents is therefore expected to be 
four (/3, /?', i^x) instead of three independent exponents for directed percolation (/?, 
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Figure 49: Test of certain scaling laws for the Manna universality class. The hyperscaling laws 7' = 
Di'± — 2/3 (upper left), a — Du± + — V\\5 (upper right), 6 = D/ z — P/v\\ — 5 (lower right), as well as 
the equality P = [3' (lower left) are checked. All scaling laws are fulfilled within the error bars. 

, i'±)- In order to check this scenario we compare the order parameter exponent f3 with 
in Figure EUl Surprisingly both values agree within the error-bars for all dimensions, 
suggesting (3 = (3' jl52j . It is worth mentioning that /? = /?' is not restricted to systems 
obeying the rapidity reversal. For example, the universality class of dynamical perco- 
lation is characterized by /3 = /?' although it violates Eq. (|3.10|) (see [HI I315| I176j and 
references therein). In the following /3 = /?' is assumed but we note that an doubtless 
verification of that equality requires direct measurements of the percolation probability. 
This remains the topic of future research. 
Finally, we consider the hyperscaling law 

8 D 
9 + ^ + 5 = —. (4.16) 
z^ii z 

Violations as well as validations of the hyperscaling law are reported in the literature 
(e.g. |285[I152] ). This confusion is caused by the fact that different forms of the hyper- 
scaling law [Eqs. 1)3.128113. 130|) ] are checked in different studies. Here, the numerically 
obtained exponents are plotted in FigureEHl according to Eq. (|4.16|) . As can be seen, 
the above scaling law is fulfilled within the error bars. 

The latter result is confirmed by independent measurements of the crossover from 
the initial particle growth to the asymptotic order parameter decay. In section 13.3.21 
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Figure 50: The universal scaling function -Rpbc.n (0, 0, oo, 1, x) describes the crossover from initial particle 
growth N^&t^ to the asymptotic order parameter decay patxt"". The data are obtained from simu- 
lations of the two-dimensional conserved threshold transfer process (CTTP). Different initial particle 
densities Pa.o = n/L^ with n = 1, 2, 5, 10, 20, ... , 20000 for L = 1024, 2048 and n = 160000 for L = 4096 
are used. Data that are affected by finite-size effects {0{t) > tpss) are skipped for a better visualization. 
The long-dashed lines correspond to the asymptotic scaling behavior of the universal scaling function 
Eq. 113.13211 . The dat a are scrutinized using the effective exponent (left) d\nx In _R(0, 0, oo, 1, a;). For 
Pa,o 1 systematic deviations from the expected scaling behavior (long dashed line, see Eq. I|3.132|l 'l 
occur. The largest deviations occur for p^ o ~ 0.02 and the corresponding curve is marked by a dashed 
line. 

we have seen that Eq. H4.16() is related to the scaling form 

p,{5p,h,L,t,p,^o) ~ (4.17) 

The hyperscaling law is valid if the scaling form is fulfilled for p^.o ~^ 0. On the other 
hand a = /3/f|| is obtained if the scaling form holds for p^.o ~^ 1- The scaling form 
can be checked by performing simulations with various initial densities p^.o- There- 
fore, an absorbing state is prepared at criticality by the natural dynamics. Then n 
different active seeds are created by n randomly selected particle movements on the 
lattice. Results of simulations of the two-dimensional conserved threshold transfer pro- 
cess are presented in Fieure lSUl As can be seen, the above scaling form is fulfilled for 
Pa.o ~^ 0) confirming that the hyperscaling law is fulfilled. But significant deviations 
from the expected scaling behavior of the universal scaling function i?pbc,n(0, 0, cxd, 1, x) 
[Eq. H3.132|) ] occur for x ^ 1. In particular, deviations are observed for y > 0.001. 
Due to the violation of the scaling form for increasing initial densities, the scaling law 
a = does not hold. This implies the violation of Eq. (|3.129|) . as observed in |285j . 
In conclusion, the hyperscaling law 14.161 holds for directed percolation as well as for 
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the Manna universahty class. But the universahty class of directed percolation fulfills 

a = ^ = ^ = 6 (4.18) 

whereas the Manna universality class is characterized by 

a / A = ^ = 5. (4.19) 

So far this scaling anomaly of the Manna class is not understood (see also jl59j l. i.e., no 
scaling law is known involving the order parameter decay exponent a. Thus the number 
of independent exponents of the Manna universality class is four if one assumes that 
13 = 13' [T^ . 



4.4 Self-organized criticality 

The concept of self-organized criticality (SOC) refers to driven-dissipative systems 
that naturally evolve to a critical state, characterized by power-law distributions of 
relaxation events jl()5| 110^ ll()7| The self-organization to the critical point dis- 

tinguishes SOC from ordinary critical phenomena, e.g. continuous phase transitions 
where a temperature-like control parameter has to be fine tuned to the critical value. 
At criticality the systems jumps among absorbing configurations via avalanche- like re- 
laxation processes. In this case the scale invariance of the system is reflected by an 
algebraic behavior of certain distribution functions, characterizing e.g. the avalanche 
size and duration. Besides of systems characterized by an extremal dynamics (see 
e.g. |316j l. the paradigms of SOC are the well-known sandpile models as the Bak-Tang- 
Wiesenfeld model (BTW) |in5j or the Manna sandpile model |290) . In particular the 
BTW model is analytically tractable and its theoretical understanding is related to its 
Abelian structure |317j . This allows to classify stable configurations as either transient 
or recurrent. The recurrent configurations can be represented by spanning trees j318j 
and many steady state properties are calculated exactly, such as the height probabili- 
ties or correlations |319| I320l I321j . On the other hand, the dynamical properties, i.e., 
the relaxation events which are often called avalanches, are not fully understood. A 
substantial progress was achieved by decomposing the avalanches into a sequence of 
relaxation events called waves |322[ I323[ I324j . Due to a spanning tree representation of 
waves and due to the equivalence between spanning trees and loop-erased random walks 
exact results of the wave exponents are available |236j . Furthermore, the intersection 
probabilities of random walks and loop-erased random walks allow to determine the up- 
per critical dimension = A |325j . confirming a numerical estimate |235j . But while 
the behavior of single waves is simple and well understood the behavior of avalanches 
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is complex and badly understood. The avalanche exponents are still not known ana- 
lytically for D = 2 and reliable numerical estimates are difficult to obtain due to a lack 
of simple finite-size scaling. Thus the interpretation of the numerical data is highly 
controversial (see e.g. IHTTIIH^ 1 . 



In contrast to the BTW model the Manna model was investigated successfully by 
numerical simulations (see for example [223 EHU IS2H1 EHSl EHS EHHl EHSl EHU M) 

but no exact results are known so far. A complete analytical solution of sandpile 
models is desirable because our understanding of SOC is still limited. In particular, 
such a solution would present a deeper insight into the fundamental mechanism of the 
self-organization to the critical point. So far it is recognized phenomenologically that 
the critical stationary state is reached through the balance of driving and dissipation 
(e.g. [3311 13321 1333j ). As pointed out by Grinstein, the critical behavior is approached 
in the slow driving limit, the so-called separation of time scales |331j . An important 
advance in the understanding of sandpile models was achieved by a mapping of sand- 
pile models to absorbing phase transitions, i.e., to a conventional critical phenomenon. 
This mapping was first discussed within a mean field approach in j334j and elabo- 
rated in a series of works [3351 12861 13331 . Both, the Manna and the BTW model 
exhibit absorbing phase transitions in closed ensembles where the particle density is 
strictly conserved (often called fixed-energy sandpiles |286j ). But numerical investi- 
gations of the BTW model are notoriously difficult since it is trapped in limit-cycles 
due to the deterministic dynamics. For example, it is still open whether the BTW 
model exhibits an absorbing phase transition of first or second order |337j . On the 
other hand, the various representatives of the Manna universality are class suitable for 
numerical investigations due to the stochastic dynamics. Therefore, we focus again 
on the Manna universality class. In particular, we consider the dynamical mean field 
theory of the conserved threshold transfer process in order to elucidate the mechanism 
of the self-organization to criticality. Numerous mean field treatments of SOC systems 
were pubhshed, e.g. E3H1 ESS ElOl IMD EH IMl ESI ESI ESI ESS|- But most 
studies focus on the determination of the avalanche exponents (e.g. by a mapping to a 
branching process) and do not explain the self-organization to criticality. The following 
mean field theory presents a complete picture of SOC in sandpile models, explaining 
how the critical state is attracted by the dynamics as well as the criticality itself. Due 
to the mapping of sandpile models to absorbing phase transitions, the SOC avalanche 
exponents can be expressed in terms of ordinary critical exponents of second order 
phase transitions jl52j . like /3, fy, and v^. This will be discussed in section 14.4.21 
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4.4.1 Sandpile models and absorbing phase transitions 

The mean field behavior of the conserved threshold transfer process with global parti- 
cle conservation is considered in section 14.1.21 In the following we consider a driven- 
dissipative modification. In addition to the dynamical rules above we introduce first an 
external drive: A particle is added from outside with probability p at an empty or single 
occupied lattice site. Obviously this perturbation breaks the global particle conserva- 
tion. Second we incorporate particle dissipation which counteracts the external drive. 
For the sake of simplicity we assume that a given target site belongs to the (dissipative) 
boundary and the particle will leave the system (An = — 1). The corresponding rate 
equations lead to the coupled differential equations 

dtp. = p(p-2pj + pJ-l + 2(l-e)p-4(l-e)p. + (l-e)Vf] (4.20) 
+ h{p-p,f, 

dtp = p{l-p.) -2ep.. (4.21) 

The latter equation reflects the competition between the particle flux in [p{l — p.)] and 
out (2epa) of the system. Global particle conservation p = const as well as Eq. 1)4. 5|1 are 
recovered for p = e = 0. For p = 0, e > 0, and h > the system tends to the empty 
lattice {p = and p. = 0) whereas it tends to the fully occupied lattice (p = 2 and 
= 1) for p > 0, e = 0. 

The steady state solution of the driven-dissipative system at zero field is given by 
_ P _ 2p2 + 6ep + e2 (2 -p)_ ^^^^^ 



2e + p' (2e + p)2 

Thus the absorbing phase {p. = 0) is not a steady state solution of the driven-dissipative 
system and the control parameter is now lost as an independent quantity. A linear 
stability analysis reveals that the above solution is an attractor of the dynamics for all 
values of e and p. In the slow driving limit {p <^ 1) the steady state solution Eq. (|4.22l) 
can be expressed in terms of the drive-dissipation ratio k = p/e 

K 2 + 6k + 2^2 

^^ = TT7' (2 + .)2 • (4-^3) 

It turns out that the parameter k, which was phenomenologically introduced by Grin- 
stein jSHB, is the appropriate parameter to describe the self-organization to the critical 
point. Eliminating the parameter k one recovers Eq. (|4.7j) . i.e., the order parameter 
behavior in the active phase of the closed systems 



pM = 2 - V5 - 2p for p>l/2 



(4.24) 
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Figure 51: The driven-dissipative conserved threshold transfer process for D = 2. The order parameter 
behavior Pa(p) of the closed {p = e — 0) and open (p > 0,e > 0) ensemble are identical (left figure). Data 
for e = IQ-'^ and p = n 10 * with n = 9, 10, 20, 30, 40, . . . 180 are shown, obtained from simulations of 
L = 128, 256, . . . , 2048. The right figure shows the order parameter and the rescaled control parameter 
as a function of the Grinstein parameter k = p/e. As in the mean field approximation, the two- 
dimensional order parameter vanishes linearly if one approaches the transition point (k 0). As 
can be seen, the control parameter behaves as (p — pc)'' 6c k. 



In this way, the control parameter in the active phase p G [Pcj2] is mapped to the 
interval of the drive-dissipative ratio k £ [0,oo[. In the strong driving limit p ^ e 
(k oo) the systems tends to the fully occupied lattice {p^ = 1, p = 2), whereas we 
find in the weak driving limit p <C e (k — > 0) in leading order 

Pa = , p = p, + K, for K^O. (4.25) 

In the limit k — > 0, the driven-dissipative system tends to the critical point of the 
absorbing phase transition [p = p^, p^ = 0). 

The crucial role of the Grinstein parameter n is even valid beyond the mean field 
approach. Simulation data of the driven-dissipative CTTP in D = 2 are presented in 
Figure lKTl As can be seen, the data of the driven-dissipative CTTP recovers the active 
phase solution of the closed CTTP, i.e., for p = e = 0. Furthermore, the critical point 
is approached for k — > 0. In leading order, the two-dimensional system behaves as 

Pa = I , P = Pc + 1^^^^ , for K ^ . (4.26) 

Surprisingly, the two-dimensional order parameter and the mean field order parameter 
behave in the same way, i.e., they vanish linearly for k — > 0. For the sake of completeness 
we just mention that the corresponding fluctuations scale as Ap^ 6c k~'^ and Apoc k. 

Now, let us focus on the dynamical behavior close to the steady state solution for 
K^O. Using Eq. (IT^ and Eq. ^1^) we find 



dptp = 0{k) , d,tp = 0{k^) , 



(4.27) 
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i.e, the control parameter varies for «; ^ only on the perturbation scale Tp^^p = 
1/p whereas it can be considered approximately as constant on the dissipation scale 
"^diss = 1/e- Therefore we set p = const in the following. Furthermore, we assume 
that the system displays an avalanche- like response for infinitesimal driving {p — > 0), 
i.e., no further perturbations take place until an absorbing configuration {p^ = 0) is 
reached. This extreme separation of time scales (k — > 0) is necessary in order to 
identify individual relaxation events (avalanches). Only in this case it is justified to 
characterize the scale invariance via avalanche distribution functions Otherwise 
the scale invariance is manifested in terms of algebraic decaying correlation functions 
similar to usual phase transitions. 

Within the approximations discussed above, the avalanche processes are determined 
by Eq. (|4.2()|) with p = 0. At criticality {p = p^) we find (neglecting again ©(pf )) 

pM = 7 • (4.28) 

(l + 4 if^ p,,o) e^* - 4 if^ p,,o 

For sufficiently small values of e the order parameter decays as 

P-M ~ -1 , — for et«l, (4.29) 

Pa,o + 4(1 - e)t 

4(1 - e) 

Thus close to the critical point (k 0) an avalanche, caused by the perturbation 
Pa.oi decays algebraically before an exponential cutoff occurs. This cutoff is caused by 
the particle dissipation and takes place at r^utofi = A pure power-law behavior is 
obtained in the limit e ^ only. Therefore the driven-dissipative CTTP self-organizes 
itself to the critical point and exhibits scale invariant avalanches in the SOC limit k — > 
and e ^ p2ll333] . 

It is customary in the literature to consider the response function 



_ dp^ 

~ dp 



2e 



that should be distinguished from the susceptibility x [Eq. (|2.14j) ]. It is often interpreted 
that the singular behavior of Xv for e ^ signals that criticality is approached in the 
SOC limit (k ^ 0, e — > 0). But this is a misleading conclusion as can be seen if one 
considers the more appropriate form 

= ^ = —777- — To- (4-32) 

op e (2 + kY 

The response function Xv exhibits for e ^ a singular behavior in the whole active 
phase (k > 0), not only at the critical point (k = 0). Thus the singular behavior of Xv 
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does not reflect the critical behavior at p^. It just reflects that an infinite response is 
obtained if one adds particles to a system with particle conservation (e = 0). 

It is worth to consider how the conjugated field affects the driven-dissipative system. 
The steady state solution of Eq. (|4.2()jl yields 

lim p{k, h = 0) = p,, lim p{k, h > 0) = . (4.33) 

K— >0 

The external field drives the system away from the critical point, i.e., SOC is destroyed 
if one applies an external field that is conjugated to the order parameter. 

In summary, we have seen that the CTTP undergoes for p = e = an absorbing 
phase transition from an active to an absorbing phase. The situation is completely 
changed if one considers an open ensemble where fluctuations of the particle density 
occur. In contrast to the closed ensemble {p = e = 0) the absorbing phase is no longer 
a solution of the driven-dissipative system (k > 0). The remaining solution of the open 
ensemble is exactly the active phase solution of the closed ensemble. It is the key to 
the understanding of SOC that the limit k ^ tunes the control parameter to the 
boundary of the active phase, i.e., to the critical point. This explains why the trivial 
limit K ^ in the open system corresponds, in the closed system, to the fine tuning of 
the control parameter to the critical point of the absorbing phase transition {p p^). 
Remarkably, this limit occurs frequently in driven-dissipative systems in nature, i.e., 
many physical processes are characterized by a large separation of times scales that 
makes k extremely small |1()7| I331| . For instance the motion of the tectonic plates that 
causes earthquakes is extremely slow compared to the time scale on which earthquakes 
proceed (see |331j ). Thus, the dissipation and driving rate are separated by many 
orders of magnitude. Although the drive-dissipation rate remains finite the system is 
very close to the critical point (k = 0) and displays a power-law behavior over several 
decades, namely the well known Gutenberg-Richter law |346j . 

4.4.2 Avalanche exponents 

In the previous section we discussed how SOC is mapped to an absorbing phase transi- 
tion. Usually self-organized critical systems are characterized by avalanche exponents 
whereas absorbing phase transitions are described by the critical exponents /3, /?', v^, 
and v±. Thus it remains to show how the avalanche exponents are connected to the 
exponents of the second order phase transition |152j . 

In the critical state of SOC systems the external driving p triggers scale invari- 
ant avalanche-like relaxation events. In particular the avalanches are characterized by 
several quantities (see for example |29H I3()8| ). e.g. the size s (number of elementary 
relaxation events), the area a (number of distinct relaxed sites), the time t (number 
of parallel updates until the configuration is stable), as well as the radius exponent r 
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(radius of gyration). In the critical steady state the corresponding probability distri- 
butions decay algebraically 

oc , (4.34) 

characterized by the avalanche exponents with x G {s,a,t,r}. Assuming that the 
size, area, etc. scale as a power of each other, 

X 6c x'^'--' (4.35) 

the scaling laws 

7..' = (4.36) 

- 1 

are obtained. The exponent 7t,. equals the dynamical exponent z, the exponent 
corresponds to the fractal dimension of the avalanches and the exponent 7^;, indi- 
cates whether multiple relaxations of a lattice site are relevant (7^,^ > 1) or irrelevant 

{%a = 1). 

These avalanche exponents are connected to the spreading exponents 5, 0, and z 
(see e.g. |3n2( 13151 1152) 1. First the survival probability -Pa(t) is simply given by the 
integrated avalanche duration 

m = E P^i^') (4-37) 
t'=t 

yielding 

Tt = 1 + 6. (4.38) 
Since jtr = z the radius exponent is given by 

= I + z5. (4.39) 

Taking into account that the avalanches of the Manna model are compact (7^,^ = D) 
below the upper critical dimension = 4 |3()8| 1309 j we find 

r. = 1 + ^ . (4.40) 

Finally the number of topplings St for an avalanche that is active at time t equals the 
integrated numbers of active sites, i.e., 

t 

s.P^it) = E ^a(t') (4.41) 
t'=o 

leading to |315j 

Ts = 1+ (4.42) 
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Figure 52: The avalanche exponents of the Manna model in various dimensions. In order to avoid 
overlaps the exponents are slightly shifted. The avalanche exponents of the SOC version of the Manna 
model (left) are obtained from glZl for D = 1 and from PTOII^ for D = 2, 3. Using the Eqs. iflTW 
14.42^ the avalanche exponents (middle) are obtained from the spreading exponents 5, 9, and z. Using 
certain hyperscaling relations it is possible to express the avalanche exponents (right) in terms of the 
exponents of the continuous absorbing phase transition (/3, /3', i'±, //y), see Eqs. 14. 4314. 461 . 



Thus the avalanche exponents of sandpile models are naturally related to the spreading 
exponents of absorbing phase transitions. Due to the valid hyperscaling law Eq. (|3.128)1 
and Eq. (|3.126j) the SOC avalanche exponents can be expressed in terms of the ordinary 
critical exponents of the corresponding absorbing phase transition [152| 

0' 

= 1 + — , (4.43) 
0' 

n = 1 + -, (4.44) 
^11 

r. = 1 + 7^, (4.45) 

Ts = 1 + ^ ^ n . • (4.46) 

i^ll + v^D — p 

In Fig. 1^21 we compare these values to the avalanche exponents obtained from SOC 
simulations of the Manna model j3()9| I29.'t1 I347j . All SOC exponents agree within the 
error-bars with the avalanche exponents derived via the above scaling laws. Thus it is 
possible to express the avalanche exponents (r^, t^,. . .) of SOC systems in terms of the 
usual critical exponents of a second order phase transition 0', I'j^, i/y). In this way, 
the critical state of SOC sandpile models is closely related to the critical state of an 
ordinary second order phase transition. 
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Crossovers between different universality classes are well known from equilibrium phase 
transitions. Classical examples are ferromagnetic systems exhibiting weak uniaxial spin 
anisotropy (see for example |17 | I348j ). Approaching the transition point Heisenberg 
critical behavior is observed. But very close to a crossover to Ising critical behavior 
occurs. In the topological language of renormalization group theory crossover effects 
may occur if more than one fixed point is embedded in the critical surface. In case of 
the example above, the crossover corresponds to a trajectory which passes close to the 
Heisenberg fixed point but is eventually driven away to the Ising fixed point. In terms of 
scaling forms, crossover phenomena are described by additional relevant scaling fields, 
characterized by a so-called crossover exponent </> |,S49j . 

Although crossover phenomena are well understood in terms of competing fixed 
points, several aspects are still open and are discussed in the literature. For example, the 
question whether the so-called effective exponents |156j fulfill certain scaling laws over 
the entire crossover region was revisited several times (see e.g. |35Ul 13511 13521 13531 1354j 
and references therein). This question is related closely to the general and important 
question whether effective exponents obey the scaling laws at all. For example, it is 
known experimentally j355j as well as theoretically [7S] that the asymptotic scaling 
behavior is often masked by corrections to scaling, so-called confluent singularities (see 
section EH- In this case it is useful to analyze the data in terms of effective exponents 
and the question above arises naturally . The validity of the scaling laws for effective 
exponents was often concluded from experimental (see e.g. j35Uj and references therein) 
and numerical investigations (see for example j353j ). In contrast, renormalization group 
approaches predicted a violation of the scaling laws for effective exponents ;.350^. This 
renormalization group result is confirmed by numerical investigations of the Ising model 
indicating a violation of the Rushbrook scaling law [352j . 

Another question of interest concerns the universality of the crossover scaling func- 
tions. The range where the universal critical scaling behavior applies is usually re- 
stricted to a small vicinity around the critical point. Therefore it is questioned whether 
the full crossover region that spans usually several decades in temperature or conjugated 
field, can be described in terms of universal scaling functions [3561 13571 13511 0551 13591 
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I352| I36n[ I36H I362| I3n3[ I363| I354j . Renormalization group approaches predicted a non- 
universal behavior if one uses finite cutoff lengths, whereas infinite cutoff lengths (which 
corresponds to unphysical vanishing molecular sizes) lead to a universal crossover scal- 
ing (see e.g. j356| 1357] ). On the other hand, also the experimental situation is unclear 
since measurements over the whole crossover region are difficult, and accurate results 
are rare (see |351| for a discussion). Thus several attempts were performed in order to 
address this question via numerical simulations [3511 13521 13591 13581 13U31 13621 1354j . But 
only one model, i.e., only one representative of a given universality class was considered 
in previous studies. Analyses that incorporate different models belonging to the same 
universality class were performed recently j3L)3ll35^ . These studies are reviewed in the 
following. 



5.1 Crossover to mean field scaling behavior 

The critical behavior of a system exhibiting a second order phase transition with non- 
mean field scaling behavior is strongly affected by the range of interactions. The longer 
the range of interactions the stronger will be the reduction of the critical fluctuations, 
and in the limit of infinite interactions the system is characterized by the mean field 
scaling behavior. As stated by the well known Ginzburg criterion j364[ I365j . mean 
field-like behavior occurs even for finite interaction ranges sufficiently far away from 
the critical point. A crossover to the non-mean field scaling behavior takes place if one 
approaches the transition point. This crossover is described by a crossover exponent (/> 
which is known exactly due to the Ginzburg criterion. Therefore the crossover from 
mean field to non-mean scaling is theoretically and numerically (see e.g. |35H I358| 
Ens ESI EHOl ESI EHl EOSl EMI ESU) well suited to examine if universality holds for 
crossovers. Following the spirit of this work, we focus on non-equilibrium crossover 
phenomena and consider the crossover scenario in the Manna universality class (in 
the so-called critical crossover limit |35H 136(1] ). But it is worth mentioning that the 
results can be applied to continuous phase transitions in general, including equilibrium 
crossover phenomena. 

According to the definitions of the CLG model, the GTTP, and the Manna model 
(see section IITT]) . particles of active sites are moved to nearest neighbors only, i.e., the 
range of interactions is i? = 1. It is straightforward to implement various ranges of 
interactions into these models |3U3j . In these modified models particles of active sites 
are moved (according to the rules of each model) to sites randomly selected within a 
radius R. The dynamics of the models considered are characterized by simple particle 
hopping processes, i.e., various interaction ranges can be implemented easily without 
affecting the simulation performance. This is an advantage compared to e.g. equilib- 
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Figure 53: The order parameter (left) and the order parameter fluctuations (right) of the two- 
dimensional conserved threshold transfer process (CTTP) for various values of the interaction range R. 
With increasing range the critical density tends to the mean field value p^.R^oo = 1/2. As can be seen, 
the power-law behavior of and Apa changes with increasing range of interactions. The data are 
obtained from simulations of various systems sizes L < 2048. 



rium systems like the Ising model where increasing interaction ranges slow down the 
dynamics. A disadvantage is that for zero field the crossover occurs only in the active 
phase, whereas equilibrium systems can be investigated below and above the phase 
transition |351l I36U] . 

Notice that the simulations have to be performed in the so-called critical crossover 
limit, i.e., the corresponding correlation length is sufficiently large in the whole crossover 
region j358l 13511 l36Uj . For any finite interaction range R the phase transition is charac- 
terized by non-mean field scaling behavior, whereas mean field scaling occurs for infinite 
interactions. The order parameter as well as the order parameter fluctuations Ap^, 
of the two-dimensional CTTP are plotted in Figure|S31 for various interaction ranges 
at zero field. As can be seen, the transition point p^ j^ depends on the range of the 
interactions and the power-law behavior of and /S.p^ changes with increasing R. 

Before examining how the increasing range of interactions affects the scaling be- 
havior it is inevitable to mention that the parameter R is not an appropriate quantity 
to describe the scaling behavior since it describes the maximum range of interactions. 
Thus the actual strength of the interactions may depend on the lattice structure. To 
avoid these lattice effects one uses the effective interaction range |359j 

^cff = Z:/, for - r.\<R (5.1) 

where denotes the number of lattice sites within a radius R. The values of the 
effective interaction ranges are listed for two- and three dimensional simple cubic lattices 
in TableEland Tabled 
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Approaching the transition point from the active phase non-mean field critical be- 
havior is observed for all finite values of R. But mean field behavior occurs away 
from the critical point if the long range interactions reduce the critical fluctuations 
sufficiently. The crossover between these two scaling regimes is described by the well 
known Ginzburg criterion |3641 1365j which states that the mean field approximation 
is self-consistent in the active phase as long as the fluctuations within a correlation 
volume are small compared to the order parameter itself 

r^Ap. « pI (5.2) 

In the steady state of absorbing phase transitions ^ corresponds to the spatial correla- 
tion length ^ that diverges at the critical point according to 

^6ci?,ff (p-p,,^)-''^.^^-. (5.3) 

Thus, despite of metric factors the mean field theory applies at zero field if 

1 <C {p - p^_^)^mf+2/3mf-x,mf^ , (5.4) 

where the usual power-laws 

Pa ex (p - Pcfi)^^^ , Ap, 6c (p - p,k)-^mf , (5.5) 

and Eq. (|5.3() are used. Obviously, Eq. 1)5. 4|) is fulfilled for p — > p^j^ if 

D > ^MF + 2/3mf ^ (5_g) 

otherwise mean field theory fails to describe the critical behavior. The marginal case 
corresponds to the upper critical dimension 

= ^MF + 2/3mf _ 

In this way the concept of an upper critical dimension is introduced phenomenologically 
by the Ginzburg criterion, independent of the renormalization group scenario. In case 
of directed percolation and the Manna universality class the mean field values /?mf = Ij 
z^±,MF = 1/2, and 7mf = lead to Z?^ = 4 as expected. It is convenient to introduce the 
exponent 

-De - -D 4 - D , , 

(p = z^±,MF = • (5.8) 

Thus mean field scaling behavior occurs as long as 1 <C R^s {p — Pc,r)'^ whereas non- 
mean field scaling behavior is observed for 1 ^ R^^ (p — Pc,r)^- The point at which 

O (P - Pc,^)'^) = 1 (5.9) 
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Figure 54: The order parameter data of the two-dimensional conserved threshold transfer process 
(CTTP). The data are rescaled according to Eq. 15.1311 . The dashed line corresponds to a power-law 
behavior of the two-dimensional order parameter with exponent /3 = Pd=2- 



is usually considered as the crossover point and (p is termed the crossover exponent |349j . 

In order to describe the crossover behavior the scaling forms have to incorporate 
the range of interactions as an additional scaling field. For example, the non-universal 
scaling form of the order parameter is given by 

pMh,R.f,) ~ X~''^^i{X{p-p^,n),hX''^\R-^X^), (5.10) 

where the scaling power of R^g is motivated by Eq. (|5.9|) . In the following, r and 9^ 
are used to denote the non-universal and universal scaling functions, describing the 
particular crossover. Obviously, the added scaling field is relevant (0 > 0) below the 
upper critical dimension whereas it is irrelevant (0 < 0) above D^. The mean field 
scaling behavior should be recovered for — > oo 

r(x,0,0) 6c x^MF foj. x^O. (5.11) 

On the other hand, the D-dimensional order parameter behavior is observed for finite 
interaction ranges 

r(x,0,i2-^) cx for x^O. (5.12) 

It is instructive to consider the amplitudes of the corresponding power-laws because 
the amplitudes display a singular dependence on the range of interactions. Setting 
R~ffX'^ = 1 in Eq. (|5.1U() yields sufficiently close to the critical point 

p.(p,o,i?eff) ~ i?:/-^/"^ t(i?L/*(p-p„,),o,i) 

6c #-'3mf)/'/' ip-p.,nf''. (5.13) 
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Figure 55: The order parameter of the two-dimensional (left, (j) — 1/2) and three dimensional (right, 
</> — 1/6) conserved threshold transfer process (CTTP) and conserved lattice gas model for various 
values of the interaction range R. The data are rescaled according to Eq. I|5.10^ . The dashed line 
corresponds to the power-law behavior of the D-dimensional models and to the mean field scaling 
behavior, respectively. The data are obtained from simulations of various system sizes {L < 2048 for 
D = 2 and L < 128 for D = 3). 



Thus this power-law behavior and the corresponding power-laws of the fluctuations, 
the susceptibility, etc. are valid only for finite interaction ranges whereas they become 
useless for infinite R, signalling the change in the universality class. In this way the sin- 
gular amplitudes reflect the breakdown of a scaling regime if a crossover is approached. 
The critical amplitude behavior can be observed in simulations. The corresponding 
data of the two-dimensional conserved threshold transfer process is presented in Fig- 
ure l54l The rescaled order parameter data tend to the same power-law behavior if one 
approaches the transition point. 

Furthermore, data obtained from simulations at zero field and various interaction 
ranges have to collapse onto a single curve by plotting R^^^^'^ as a function of 
{p — Pc,h) Rts' ■ The respective plots of the two- and three-dimensional CTTP and 
the CLG model are shown in Figure l55l A good collapse is observed over the entire 
crossover which spans roughly 8 decades. The asymptotic power-law behavior of the 
D-dimensional scaling behavior as well as of the mean field scaling behavior are recov- 
ered well. But as can be seen for D = 2, each model is characterized by its own scaling 
function r. In order to obtain the universal crossover scaling function non-universal 
metric factors have to be introduced. This will be done in the next section where we 
focus on the control parameter driven crossover (temperature driven crossover in equi- 
librium). The subsequent section is devoted to the field driven crossover. In particular, 
the data of the field driven crossover allows to check the validity of the Widom scaling 
law over the entire crossover region. 
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5.1.1 Control parameter driven crossover 

In the following we address the question whether the crossover scaling functions, which 
spans several decades, can be described in terms of universal scaling functions. There- 
fore, non-universal metric factors have to be incorporated in the scaling forms, i.e., we 
make the phenomenological ansatz 



p^{p,h,R^s) ~ A-^^^- ^(a,(/9-p,,«)A,a,/iA"^^a«i^;fllA'^), (5.14) 
The mean field scaling function should be recovered for i? — > oo, thus 

^(x,y,0) = RMF{x,y) (5.15) 

and, therefore, 

^(1,0,0) = i?MF(l,0) = 1, ^(0,1,0) = i?MF(0,l) = 1, (5.16) 

which imphes 

Op = , 0/1 = ah.fl^oo • (5- 17) 

These metric factors were already determined in independent studies where absorbing 
phase transitions with infinite particle hopping were investigated l,S()5j . 

The non- universal metric factor has to be determined by a third condition. 
Several ways are possible (e.g. 9^(0, 0, 1) = 1) but it is convenient to enforce that 9^ 
scales as 

$n(3;, 0, 1) ~ x^^ , for x^O. (5.18) 
Setting a^^R^fi^ A* = 1 in Eq. (IQHl yields at zero field 

pM 0, R.s) ~ (a« ii.ft )-^^^/'^ ^{a,{p - p^,n) a"!"^ RIL^ , 0, 1) . (5.19) 

Taking into account that the D-dimensional scaling behavior is recovered for i? = 1 we 
find 

P^apj^\ _ ^^20) 

According to the scaling form above, we plot in FigurelSBl the rescaled order parameter 
as a function of the rescaled control parameter for the two-dimensional {(p = 1/2) 
models. The data of the metric factors are obtained from directed measurements of the 
corresponding amplitudes (see |141j ). An excellent data-collapse is observed over the 
entire range of the crossover confirming the phenomenological ansatz. 

Since the full crossover region covers several decades it could be difficult to ob- 
serve small but systematic differences between the scaling functions of both models. 
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Figure 56: The universal crossover scaling function of order parameter at zero field for D — 2 (left). 
The metric factors are given hy c — OpO^ and d — a%. The data of all models display an excellent 
collapse onto the universal crossover scaling function ^R{x,0,l). The dashed lines correspond to the 
asymptotic behavior of the two-dimensional system {Pd=2 ~ 0.639) and of the mean field behavior 
(/3mf = 1). The right figure displays the corresponding effective exponent /3eff. 



Therefore, it is customary to scrutinize the crossover via the so-called effective expo- 
nents nsni Esa M M 

/Soft = li In ^(x, 0,1). (5.21) 

omx 

The corresponding data are shown in FigureEHl The effective exponent changes mono- 
tonically from its mean field value to its two-dimensional value. Non-monotonic crossover 
behavior may also occur, e.g. in equilibrium systems j^-i51j and even in the Manna uni- 
versality class |296j . The excellent collapse over more than 7 decades supports strongly 
the hypothesis that the crossover function is a universal function. It is worth mention- 
ing that the collapse includes the data for small interaction ranges, in particular R = 1 
and R = 2. Renormalization group analyses of the same crossover scenario in spin 
models reveal non-universal corrections to the universal scaling functions |361j . These 
corrections scales as R~^ (in two loop order), i.e., significant corrections are expected 
to occur for small interaction ranges R. 

We now consider the order parameter fluctuations. Analogous to the order param- 
eter we use the scaling form 

0^ Ap^ip, h, i?,ff) ~ A^MF 2)(a^ {p - p^ ^) A, a, h X^^^,a~^ i?"! ) . (5.22) 

Again the mean field behavior should be recovered for R ^ oo, implying 

S)(x,y,0) = D^Ax,y) (5.23) 

and thus 

2)(1,0,0) = ^mf(1,0) = 2, 3(0,1,0) = 5mf(0,1) = 1, (5.24) 
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Figure 57: The universal crossover scaling function of the order parameter fluctuations at zero field for 
D = 2 (left). The metric factor is given by c = apa%. The data of both models display a good collapse 
to the universal crossover scaling function 'Z>{x,0, 1). The dashed lines correspond to the asymptotic 
behavior of the two-dimensional system (7^=2 = 0.381) and of the mean field behavior (7mf = 0). 
The universal amplitude is given by ttia.p = 1.28. The right figure displays the corresponding effective 
exponent 7^ff. 



as well as 

Setting a-^R-g^ A"^ = 1 in Eq. SKl^ yields at zero field and for 'yl^j, = 
Ap,(p, 0, ~ 2)(a, (p - aJ'^Rll'^, 0, 1). 



(5.25) 



(5.26) 



The fluctuations diverge at the critical point, i.e., the universal function D scales as 

2)(j;,0,l) ~ mA,pX"'^D , for x^O. (5.27) 
The power-law amplitude m^^p is determined in the following way j3()3j : The scaling 



(5.28) 



form Eq. (|5.22)) has to equal for R = 1 the D-dimensional scaling behavior 

(p-/Oc,fl=l)\"^^ 



Apa ~ a^^^^, Do (1,0) iap^ii=i 



The value of the universal scaling function D£,(1,0) can be obtained from a direct 
measurement of the corresponding Z)-dimensional system (see TableOJ and the universal 
amplitude of the crossover fluctuation function is therefore given by 



m 



\ 7d/3mf/(/3mf-/3d) 



(5.29) 



According to the scaling form Eq. (|5.2(i)) we plot in Figure !^ the rescaled fluctuations 
as a function of the rescaled control parameter for the two-dimensional CLG model 
as well as for the CTTP. We observe again a good data-collapse of the data over the 
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entire region of the crossover. Furthermore, the asymptotic behavior is recovered for 
both power-laws, confirming the scaling ansatz Eq. H5.22() . The corresponding effective 
exponent 

= -|j^ln2)(x,0,l) (5.30) 

is displayed in Figure lFTl Although the data of the effective exponents are suffering 
from statistical fluctuations one can see that both models are characterized by the 
same universal behavior. 

5.1.2 Field driven crossover: Violation of Widom scaling law 

In this section we focus our attention on the field driven crossover, i.e., we consider the 
scaling behavior at the critical density p^ ji as a function of the conjugated field. The 
scaling form at the critical point is given by |354j 

~ (a«i?.«.)-^MF/0 ^(o,a,/iar"/'^i?:r/^l). (5.31) 

For vanishing field the universal function scales as 

^(0, X, 1) ~ m,,ft x^^/"^^ , for x^O. (5.32) 

The universal amplitude m^^h can be determined similar to the fluctuation ampli- 
tude mA,p, yielding 

We plot in FigureEHl the rescaled order parameter {auR^sf' as a function of the 
rescaled field ahh^anR^s)'^ . Again the data cover the entire crossover region and both 
asymptotic power-laws are clearly recovered. The data-collapse of the effective exponent 

^ ln$n(0,x,l) (5.34) 



o-/eff 9 In: 

confirms again the universality of the crossover scaling function fR. 

Next the order parameter susceptibility is considered. The scaling form of the 
susceptibility is given by |354j 

X(P, h, i?,ff) ~ A^^F X{a,{p - p,,h) a, o, h A"^l^ o-^iJ-^ A*) . (5.35) 

The mean field behavior is recovered for R ^ oo, i.e., 

X(x,y,0) = XMFix,y), (5.36) 

implying 

X(1,0,0) = Xmf(1,0) = 1, X(0,1,0) = Xmf(0,1) = ^. (5.37) 
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Figure 58: The universal crossover scaling function of the order parameter at the critical density for 
D = 2 (left). The metric factors are given by e = 0^0^ and d = a%. The dashed lines correspond 
to the asymptotic behavior of the two-dimensional system {(50=2/(^0=2 = 0.287) and of the mean field 
behavior (/3mf/o"mf = 1/2). The universal amplitude is given by m^,h = 0.681. The right figure displays 
the corresponding effective exponent {f}/a)^([. 

In order to obtain a data-collapse we set a^^R^g A*^ = 1 in Eq. (|5.35|1 . yielding at the 
critical density 

c^h' x(Pc,«, h, - (a«i?,ff)^^'^/<^ X(0, a, h {aM'''^'^/^, 1) . (5.38) 
Approaching the transition point the susceptibility is expected to scale as 

X(0,x, 1) ~ m^,fta;-T°/'^^, for x^O. (5.39) 
On the other hand, the susceptibility obeys for R = 1 

X ~ a,,n=i X^(0,1) (a,,«.i/i)-^^/'^« . (5.40) 
Thus the universal power-law amplitude is given by 



a 



„-7mf/<P+o"mf7d/o-d0 I-'d _ HP 



m^,^ = ^-^iM^,'P^(^MviD,^D9 s:^ ^ m,,,^, (5.41) 

where Eq. (|3.52j) is used. The rescaled susceptibility is shown in Figurel^ Over the 
entire crossover region we got an excellent data-collapse including both asymptotic 
scaling regimes. The right figure displays the effective exponent 

lnje(0,x,l) (5.42) 



aj^fi dlnx 



which exhibits again a monotonic crossover from the two-dimensional scaling regime to 
the mean field scaling behavior. 
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Figure 59: The universal crossover scaling function of the susceptibility at the critical density for 
D = 2 (left). The metric factors are given by e = aha^ and b — aj^^a%. The dashed lines correspond 
to the asymptotic behavior of the two-dimensional system {;y 0=2/(^0=2 ~ 0.713) and of the mean field 
behavior (7mf/o"mf = 1/2). The universal amplitude is given by m^^h = 0.208 which agrees within the 
error-bars with ma,hPn/o-n ~ 0.195±0.013 [see Eq. 115. 4H ]. The right figure displays the corresponding 
effective exponent {"f/a)^it. 



In this way we have obtained the effective exponents {P/a)^^ and {'y/a)^^ for the 
field driven crossover from mean field to non-mean field behavior. Thus we are able to 
check the corresponding Widom scaling law [Eq. (|3.5fl|) ] 

I) • f^'^^' 

for the whole crossover region. The corresponding data are shown in Figure l^ As 
can be seen the Widom scaling law is fulfilled for the asymptotic regimes [D = 2 and 
mean field scaling behavior) but it is clearly violated for the intermediate crossover 
region. This result is not surprising if one notices that the above Widom law Eq. (|5.43jl 
corresponds to the differential equation [see Eq. ()5.34|) and Eq. (|5.42l) ] 

5 In (9 ~ (9 In ~ 

- TH— TT ^(0' ^' 1) = 1 - 1^ ^(0' ^' 1) ■ (5-44) 
omx ox omx 

Using 1 = dlaax/ dlnx we get 

-ln5,^(0,x,l) = Inax - ln^(0,x, 1) + c, (5.45) 

where c is some constant. It is straightforward to show that this differential equation is 
solved by simple power-laws [^(0, x,l) = CqX^^ with Ci = 1/aexpc]. Thus the Widom 
scaling law is fulfilled in the asymptotic regimes only. In the case that the scaling 
behavior is affected by crossovers, confluent singularities, etc. no pure power-laws occur 
and the scaling laws do not hold for the corresponding effective exponents. 
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Figure 60: The violation of the Widom scahng law Eq. H5.43^ in the crossover regime. Data of the 
Manna model and the conserved threshold transfer process are plotted. The dashed lines correspond 
to the asymptotic regimes (two-dimensional behavior and mean field scaling behavior). 



In conclusion, the numerically investigated crossover from mean field to non-mean 
field scaling behavior within the Manna universality class reveals that the corresponding 
crossover scaling functions are universal. Although the systems considered exhibit non- 
equilibrium phase transitions, the results can be applied to continuous phase transitions 
in general, including equilibrium systems. All effective exponents considered change 
monotonically from their mean field values to their D-dimensional values. In contrast, 
the crossover from the stochastic Manna universality class to the universality class of 
deterministic particle distribution is characterized by a non-monotonic behavior |296j . 
Regardless of a monotonic or non-monotonic behavior, the crossover scaling functions 
are universal. This is of experimental relevance since measurements of critical behavior 
are performed often within a crossover region. 
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Similar to equilibrium critical phenomena, the great variety of non-equilibrium phase 
transitions can be grouped into different universality classes. Each universality class 
is characterized by a certain symmetry which is often masked within the Langevin 
equation approach, but it is expressed clearly within the corresponding path integral 
formulation (see e.g. |H1 1176j and references therein). For example, non-equilibrium 
critical systems belong to the directed percolation universality class if the associated 
absorbing phase transition is described by a single component order parameter and if 
the corresponding course grained system obeys the rapidity reversal symmetry asymp- 
totically. Different universality classes occur if the rapidity reversal is broken, e.g. by 
quenched disorder or by a coupling of the order parameter to a non-diffusive conserved 
field as in the case of the Manna universality class. Since no complete classification 
scheme is known so far, a phenomenological identification of non-equilibrium univer- 
sality classes is needed. Beyond the values of the critical exponents, universal scaling 
functions are a very sensitive and accurate tool to verify a system's universality class. 
Therefore, universal scaling plots reflect the robustness of a given universality class im- 
pressively. For example, a striking demonstration of the universality class of directed 
percolation is presented in Figure[^ where a universal data-collapse of five different 
models is shown. 

Additionally to the manifestation of universality classes, universal scaling functions 
are useful in order to check renormalization group results quantitatively. Since the 
renormalization group theory is the basis of our understanding of critical phenomena 
it is of fundamental interest to examine the obtained results. Due to the continuing 
improvement of computer hardware, accurate numerical data have become available 
in the last years, resulting in a fruitful and instructive interplay between numerical 
investigations and renormalization group analyzes. Incorporating the conjugated field, 
a series of opportunities is offered to compare renormalization group results to those 
of numerical investigations. Simulations performed for non-zero field include the mea- 
surements of the equation of state, of the susceptibility, as well as a modified finite-size 
scaling analysis appropriate for absorbing phase transitions. 

In this work, two non-equilibrium universality classes are considered in detail. 
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namely the directed percolation and the Manna universality class. Although steady 
state, dynamical, as well as finite-size scaling functions are determined future work is 
needed to complete the gallery of universal scaling functions. For example, the univer- 
sal scaling behavior of correlation functions has not been considered so far. Another 
example is the scaling behavior of persistence distributions which has attracted a lot 
of research interests in the recent years (see e.g. f^FTH IHHHl IHHHl 117701 lITm IITT^ r?7^ 
13751 13761 13771 13781 1379j ) . Persistence distributions are defined as the probability that 
a certain physical quantity does not change its state during the stochastic evolution. 
Various distributions are investigated, associated with the persistence behavior of the 
local as well as global order parameter. In particular, the global distributions are ex- 
pected to be universal, i.e., they are suitable to reflect the corresponding universality 
class. So far, most work has focused on the determination of the exponents, neglecting 
the corresponding universal scaling functions. 

Furthermore, it is interesting to consider crossover phenomena between the different 
universality classes. The threshold transfer process provides an opportunity to examine 
the crossover from directed percolation to the Manna universality class. Assuming that 
particle creation/annihilation processes occur with rate k, a crossover to the conserved 
threshold transfer process takes place for k — > 0. 

Within this work, we focus our attention on short range interacting systems. Long- 
range interactions may affect the scaling behavior significantly. It is known from equilib- 
rium that slowly decaying long-range interactions reduce the value of the upper critical 
dimension |2691 12701 12711 1272j . Since lower values of the upper critical dimension 
allow to simulate larger system sizes, it should be possible to investigate the logarithmic 
corrections at with a higher accuracy in that way. 

Additionally to the considered directed percolation and Manna universality class, it 
is desirable to determine the universal scaling functions of other non-equilibrium phase 
transitions. According to field theories as well as numerical investigations, different 
universality classes emerge from disorder effects 13 EHUl EHD EHl EHSl EHl • The 
different universal behavior is of particular interest since disorder effects are suspected 
to destroy possible experimental realizations of directed percolation (see e.g. |28()j ). 
If this explanation is valid, the disorder associated scaling behavior will also occur 
experimentally and universal scaling functions will be useful for an accurate identi- 
fication. Another established universality class is related to the parity conservation 
process |283| 128^ . So far, no systematic analysis of the corresponding universal scaling 
functions has been performed. 

Pinning transitions in driven-disordered media are another well known non-equili- 
brium critical phenomenon (see e.g. |90 | 1385j ) . These transitions occur at zero tempera- 
ture due to the competition between an external driving force and an internal quenched 
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disorder. Taking thermal fluctuations into account, universal scaling functions become 
available in numerical simulations [3861 13871 13881 1237j and can be compared to exper- 
imental data (see e.g. j389j ). Further numerical, theoretical as well as experimental 
research is needed to provide a complete picture of the scaling behavior of depinning 
transitions. 

Eventually, another open point and perhaps the most important one is the formu- 
lation of a universality hypothesis of non-equilibrium phase transitions. Due to a lack 
of a unifying theoretical framework, a classification scheme similar to Kadanoff 's hy- 
pothesis in equilibrium is still unknown. Thus, the concept of universality remains the 
major tool to order the great variety of critical phenomena. Therefore, we hope that 
the presented picture gallery of universal scaling functions will be useful for future re- 
search, because universal scaling functions are the most impressive and most beautiful 
manifestations of universality at all. 
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A.l Branching process 

In probability theory (see e.g. |39L)ll39l] ^. a branching process models the evolution of a 
population. Each individual in generation n (e.g. particle at time t) creates in the next 
generation (time t+1) k offsprings (particles) with probability p^. It is assumed that the 
probabilities Pk do not vary from individual to individual. A central question is whether 
the branching process leads to ultimate extinction, i.e., no individuals exist after some 
finite number of generations. The expected number of offsprings of each generation 
is denoted by /x. Starting in generation zero with one individual, the expected size 
of generation n is given by /i". Thus three cases are distinguished: the subcritical 
case < 1, the critical case /i = 1, and the supercritical case > 1. Ultimate 
extinction is certain for /x < 1 whereas the probability of ultimate extinction is less than 
one for /x > 1. At criticality the survival probability P{n), i.e., the probability that 
the branching process is still alive after n generations is known to scale as P{n) 6c n~^. 
Furthermore, the probability P{s) that a critical branching process creates exactly s 
individuals scales for large s as P{s) cxs~^/^ j392j . 

Considering activity spreading simulations, the mean field scaling behavior at crit- 
icality equals that of a critical branching process. Thus the mean field values of the 
spreading exponents are given by 9mf = 0, 5mf = 1, as well = 3/2. 



A. 2 Scaling laws 

In general, absorbing phase transitions are characterized by four independent expo- 
nents, for example /3, /?', z^x, and i/y (see e.g. 7 J. All other exponents can be expressed 
in terms of these exponents. The spreading exponents 6, 9, and the dynamical expo- 
nent z are given by [see Eq. ()3.126|) and Eq. 1)3. 128(1 ] 

f|| Z Z/|| f|| Uj^ 

The steady state critical exponents 7, 7' and a are related to (/?, /?', u^, i^n) via 

7 = cj-/?, CT = D i/x + i^ii - /?' , i = Dv^ - 2(3. (A.2) 
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Furthermore, the fractal dimension of the spreading clusters [see Eq. (|2.47|) ] is given by 

Df = z{d + 8) = D-z6 = D- — . (A.3) 

Simple dimensional analysis offers a convenient way to derive these scaling laws. There- 
fore, we remind that the various quantities of interest enter the scaling forms in com- 
binations with their scaling powers, e.g., 

X6p, Xfp^, X'^h, A-^Xa, X-^'Ap^, X-^^x, X-^^L, 

A-^^i?, X-'^H, A^'^^-^iiVo, A-^'^iiiV,, A^^iiP., .... (A.4) 

In the language of real-space renormalization, the rescaling is usually related to the 
transformation x i— > bx, corresponding to A = b^^'^^. Enforcing scale invariance, the 
definition of e.g. the susceptibility leads to 

Xa = |^ =^ I = a-I5. (A.5) 

Similar, the hyperscaling law of the fluctuation exponent is obtained from 

Ap. = {{pD - {p^f) ^ V = Du, - 2(3 . (A.6) 

Taking into account that an initial homogenous particle density p^ ^ may be represented 
by the external field we find 

h = /Oa,o S{t) =^ cr = Du^ + I'll - . (A.7) 

Considering activity spreading, the average number of active sites (usually averaged 
over all runs) scales as 

A^, = const P° =^ 61 = Dv^ - (3 - u^^5 . (A.8) 
Eventually, we mention that the correlation function of active sites 

r(z:„rj = ( (p.(rj - (p.(rj)) (/,.(rj - (p.(r^,))) ) (A.9) 
is expected to decay for translationary invariant systems {l = Li~ Lj) as 

r(r,0,0) « r-^+2-r,x^ (A.IO) 

The correlation function is related to the order parameter fluctuations 

l\p^{dp,h) = J2 nL,Sp,h), (A.ll) 

r 

yielding the scaling power of the correlation function 2/3. Thus, the correlation function 
exponent rj^ obeys the scaling law 

(2 - 7]^) = Dv^ - 2(3 = i (A.12) 

which corresponds to the Fisher scaling law of equilibrium phase transitions. 
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A. 3 Critical exponents, universal amplitude combinations 
and critical parameters 

A. 3.1 Directed percolation 

Table 1: The critical exponents and various universal amplitude combinations of di- 
rected percolation for various dimensions D. In D = 1, the exponents 7, i^x, and z/|| 
are obtained from a series expansion by Jensen |132j . For D = 2 and D = 3 activity 
spreading simulations are performed yielding 6, 9, as well as z jl51lll49] . Additionally, 
the exponent i^y is determined jl93| I149j in order to estimate the full set of exponents 
via scaling laws. 





7^ 1 H 0*M 
U = I llozj 


u — I jioil liyoj 


7~> Q |-j /jAI 
U = jl4ijj 


Mean field 


p — p 


U.z /d4o0^o ) 


U.0oo4 zL U.UUoU 


u.oio =t u.uuy 


i 




1.096854(4) 


0.7333 ± 0.0075 


0.584 ± 0.005 


1/2 




1.733847(6) 


1.2950 ±0.0060 


1.110 ±0.010 


1 


a 


2.554216(13) 


2.1782 ±0.0171 


2.049 ± 0.026 


2 


i 


0.543882(16) 


0.2998 ±0.0162 


0.126 ±0.023 





7 


2.277730(5) 


1.5948 ±0.0184 


1.237 ±0.023 


1 




1.504144(19) 


1.5912 ±0.0148 


1.783 ±0.016 


2 


5 = a 


0.159464(6) 


0.4505 ±0.0010 


0.732 ± 0.004 


1 


e 


0.313686(8) 


0.2295 ±0.0010 


0.114 ±0.004 





z 


1.580745(10) 


1.7660 ±0.0016 


1.901 ±0.005 


2 


^(1,0) [nu 


1.46 ±0.12 


1.65 ±0.09 


1.83 ±0.11 


2 


x(+i,o) ™ 

X(-1,0) >^ 


0.033 ± 0.004 


0.25 ±0.01 


0.65 ±0.03 


1 


Rx ma 


0.60 ±0.04 


0.72 ± 0.04 


0.86 ± 0.08 


1 


u [2SH1 


0.833 ±0.011 


0.704 ± 0.013 


0.61 ±0.02 


1/2 



A. 3. 2 Compact directed percolation 

Table 2: The critical exponents of the compact directed percolation universality class 
{D = 1). The critical behavior equals that of the 1 ± 1-dimensional voter model jl77j 
and is tractable analytically due to the mapping to random walks [1821 [T9T]. 



(3' 


(3 






z 


5 


6 


1 





1 


2 


2 


1/2 
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A. 3. 3 Manna universality class 

Table 3: The critical exponents and various universal amplitude combinations of the 
Manna universality class below the upper critical dimension = 4. The data of the 
exponents P, a, Uj^, and 7' are obtained from steady state simulations jl39lll52j . whereas 
activity spreading reveals the values of a, 5, 6, and z |152j . The exponents P' and are 
determined via scaling laws. In particular, the values of fy are in good agreement with 
those of direct measurements of the order parameter persistence distribution |379j . In 
case of the one-dimensional models we observe a splitting of the universality class. 





D = l 


D = 2 


L> = 3 


Mean field 




0.382 ± 0.019 


0.639 ± 0.009 


0.840 ±0.012 


1 


/?' 


0.319 ± 0.052Manna 


0.624 ±0.029 


0.827 ±0.034 


1 




1.347 ± 0.09lManna 

1.760 ± O.OeOcTTP 


0.799 ±0.014 


0.593 ±0.013 


1/2 




1.876 ±0.135Manna 

2.452 ± 0.106cTTP 


1.225 ±0.029 


1.081 ±0.027 


1 


a 


2.710 ± 0.040Ma„na 

1.770 ± 0.058cTTP 


2.229 ± 0.032 


2.069 ± 0.043 


2 


i 


0.550 ± 0.040Manna 

0.670 ± 0.040cTTP 


0.367 ±0.019 


0.152 ±0.017 





7 


2.328 ± 0.044Manna 

1.388 ± 0.040cTTP 


1.590 ±0.033 


1.229 ±0.045 


1 


v± 


1.592 ± 0.040Manna 

1.619 ± 0.028cTTP 


1.541 ±0.025 


1.744 ±0.029 


2 


a 


0.141 ±0.024 


0.419 ±0.015 


0.745 ± 0.017 


1 


5 


0.170 ±0.025 


0.510 ±0.020 


0.765 ± 0.025 


1 


e 


0.350 ± 0.030 


0.310 ±0.030 


0.140 ±0.030 





z 


1.393 ± 0.037 


1.533 ± 0.024 


1.823 ±0.023 


2 


^(1,0) 




1.81 ±0.03 


1.91 ±0.08 


2 


X(+1,0) 
X(-1,0) 




0.41 ±0.05 


0.8 ± 0.05 


1 
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A. 3. 4 Critical parameters 



Table 4: Estimates of the critical probabilities {p^ and A^.) for the directed percolation 
universality class. 





DPbcc 


CPsc 


PCP,, 


D 


= 1 


0.70548515(20) [M] 


3.297848(22) [M] 


0.077093(3) [nn] 


D 


= 2 


0.34457(1) dnS] 


1.64877(3) |M] 


0.20053(9) [US] 


D 


= 3 


0.160958(6) [Sni 


1.31686(1) innsi 


0.25803(13) [US] 


D 


= 4 


0.0755850(3) [211] 


1.19505(15) [US] 




D 


= 5 


0.0359725(2)W 


1.13846(11) [US] 


0.29874(15)111] 



Table 5: Estimates of the critical densities for the Manna universality class. 





Mannasc 


CLGgc 


CTTP,, 


D 


= 1 


0.89199(5)1111] 


1/2 


0.96929(3) [ISg 


D 


= 2 


0.68333(3) [HI] 


0.34494(3) [ISIj 


0.69392(1) [ISg 


D 


= 3 


0.60018(4) [HI] 


0.21791 (9) [HI] 


0.60489(2) [HI] 


D 


= 4 


0.56451(7)1111] 


0. 15705(10) [m] 


0.56705(3) [HI] 


D 


= 5 


0.54704(9) [HI] 


0.12298(15) [HI] 


0.54864(5) [HI] 
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A. 3. 5 Crossover parameters 

Table 6: The range of interactions R, the corresponding number of interacting neigh- 
bors on a two-dimensional square lattice and the effective range of interactions i^^tf 
for which simulations have been performed. Additionally, the values of the critical 
densities arc listed. 



R 






^Manna 


„CLG 


„CTTP 


1 


4 


1 


0.68333(3) 


0.34494(3) 


0.69392(1) 


2 


12 


7 

3 


0.60595(4) 


0.22432(4) 


0.63649(2) 


4 


48 


8 


0.54378(6) 


0.16802(7) 


0.55005(3) 


8 


196 


1546 
49 


0.51549(8) 


0.14050(9) 


0.51688(4) 


16 


796 


25274 
199 


0.50481(9) 


0.12977(10) 


0.50552(6) 


32 


3208 


204875 
401 


0.50146(10) 


0.12598(11) 


0.50161(7) 


64 


12852 


13146247 
6426 


0.50045(12) 


0.12499(16) 


0.50046(8) 


128 


51432 


105255421 
12858 




0.12465(19) 


0.50019(9) 



Table 7: The range of interactions i?, the corresponding number of interacting neigh- 
bors on a three-dimensional square lattice and the effective range of interactions R^s 
for which simulations have been performed. Additionally, the values of the critical den- 
sities are listed. 



R 






^CTTP 


1 


6 


1 


0.60483(1) 


V2 


18 


5 
3 


0.56770(2) 


2 


32 


39 
16 


0.54406(3) 


2^/2 


92 


219 
46 


0.51720(4) 


4 


256 


597 
64 


0.50694(6) 


4^2 


750 


2383 
125 


0.50260(7) 


8 


2108 


20016 
527 


0.50136(8) 


16 


17076 


217886 
1423 


0.50044(8) 
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